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MAGNETIC SUSCEPTIBILITIES OF IRON GROUP SALTS 
AT LOW TEMPERATURES! 


By A. F. JoHNson? AND H. GrAyson-SmITH® 


Abstract 


A magnetic balance of the Sucksmith ring type has been constructed. It is 
capable of a precision of 0.5 to 1% in measurements of displacement. The esti- 
mated maximum error in susceptibility measurements is 2%. This apparatus 
has been used to measure the temperature variation of the susceptibilities of 
chromium sulphate, hydrated and anhydrous, and of the nitrates of chromium, 
cobalt, and nickel. Measurements were made at several values of the magnetic 
field strength, so that a correction for small amounts of ferromagnetic impurity 
could be obtained. The reciprocals of the corrected specific susceptibilities were 
plotted against the absolute temperature. The intercepts on the temperature 
axis of the resulting straight lines gave the values of the Curie-Weiss constant. 
These values have been compared with those for other salts of the same ions. 
The effective Bohr magneton numbers were calculated and compared with theo- 
retical values given by van Vleck. 


Introduction 


The salts of the transition metals, chromium to nickel, are all paramagnetic 
owing to the metallic ions, which have incomplete d-shells, and so have mag- 
netic moments. According to the simple Langevin—Debye theory (1, 18), the 
susceptibility of a paramagnetic material should depend upon temperature 
according to Curie’s law, 


“Ve OF Es (1) 


However, for most paramagnetic solids, a better description of the experi- 
mental behavior is given by the Curie-Weiss law 


x = C/(T — 8), (2) 
where @ is an empirical constant, positive or negative. 


The Curie constant, C in Equations (1) or (2), is determined by the magnetic 
moment of the ion. It is usual to define the effective magnetic moment per 
mole, ness, and the effective magneton number, p.;;, by the relations 


tery = V3BRC = V3Rxm(T — ©) = Persp, (3) 


where xm is the molar susceptibility, and ug = 5565 c.g.s. units per mole is 
the Bohr magneton. The effective magnetic moment is the value which would 


1 Manuscript received October 15, 1949. 
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be required according to classical theory in order to produce the observed 
susceptibility. It can be calculated by quantum methods (7, 18) for any free 
ion of known electronic configuration, but may be modified in various ways for 
bound ions in a crystal. These modifications will be discussed later in con- 
nection with the empirical values of ueyz determined by the authors. 


As far as the magnetic moment is concerned the theory gives a satisfactory 
explanation of the facts. Theoretical explanation of the Weiss constant, 0 
in Equation (2), rests upon much less certain ground. Equation (2) also holds 
for ferromagnetic materials in the high temperature paramagnetic state, with 
98 approximately equal to the Curie temperature. According to Heisenberg’s 
theory of ferromagnetism (1, 6) the behavior in this case can be explained by 
assuming a strong exchange interaction bétween parallel aligned electron spins. 
A similar effect, involving interactions between neighboring ions, could also 
explain Equation (2) in cases where @ is positive and of the order of 0 to 50°. 
This is almost certainly the case for certain of the transition group chlorides, 
which are known to become weakly ferromagnetic at very low temperatures 
(5, 12). On the other hand, some of the manganese compounds are known to 
become antiferromagnetic (2,14), probably owing to stability of the anti- 
parallel alignment of magnetic moments. The internal electric field of the 
crystal can also cause a tendency for the ions to align themselves, and so cause 
a temperature variation of x which simulates the Weiss law over a fairly wide 
range of temperature (1). 


In any case, one would expect a progressive change in the values of 0 for dif- 
ferent salts of the same ion, with different inter-ion distances and different 
anions. It is therefore interesting to compare the paramagnetic constants of 
as many different salts of the transition group elements as possible. Such 
comparisons have been given by Van Vleck (18) and by Burton et al. (3), but 
the number of salts hitherto studied is not sufficient to make the evidence 
conclusive. 


The principal purpose of the present investigation was to add to the experi- 
mental evidence on this point by studying several more simple salts of the 
ions in question. It was hoped that this might make it possible to decide 
whether or not the interactions between ions are the dominant factor. This 
hope has been only partially justified. However, the experiments have given 
further evidence, both on this point and on factors which modify the effective 
magnetic moment. 


Apparatus 


The magnetic balance constructed for this investigation was based upon the 
design of Sucksmith (16, 17) as adapted for low temperature measurements by 
Jackson (8). The essential components are shown in Fig. 1. The powdered 
sample was contained in a hard gelatin capsule C, which was sealed with 
glyptal laquer and then cemented to a quartz rod CF,F:2, connected to a brass 
rod at F,. The brass rod was attached by means of a coupling J to a phosphor- 
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bronze ring R, which carried two small mirrors M,M:. The ring was adjusted 
so as to be nearly circular under the weight of the suspension and sample. 
Then, when the magnetic field was applied, the force on the sample due to 
the field deformed the ring. The mirrors were attached at the points where 


Te 


Detail at J’ 


Brass Rod 


0.75 mm. 
Sheet Brass 





Detail at "C” 


,: 


Fig. 1. Sucksmith type magnetic balance. 


~” 


Nw 


the angular deflection due to this deformation was a maximum. Light from 
a distant source passed through the lens L, was reflected by M, and M, and 
was focused by a microscope objective. By this device the displacement of 
the sample was magnified approximately 100 times. 
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The position of the suspended system could be adjusted by means of the 
screw adjustments and the sylphon bellows A. If necessary, the shape of the 
ring could be adjusted by adding small weights w, which could be manipulated 
by means of a hook passing through a sylphon bellows at the back of the case. 
The front of the case was sealed by means of a thick plate glass window. 


The rings were made out of strips cut from phosphor-bronze sheet, bent into 
shape and soft-soldered together at the point of support. Several rings were 
made, with diameters 5 to 6 cm., width 2 to 3 mm., and thicknesses 0.05 to 
0.125 mm. Each ring was shaped by trial to an ellipse of such eccentricity 
that it was pulled approximately circular by the load of the suspension and 
specimen. The ring finally selected had a sensitivity of 5.018 cm. per gm. ata 
scale distance of 65 cm. The deflection was linear for weights up to 500 mgm. 
By means of a vernier the vertical travel of the microscope could be read to 
0.01 mm. Over the range of displacements observed in the experiment this 
gave a precision of 0.5 to 1% in the measured deflections, and an estimated 
maximum error in the susceptibility measurements of 2%. 


The phosphor-bronze ring with its attached mirrors was very susceptible to 
vibration, and some form of damping was necessary. To isolate the system 
from external vibrations, the cryostat with the balance attached, the light 
source, and the microscope were all mounted on a marble slab E, resting on 
two strips of sponge rubber. To damp internal vibrations, the rod by which 
the specimen was suspended carried a hollow monel metal ring r which was 
immersed in an annular oil bath D. A flat spiral spring S of thin phosphor- 
bronze prevented the specimen support from swaying, or from being pulled 
sideways towards the pole-pieces of the magnet. 


The specimen hung inside a long glass tube in a specially shaped Dewar 
flask. Various refrigerants were introduced into the Dewar flask through a 
double-walled glass tube at 7), and the gas from the evaporating liquid was 
pumped off through 72. The thin metal baffle B reduced radiation down the 
tube. By means of the ground bronze joint G, the surrounding tube and the 
Dewar flask could be removed to change the specimen. After the specimen 
was mounted, the tube in which it hung and the case containing the balance 
were filled with nitrogen at atmospheric pressure. Temperatures were indicated 
roughly by the vapor pressure of the cryostat liquid, but it was found neces- 
sary to measure them more accurately by means of a platinum resistance 
thermometer, which was wound around the lower part of the glass tube and 


cemented in place with glyptal. 


For small specimens, where the magnetic field may be assumed constant 
over the volume of the specimen, the vertical component of the force is given 


by F, = mxH,aH,/ax, (4) 


where m is the mass of the specimen, X is the specific susceptibility, and the 
pole-pieces of the magnet lie in the y-direction. The pole-pieces were 9 cm. 
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in diameter with a gap of 3.9cm. The field between these poles was explored 
carefully, and the specimen was placed in the position where HOH /dx was a 
maximum, and was changing least rapidly. This quantity was then very 
nearly constant over the volume of the specimen, and did not vary appreciably 
over the small displacement which occured during a measurement. It was 
therefore unnecesary to use a restoring force to bring the specimen back into 
the original position in the field. 


In order to calibrate the magnet and to find the absolute values of HdH/dx 
for different current strengths, a series of measurements was made with man- 
ganous sulphate. The susceptibility of this material has been well established 
over the temperature range concerned by Onnes and Oosterhuis (11). It was 
found that the calibration factor varied slightly with the temperature of the 
cryostat. The appropriate factor was therefore used for each of the four 
temperatures at which measurements were subsequently made. 


Correction for Ferromagnetic Impurities 


A serious source of error in the measurement of small susceptibilities is the 
presence of ferromagnetic impurities, either in the suspension system or in the 
sample. To test the suspension system a series of measurements was made 
with an empty gelatine capsule. Using the largest magnetic field which could 
be obtained from the magnet, the deflection was too small to observe, either 
at room temperature or at the temperature of liquid nitrogen. Therefore 
no correction was necessary for the suspension system. 


In order to eliminate ferromagnetism in the samples, measurements were 
made at several different field strengths at each temperature. It can be assumed 
that any ferromagnetic material present would be magnetized to saturation at 
the field strengths used, and would have a constant magnetic moment M. The 
total force on the specimen is then + 


F = MoH /dx + mxHdH/dx = a\H\| + 6\H\?, (5) 


since 0/1/dx is everywhere proportional to the absolute value |H| of the field 
at the center of the magnet. Therefore, if there is a small amount of ferro- 
magnetic impurity, graphs of F against either |H|* or Hd@H/dx will deviate 
slightly from the ideal straight line through the origin (Equation (4)). For 
large field strengths the slope of the curve of F against HOH /dx should approach 





the value mx. 


For the manganous sulphate standard a graph of balance displacement 
against |//|? was linear within experimental error, and no correction was neces- 
sary in using this standard to calibrate the magnet. In other cases a small 
deviation from linearity was observed. This is illustrated in Fig. 2 for the 
case of cobalt nitrate, which is typical of the various salts measured. In this 
figure, balance displacement is plotted against HdH/dx for measurements 
taken at four different temperatures, and the small but distinct deviation can 
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be easily seen. In all cases the curves were very nearly linear for the larger 
field strengths, and corrected values of mx could be readily determined from 


the slopes. 














Displacement (cm.) 























H 9H , 107° 
ox 


Fig. 2. Balance displacement vs. HOH /Ox for cobalt nitrate, at four different temperatures, 
to illustrate the small deviation from proportionality due to ferromagnetic impurity. 


Experimental Results 


New measurements were made on the salts listed in Table I, at room tem- 
perature, in liquid ethylene, liquid oxygen, and liquid nitrogen. The materials 
were obtained from British Drug Houses (Canada) Ltd., and were specified as 
having amounts of impurity less than those given in the list. The anhydrous 
chromium sulphate was prepared from the hydrate by heating it at 400°C. 
until the weight was constant. The loss of weight was 1 to 2% greater than 
that which would be accounted for by the removal of 18 H,O per molecule. 


The susceptibilities of these salts were derived from graphs of F against 
HdH/dx, as described above. They were then corrected for the small dia- 
magnetism of the water of crystallization and of the anions, using values given 
by van Vleck (18). The final results are collected in Table II. For each 
specimen the corrected values of 1/x were plotted against the absolute tem- 
perature (Figs. 3 to 7). In all cases the resulting graphs were linear, showing 


- 


that the Curie-Weiss law (Equation (2)) is accurately obeyed down to 77°K 
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TABLE I 
SPECIFICATIONS OF MATERIALS USED 


Chromium sulphate, Cr2(SO4)3.18 H2O 


CaO os, Sok 2g tos Coe ey a ao ee 0.003% 
RPIMISEO Re os e835 Oa eae. oes Valbaed EDER DARK aN 0.01% 
PRAM Oe oe PR cicy Cah loka Cee a eee 0.5% 


Chromium sulphate anhydrous, Cr2(SO4); 


Chromium nitrate, Cr(NO3)3.9 H:O 


Chloride (Cl)....... Pa So ek teca tone 0.003% 
Sulphate (SO,).... ce ie See eee eee 0.01% 

PROM OGD Toe 8 ss desta neo Gkra ee ae bee a ween 0.01% 

ANG 5. sik Kee ee eS CEE eon Sata ale 0.5% 


Cobalt nitrate, Co(NO3)2.6 H2O 


Chloride (Cl)....... ceed wea py BAe ale hae cunet orale 0.001% 
SERIE CIT wscieiarg ie mv sa acon ob aig airings 5 Sra salle oat 
WeaeclR Oh. ose. ss Peek OM Se PE ee ferret e, 0.003% 
Nickel (Ni)... .. ph eee Roe ate ce 0.002% 
Alkalis and alkaline earths (Na)................-. 0.03% 
MNS cox 5 oes ala SOR aoe wr oe wera «IO DOROCION 


Nickel nitrate, Ni(NO3)2.6 H:O 


rte WE ee aie I a Ae ee tain yea ee 0.002% 
Sulsteite Css), 2... sc seas pest Recess ale aleae 0.01% 
Cobalt (Co)....... Re ENTS aewiltes didia bee Sipe ...---0.0005% 
Fron (Fe)... earl Ste ng wave sii eos on OO 
Alkalis and alkaline earths (Na).................. 0.03% 


TABLE II 
CORRECTED PARAMAGNETIC SUSCEPTIBILITIES 


























Mass of Diamagnetic | Temperature, | Corrected 
Salt sample, correction eK | susceptibility 

gm. (per gm.) | (per gm.) 

Cr2(SO4)318 H2O 0.1171 0.50 X 10-8 293.8 14.8 & 10-* 
168.9 | 26.7 
90.5 53.9 
77.3 67.9 

Cr2(SOx)s 0. 1529 0.31 X 10-8 293.0 26.8 X 10-6 

168.9 | 46.6 
90.3 | 98.1 
Pee {108.0 

Cr(NO3;)3.9 HO 0. 2075 0.43 X 10-° 294.5 13.3 X 10-8 
168.9 | 23.1 
90. 4 | 50.3 
Thom 59.1 

Co(NOs3)2.6 H.O 0. 1393 0.40 « 10-6 294.8 | 31.1 X 10-* 
168.8 54.5 
90.4 109.2 
77.0 127.6 

Ni(NOs3)2.6H2O 0. 1268 0.40 x 10-5 295.8 | 12.7 X 10-* 
169.1 | 23.7 
90.3 49.2 
77.2 61.2 
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Cr2 (SO.),-18 HO 








6=19 Temperature , °K 


Fig. 3. Reciprocal susceptibility, 1/X, for chromium sulphate hydrate, showing determina- 
tion of the Weiss constant, 9. 





50 100 150 200 250 300 


6-1 Temperature , °K. 


Fig. 4. Reciprocal susceptibility, 1/X, for anhydrous chromium sulphate, showing deter- 
mination of the Weiss constant, 
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Fig. 5. Reciprocal susceptibility, 1/x, for chromium nitrate, showing determination of the 
Weiss constant, 9. 
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Fig.6. Reciprocal susceptibility, 1/x, for cobalt nitrate, showing determination of the 
Weiss constant, 9. 
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6= 2! ‘Temperature ,°K 
Fig.7. Reciprocal susceptibility, 1/X, for nickel nitrate, showing determination of the Weiss 
constant, 8. 


Effective Magnetic Numbers 


Experimental values of the effective magneton number, p.;;, were calculated 
for each salt, by means of Equation (3), and are given in the last column of 
Table III. The only previous measurements below room temperature on any 
of these salts are those of Serres (13) on chromium sulphate hydrate, and those 
of Janes (9) on nickel nitrate. In each case, measurements were made at only 
one value of the field strength, and no correction was made for ferromagnetic 
impurity. Serres gave a value of p.yzy = 3.85 for Crt**, which is so much 
smaller than the present value as to suggest some radical difference in the ma- 
terials measured. Janes used a different method of comparing experiment 
with theory, and did not quote values of pers. His values for the susceptibilities 
are somewhat higher than those found by the present authors. Chatillon 
(4) has shown that the values of p.s; for Co** depend on the mode of prepara- 
tion and history of the material. It is possible that variations can occur in the 
arrangement of the water of crystallization molecules, and that these can 
cause variations in the internal field. 


According to the Langevin—Debye theory of paramagnetism (1) the molar 
susceptibility is given by 


Xm = N?g?J(J + 1) up’/3RT, (6) 
which gives eae 
Persp = EV I(J +1), (7) 


where J is the resultant angular momentum quantum number of the magnetic 
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ion, and g is the corresponding Landé factor. This theory is based upon the 
assumption that every ion is in the same quantum state, and it is expected to 
be a good approximation when the energy separation in the low-lying multiplet 
is large (hy >> kT). In the case of narrow multiplets and Russell—-Saunders 
coupling van Vleck (18) has shown that 


bey = VL(L + 1) + 4S(S + D). (8) 


The complicated case of intermediate multiplet widths has also been discussed 
by Laporte (10) and van Vleck (18). _In this case the effective magneton num- 
ber may vary with the temperature. 


In the case of the salts of the iron group neither formula fits the experimental 
results very well. Stoner (15) first suggested that the internal electric field 
due to neighboring atoms in the crystal might quench the contribution of the 
orbital angular momentum to the magnetization. The expression for the 
magneton number then becomes 


Pert = V4S(S + 1). (9) 


In general where quenching is incomplete the value of p.sz should lie between 
those given by Formulas (8) and (9), as is the case for most of the salts of the 
iron group. 


Table III shows the experimental values obtained in this work, compared 
‘ with the different theoretical values. Column 7 gives the theoretical value 
for multiplet width comparable with kT, as calculated by van Vleck (18). The 
experimental values obtained for the Cr*** ion in both the hydrated and an- 
hydrous sulphates are somewhat higher than those previously obtained for 
chromium salts. These values are only slightly less than the theoretical value 
for narrow multiplets, and indicate very little quenching of the orbital moment. 
In the nitrate the value for the chromium ion indicates complete quenching of 
the orbital moment and partial quenching of the spin moment. The result for 
Cot* indicates partial quenching of the orbital moment, while that for Ni*t* 
corresponds to almost complete quenching of the orbital angular momentum. 


TABLE III 
EFFECTIVE BOHR MAGNETON NUMBERS 


cake Theoretical pesf 
Lowest] Expt? to 5 i a eee 





Ton Salt state eff * ] i ; > jf 
hv>>kT | hv<<kT | hv~kT | Spin only 
del acctnagaecs ‘oh abe et at | aig, 5a 
| Cro(SO.)s 4 F 5/2 1.89 |; 5.20 Ww 1. Sat. ae 
| criNO.)s9 H;O 3.43 | | 
Cot+ | Co(NOs)2.6 H20 | ‘Fo/2 4.58 5.20 6.64 | 6.56 3.87 


Nit+ | Ni(NOs)2.6 H.0 | *F, 2.86 | 4.47 | 5.50 | 5.56 | 2.88 
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The Weiss Constant for Different Salts 


In Table IV, values of the Weiss constant 6 have been collected for as many 
as possible of the simple salts of the iron group. Values marked with an asterisk 
are from the present investigation, and the remainder are from a comparison 
table given by Burton et al. (3). 

TABLE IV 
CuRIE-WEISS CONSTANT FOR IRON GROUP SALTS 


| 
lon | Salt | © | 


— — —— — — —— = 
| sat | 0 | Salt | @ | Sit =| © 


i 
— | ————— | | | 
Cr+++| CrCls peer. 11 | Cr2(SO4)s.18 H,O* | 19 | Cr(NO,)3.9 H:0*| 20 








Mn++! MnCl, | 0 | MnSO, 24 | MnSO..4H:0 | 0| 
Fet+ | FeCl. | 20 | FeSO, le | FeSO,.7 H:O ee | 
| | 
Cot+ | CoCle | 20} CoSO, |—45 | CoSO,.7 H:0 |-14 | Co(NO;)2.6 H:0*) 8 
Nit+ | NiCl, | 67 | NiSO, |—79 | NiSO,.7 H:0 | 59 | Ni(NOs)2.6 H.0*| 21 
Conclusion 


As indicated in the introduction, it was hoped that comparison between 
different simple salts would show at least a qualitative correlation between 
the complexity of the salt and the value of the Weiss constant, 0. This hope 
has been only partially realized. 


The cases of hydrated and anhydrous chromium sulphate may be significant. 
These salts were found to have nearly the same effective magneton numbers, 
showing nearly the same degree of quenching of orbital momentum by the 
crystalline field. Then, if the necessity for using the Curie-Weiss law (Equa- 
tion (2)) is due to interaction between ions, one would expect the value of 0 
for the hydrated salt to be less than that for the anhydrous. On the other 
hand the presence of the polar H.O groups in the hydrated salt makes the in- 
ternal electric field larger than in the anhydrous salt. The fact that the ob- 
served value of 6 was larger for the hydrate suggests that the crystal field is 
the predominating factor in chromium sulphate. Measurements on these 
two salts in particular should be extended to still lower temperatures, where 
deviations from the Curie-Weiss must appear. 


The only other correlation noticeable in Table IV is that the large values 
of 9 for nickel salts, considerably larger than the corresponding cobalt salts, 
extend to the newly measured nitrates. 
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THE NEUTRON DENSITY NEAR A PLANE SURFACE IN A 
CAPTURING MEDIUM! 


By JEANNE LECAINE? 


Abstract 


The variational method is applied to the Milne problem with capture, to 
obtain, for the density and for the angular distribution, expressions which give 
quite accurate values and which have a relatively simple analytical form. Tables 
of the density and angular distribution are included. The extrapolation distances 
determined by the Variational Method agree to high accuracy with the exact 
values. The angular distribution determined by the Variational Method is 
shown to be accurate within 0.05% in all cases where exact results are available. 
Tables are also given for the density and angular distribution in the Milne pro- 
blem with capture and constant production. 


Introduction 


The formal solution of the Milne problem with capture has been obtained 
by Adler and Mark (1) using the Wiener—Hopf technique. Starting from the 
integral equation 


f(z) = 





=| Fe) E,(\z — 2'|) de’ (1) 


they have shown that the asymptotic solution for the density, expressed as a 
function of distance from the boundary, is 


— 


a(v? — a) 





fas(2 sinh v (z + 20(a)), (2) 


where v(a) is the rea! root of the transcendental equation 





Pn (3) 


v 


and Zo(a) is given by 


_ l= = ie arth (vp) a (4) 
tiie 3 [1-=«) warthaf +[70 59] 


This integral has been evaluated by numerical integration for a = 0, 0.1, 0.2, 
0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9. The values of zo(a) and »(a) for these values of a 
are tabulated in the paper by Adler and Mark. As a— 1, 20(a) becomes infinite, 
and, for a near 1, 20(a) is given by 








Zo(a) = —+5 + —log (1 — a) — 5 lo 2 a * log 2+. (5) 


ie 


1 Manuscript received October 27, 1949 
Contribution of the Chalk River Laboratories of the National Research Council of Canada. 
Issued as N.R.C. No. 2180. This work was completed April 17, 1945. 
2 Now at Oklahoma Agricultural & Mechanical College, Stillwater, Oklahoma. 
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Adler and Mark and Davison (2) have also shown that the density at the 
boundary is nme 
. v- ° a 
fo) = 4/ © - jo (6) 


They have obtained the angular distribution at the boundary, and the devia- 
tion of the true density from the asymptotic density, as rather complicated 
integrals. A few isolated values of the angular distribution have been obtained 
by numerical integration. 

In the present paper the application of the Variational Method to this 
problem is discussed. It is intended to complement the rigorous work done 
and to give a simpler way of obtaining quite accurate values of the density 
and angular distribution. 


The application of variational technique to the Milne problem without 
capture has already been discussed (5). It has been shown that a good approxi- 
mation to Zo can be obtained with a crude trial function, and that by choosing 
a trial function the behavior of which is similar to that of the true function 
it is possible to obtain an expression for the density which has a much simpler 
analytical form than that of the true density and which is a very good approxi- 
mation to the true density over the whole range.* 


In the following, a similar discussion is carried out for the Milne problem 
with capture. The trial function used is the logical generalization of that used 
‘in the case a = 0. Fora > 0 the true density has not been calculated but from 
physical considerations the accuracy of the approximation would be expected 
to improve as a increases, since, as a increases, the distinction between the two 
semi-infinite media becomes less marked. In Section 6, Tables VII, VIII, some 
figures are given indicating the accuracy of approximation to the angular 
distribution, for a = 0.1, 0.2, 0.4, 0.6. In every case where the true value was 
determined it was found to be within 0.05% of the approximate value. cs 


Notation 


The variational technique used is based on the fact that (3, 8) for the non- 
homogeneous integral equation 


q(z) = f q(z')K(z,2')dz’ + h(z), (7) 
q(z) being the unknown function, and K (z, 2’) a symmetric kernel, the func- 
tome! F(q) = fa(z) [q(z) —fq (2) K (z, 2") dz’|dz (8) 
[fh(z)q(z)ds}? 
is a minimum when q(z) is the solution of (7). F(q) is then equal to 
1 
fh@)q(2)dz 


* A comparison of the true density and the approximation to it obtained by Variational Technique 
has been given by C. Mark (7). 
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We introduce the notation, 


wee 


Gav, 2) = 1" {e*B, [(1 + ve] + e-" Ex [(1 — »)2] — 2Z,(@)} 9) 


© e 
=(1-— ar] — du 


1 u(u? — pv?) 











, 1 — pal , ete 
Gov, z) = ae {e"E,[(1 a v)z|—e *F,[(1 _ v)z]} (10) 
v 
=—(1 -«)| aie, 
; eo} 


The function G;(y, z) is finite, positive and monotone decreasing in g; 


G;(v, 0) = ae 
: 


and asz—>o, a a 
. —alyve ~ 
G3(y, 2) ~ - ————- —— + 

l—, gs 


The function G2(v, z) is finite, negative and monotone increasing in 3; Go(v, 0) 
— (1 — a)e~* 


The slope of G2(v, z), asz—- 0, 
(1 — v*)z 


= — landasz—~,G2(v,2) ~ 
becomes infinite logarithmically. Also, 


: Bek se ay 
lim —G;(v, z) = E3(z) 
v>o vp 





(11) 
lim Go(v, z) = E2(z). 
v>od0 
We note also that 
<< | e"Ex(\z — 2'|)d2’ = e+ we ans (12) 
and 
= =| e "E,(\z — 2’|)d2’ =e” — Gul, 9) — Siler) — Gal, 8) . (13) 
2 0 2 


To obtain an approximation to f(z) we can proceed by two slightly different 
methods. We can write 


f(z) = Ky [sinh v(z + 20) + qi(z)], (14) 
9,2 | 

where qi (z) is O(e~*) asz— ©, k > vand K, = 4/2 if j(0) = 1. 
a(v? — a) 


Substituting this in (1) and using (12) and (13) we obtain for qi(z) the nonho- 
mogeneous equation 


sinh vz 


G3(», 2) iL a G2(r, 2), (15) 


cosh vZo 
2 





l1-—af~ ,, fa pea 
qi(z) = —> | a \Ex(|\z — 2’|)ds’ + 
0 








oor 
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from which we can obtain by variational technique an approximation to 29 
and qi(z). 


Or we can write f(z) = Kz [e* + q2(z)] , (16) 


vSo 
where K2= KS , and go(z) =e" e”* + 2e~”” qu(z). (17) 
Substituting (16) in (1) and using (12), we obtain the equation 


l-a ¥ : G3(v,z) + Ga(v, 2 
qe(z) = —— | qo(2’)Ex(|2 — 2’|)d2’ + Ga(y, 8) Ca, 5) ‘ 
2 J 0 2 
from which an approximation to 2» and g2(z) can be obtained. 


(18) 


Method I 
Let us consider first Equation (15). Here the nonhomogeneous term is 


cosh vZ0 sinh VZ0 


G3 (v, ) + ——— Gal, 2). (19) 


2 





h(z) = 
h(z)qi(z)dz has been evaluated by Marshak as follows. Multiply (15) by 


0 
sinh v(z + zo), and integrate with respect to z, over the range (0, ~). 


co? ae Pe © 
| sinh p(s + 20)qi(z)dz = +— - | sinh v(z + Zo) [ atone (ls — 2’|) dz’dz 


0 = ~0 0 





ro rc a 
+ | sinh v(z + 2) E Sc 3(v,3) + 2) += = Go(y, 2) fas. 


0 


By interchanging the order of integration and using (12) and (13), this becomes 


rs. = vZo sinh vZ0 
qi (2) G3 (v, 2) a go RH) dz 
0 





= — (ss v(z + 20) l° os s G3(v,z) + m = Ge (yr, 2) dz 
1—- by 2y? 
= — —_— *| (cosh 2vzo) log (1 — v?) + log (1 — »*) + - “|. (20) 
8v? 1-?r 


An appropriate choice of trial function for qi(z) is 
qi(z) = AG;(v, 2) + BGo(», 2) (21) 


The function (21) is suitable since it has the correct asymptotic behavior as 
z— o, and since its slope becomes infinite logarithmically as z — 0, which is 
easily seen to be true of the true function qi(z). 


With 9:(z) given by (21) the functional (8) becomes 


K, + Kia + x,(2Y 


(«+08) 
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where Ky, Ke, K3, Li, Le are integrals involving G2, G;, E1 which can be eval- 


. a 
uated analytically.* Only the ratio = can be determined from the extremum 





property of F(g:). We have a = 0 
(4) 
ss B _ 2LaKi — LK, (23) 
A 211 K3 — L2Ke 
A good approximation to 2% can be obtained by setting F(q:) "aes 
2)q(z)dz 


where ; for F(q:) is given by (23) and fh(z)q(z)dz is given by (20), so that 
we have the transcendental equation for 20: 


4KiK;— K} —1 . (24) 
(L?Ks + L2K1 — LiL:K3] » [cosh ad, 3 = | as) 
er. vy 














Thus we obtain an approximation to f(z): 
f%® = K,| sinh v(z+ 20) +A (G0, 2) +5 G2(r, :))| (25) 


in which 29 and “ are known, and A is still to be determined. In order to 
determine A we require some additional property of f(z). For example, we 
can use (6), f(0) = j(0) 4/=. which implies, from (2), the choice of normal- 
a 
ization 
Ki = j(o) 5 = (1 — +’) 
a(v?— i 


So that, making f(z) correct at the boundary we get the equation 








ee ‘ Baw 26 
2(1—r*) {Gr(r.0) +2 Gatoo)} ~ 
and ye — sinh vz 

A selec MR eas kali (27) 


ie B 
log } aia) ae es 
v g ( . A 





Method II 


Method I has the disadvantage that some additional information regarding 
f(z), e.g. f(0), is necessary before f(z) can be determined. If we consider 
* The expressions for Ki, Ko, Ks, Li, L2 are given in the Appendix. 





APSARA TNS nr 8 
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Equation (18), an approximation to f(z) can be completely determined from 
the extremum property of F(q). For Equation (18) we have 


Gs (v8) + Gals), 


h 
(z) = ‘ 
The value of | h(z)q2(z)dz, in terms of 2, can be obtained from (20), by 
0 


making use of the relation (17) between g:(z) and g2(z) and also of the equation 


[a (2) a oa = [eTicee 8 + Galea dz 





which can be obtained from Equation (18) by multiplying by e~” and inte- 
grating with respect to z over the range, (0, ~). We have 
| an(s) 22%) + Go(v, 2) 7 __ — (l~a) 
0 2 2v 
log (1 — v?) ( 2 2 ) _2vs | 
fa gee Sicieapivie 4 aie e an a8 
[ 2 - 1—rv Il-a a 
To obtain an approximation to f(z), we start from the trial function 
G2(z) = k[e~* + AG;(v, z) + BG2(v,2)], (29) 


which corresponds to (21). If k = e~*””, e’* + Ge(z) has the correct asymptotic 
behavior. In this case F(g2) is a rather lengthy expression which will not be 
given here*(6). By minimizing it, one can determine A and B. 


The value of zo given by this approximation is found from 


—1 











F , 
(9): tales a cn (2 2 2 )] 
iced that aatiiatnenitne a canta 
2 2v l1-—-v l-a 
that is, 
i au5e log (1 — »*) — - : 
ye he eee (30) 
2v 














i=8( e  %. -) 
v I1-rv Il—-a 


No numerical work has been attempted with the trial function (29). For the 
simpler trial function g(z) = e~”, one obtains an analytical expression for 20: 








1 1 2 1 
— —log (1 — logs (1 - »*)| 
1 [4 “tog (1 +») | ze" a) 
Zo = — — log 


2y 1 1 1 a 
v 3 - pK) eet - 


As will be seen below, this expression represents 2 within 0.3%. 








* The present paper is based on Montreal Report M.T. 119, which contains the full expression 
for F(@:2). Anyone interested should write to the Theoretical Physics Branch, Chalk River Labor- 
atories, National Research Council, Chalk River, Ont. 
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Angular Distribution of Outgoing Neutrons 
The angular distribution is related to f(z) by 


hid tee | f(2')e* "dz, (32) 
2u 0 


where yu is the cosine of the angle between the direction of the neutron and 
the outward normal. Thus the angular distribution corresponding to 











fe) = A/ 2 U—”) sinh » (¢ + 20) + AGs(v, 2) + BGs (y, 2)] 
a(v- — a) 
is Ne ees 
, t= a fai — A) “35 en 
in) = 128 4/20=*)| _ _€ 
- + a(’?—a) Ll—w 1l+p 
(1 — «)A Fo log (1 — v*?) — ll + 2p’? log(1 ++) 
(1 — pv*)v | ‘ l—-« B/, 
- Ake he log (1 — v*) pv + Quy log(1 + +) 
(1 — p?y?)p { l—s bb 
and (33) 
l—«* clio | 2 o 
vn(u)du = —— = ‘| — e log (1 — v) — e~”* log (1 + v) 
l, 4 V aly —-a) 
A(l—a)f | 2p z (- weGey} 
——__— + — 4v log 2 + —— log (1 — — 2) +———__ [| — 
eo ee ee SF 1+» 
peasy. (1 + v) log (1 — |. 
v 


Numerical Results 


It has already been noted that the variational method gives a better approxi- 
mation to zo than to the density function. Table I gives a comparison of the 











TABLE I 
a | From (81) | From (24) | By Integration 
0 | 0.7083 0.7104 | 0.7104 
0.1 | 0.7095 | 0.7106 0.7106 
0.2 | 0.7105 | 0.7113 0.7113 
0.3 | 0.7121 0.7127 0.7127 
0.4 | 0.7150 0.7154 0.7154 
0.5 | 0.7202 0.7204 0.721 
0.6 0.7298 0.7300 0.730 
0.7 0.7484 0.7484 0.749 
0.8 | 0.7848 0.7847 0.785 
0.9 | 0.8538 0.8539 0.85 
1 | 1.0000 1.0000 1.0 

















A nn 
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values for (1 — a) 2 obtained from the simple expression (31) and the more 
accurate approximation (24), with the values obtained from numerical inte- 
gration of the exact expression (4). 


The simple expression (31) is thus seen to give a good approximation to 2. 
Its maximum error (0.3%) occurs for a = 0 and the accuracy increases with 
increasing a. Equation (24) is of course a much better approximation. For 
a > 5, when the numerical integration is difficult, (24) and even (31) seems 
to be more accurate than the result of numerical integration. 


Method I has been applied to determine f(z) and ¢(u) for a = 0.1, 0.2, 0.3, 
0.4, 0.5, 0.6, 0.7, 0.8, 0.9. We get 
f(z) = sinh v(z + 20) + AG;(v, 2) + BG2(y, 2) 
= sith xis +404 - = {(A + Boe*Ex((1 + »)2) 
+(A — B)e"E,((1 — v)z) — 2AEx(z)}, (84) 


where the constants v, zo, A, B, depend on a, and are given in Table II. 

















TABLE II 

a | v | Zo | B+A | A | B 

0.1 | 0.5254 | 0.7896 | 0.5794 | 0.2581 | 0.1496 
0.2 | 0.7104 | 0.8891 | 0.7883 | 0.3052 | 0.2406 
0.3 | 0.8286 | 1.0181 | 0.9219 | 0.3748 | 0.3456 
0.4 | 0.9073 | 1.1923 | 1.0070 | 0.4868 | 0.4902 
0.5 0.9575 | 1.4408 | 1.0531 0.6886 | 0.7251 
0.6 | 0.9856 | 1.8249 | 1.0649 | 1.1281 | 1.2014 
0.7 | 0.9974 | 2.4947 | 1.0490 | 2.4727 | 2.5938 
0.8 | 0.99997/ 3.9238 | 1.021 | 11.48 11.72 
0.9 | | 1243 


ee 8.5394 1.004 1238 





¢(u) is given by (33) with the constants as in Table II. 
For a near zero, v approaches zero, and 
sinh v(z + 20) = v(z + Zo) 
G3(», Zo) = vE;(z) 


G2(v, 3) = — E2(z) 
= = 1.0915». 
A 


Thus for a = 0, the expression for f(z) determined by Method I is 
f(z) =3{ s + 29 + 0.2248 Es(z) — 0.2455 E,(z)}. (35) 
Method II applied to the case a = 0 gives 
, f(z) = 3 {2 + zo + 0.2244 E;(z) — 0.2436 Ex(z)}. 


In the limiting case a and » near 1, the ratio — also approaches 1. 
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1.1835 
log? (1 — v) 


2s 
a 


We have also 
iad Vl—-a 


4/2 


el-a 








1 3. Vl—a m=a( l-—a ) 
sinh v(z + 29) = e” - —el-a)] + log 2 
0) 2/9 9 g 


and 











ee Be oe ee 
yee 2 REE fp te \ 
av?—a)  Vi-q 2 “< 


Thus in the limit for a = 1, f(z) is 





f(s) = + — (1 + log 2).. \ 2 Z 21 EG) - e*Ey(22)) (36) 


/ 


6 


f(z) / f(0) 
ol 











0 0-5 1°O 1-5 


z 


Fic. 1. Neutron density f(z) /f(o) as a function of distance z from the boundary. Unit of length: 
total mean free path. 





ae er er a CE ER IN 
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The angular distribution y,(u), given by (33) becomes for a = 1, 
l-—a l-—a 1 \ 
sh octet et D ata log(1 ++)+ sauye 
¥n(u) = TF 2-4 r ili ; 
These forms agree with those obtained by the rigorous treatment for a 
For a = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1, the density and the 
asymptotic density are tabulated, in Tables III and IV, asa function of z, the 


distance from tne boundary measured in units of total mean free path. In 
Table III the normalization used is unit outgoing current at the boundary, 


= 1 


so that oo 
f(z) = MA — { sinh »(z + go) + AGa(0, 2) + BG2(v, 2)}. (37) 


aj? —a 


In Table IV the normalization used is unit density at the boundary: 


7) = 4/ £0. (38) 


In each table column (a) gives the density, column (0) the asymptotic density, 
and column (c) the deviation of the density from the asymptotic density. The 
figures used for a = 0 are those obtained by numerical integration (7). 


Fig. 1 isa plot of density normalized to unit density at the boundary (column 
(a), Table IV) over the region within one and a half mean free paths of the 
boundary, for a = 0, 0.1, 0.2, 0.4, 0.6, 0.8, 1. 


0-8 
0°6 
0°4 
0-2 


0 





Fic. 2. Deviation of neutron density from the asymptotic density, normalized to unit deviation 
at the boundary, as a function of distance z from the boundary. Unit of length: total mean free path. 
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The deviation of the true density from the asymptotic density is investi- 
gated further in Table V and Fig. 2. Table V gives the deviation normalized 


to unit deviation at the boundary. The values used for a = 0 in this table were 
obtained from (35). Although, as seen from Tables II and III the deviation 


from the asymptotic is identically zero for a 


1, it is seen from (36) that the 


deviation normalized to unit deviation at the boundary, approaches, as a > 1, 
e?E,(2z) — Ey(z) ; 
— log 2 
Fig. 2 is a plot of this deviation, for a = 0, 0.1, 0.3, 0.5 0.7, 0.8,1. Thecurve 


the limit 


for a = 0 is lowest, the curves rise as a increases until a 
for larger a, so that the limiting curve a 


a = 0.4. 


0.7, and then fall again 
1, lies slightly below the curve for 


The angular distribution of emerging neutrons, ¢(u) is tabulated in Table VI. 


The values for a 


0 are obtained by numerical integration and the values for 


a = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, are obtained from (33). The normal- 
ization used is 


| 


1 
Ko = 1. 


TABLE III 


DENSITY NORMALIZED TO UNIT CURRENT 


7? 

















a=0 a = 0.1 a = 0.2 
(a) | (by ik yo SB) ot, Se) (a) | 6) | © 

0 | 1.7321 | 2.1312 | 0.3991 | 1.6619 | 2.0340 | 0.3721 | 1.5889 | 1.9320 | 0.3431 
0.1 | 2.1841 | 2.4313 | 0.2472 | 2.0693 | 2.3089 | 0.2396 | 1.9566 | 2.1823 | 0.2257 
0.2 | 2.5485 | 2.7312 | 0.1827 | 2.4114 2.5905 | 0.1791 | 2.2731 | 2.4441 | 0.1710 
- 0.3 | 2.8899 3.0312 | 0.1413 | 2.7396 | 2.8792 | 0.1396 | 2.5834 | 2.7180 | 0.1346 
0.4 | 3.2190 | 3.3312 | 0.1122 | 3.0646 | 3.1761 | 0.1115 | 2.8973 | 3.0056 | 0.1083 
0.5 | 3.5406 | 3.6312 | 0.0906 | 3.3905 | 3.4810 | 0.0905 | 3.2198 | 3.3083 | 0.0885 
1.0 | 5.0955 | 5.1312 0.0357 | 5.1338 | 5.1701 | 0.0363 | 5.0702 | 5.1066 | 0.0364 
1.5 | 6.6153 | 6.6312 | 0.0159 7.2017 | 7.2179 | 0.0162 | 7.5392 7.5559 | 0.0167 
2.0 | 8.1237 | 8.1312 | 0.0075 | 9.7593 | 9.7669 | 0.0076 |10.9604 |10.9684 0.0080 
3.0 }11.1294 |11 1229 | 0.0019 |21.7871 0.0020 





1312 | 0.0018 |17.1210 |17. 


e 7891 



















































a = 0.3 a = 0.4 a = 0.5 
(a) -}) 9) (c) (a) (db) (c) (a) (db) 

0 | 1.5130 | 1.8246 | 0.3116 | 1.4347 | 1.7116 0.2769 | 1.3543 | 1.5921 
0.1 | 1.8420 | 2.0512 | 0.2092 | 1.7249 | 1.9146 | 0.1897 | 1.6067 | 1.7727 
0.2) 2.1311 | 2.2915 | 0.1604 | 1.9860 | 2.1333 | 0.1473 | 1.8393 | 1.9696 
0.3 | 2.4202 | 2.5477 | 0.1275 | 2.2514 | 2.3696 | 0.1182 | 2.0789 | 2.1845 
0.4 | 2.7178 | 2.8214 | 0.1036 | 2.5284 2.6254 | 0.0970 | 2.3323 | 2.4196 
0.5 | 3.0293 | 3.1146 | 0.0853 | 2.8228 | 2.9031 | 0.0803 | 2.6037 | 2.6767 
1.0 | 4.9103 | 4.9468 | 0.0365 | 4.6665 | 4.7018 | 0.0353 | 4.3538 | 4.3868 
1.5 7.6231 | 7.6403 | 0.0172 | 7.4682 | 7.4852 | 0.0170 | 7.1053 | 7.1216 
2.0 |11.6557 {11.6642 | 0.0085 11.8271 [11.8358 | 0.0087 |11.5116 |11.5201 
3.0 (26.8575 |26.8599 | 0.0024 |29.3982 |29.4006 | 0.0024 |30.0449 (30.0473 


(c) 


0.2378 


0.1660 
0.1303 
0.1056 
0.0873 
0.0730 
0.0330 
0.0163 
0.0085 
0.0024 
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TABLE I1I1—Concluded 
DENSITY NORMALIZED TO UNIT CURRENT—concluded 



































| a = 0.6 a = 0.7 a = 0.8 
e | | | | | | 
| ®@ | o | oe] @ |] @ | © |@ | |] © 
0 | 1.2724 | 1.4656 | 0.1932 | 1.1922 | 1.3341 | 0.1419 | 1.1179 | 1.2024 | 0.0845 
0.1 | 1.4889 | 1.6256 | 0.1367 | 1.3750 | 1.4760 | 0.1010 | 1.2695 | 1.3290 | 0.0595 
0.2 | 1.6931 | 1.8014 | 0.1083 | 1.5522 | 1.6324 | 0.0802 | 1.4217 | 1.4688 | 0.0471 
0.3 | 1.9061 | 1.9946 | 0.0885 | 1.7394 | 1.8052 | 0.0658 | 1.5849 | 1.6233 | 0.0384 
0.4 | 2.1337 | 2.2073 | 0.0736 | 1.9410 | 1.9959 | 0.0549 | 1.7622 | 1.7941 | 0.0319 
0.5 2.3795 2.4415 | 0.0620 | 2.1601 | 2.2065 | 0.0464 | 1.9560 | 1.9829 | 0.0269 
1.0 | 3.9963 | 4.0224 | 0.0288 | 3.6171 | 3.6290 | 0.0219 | 3.2568 | 3.2693 | 0.0125 
1.5 | 6.5855 | 6.6000 | 0.0145 | 5.9843 | 5.9955 | 0.0112 | 5.3838 | 5.3901 | 0.0063 
2.0 |10.8057 10.8133 | 0.0076 | 9.8684 | 9.8743 | 0.0059 | 8.8832 | 8.8866 | 0.0034 
3.0 [28. 9852 |28.9874 | 0.0022 26.7729 26.7747 | 0.0018 |24.1531 2a! 1541 0.0010 
| | | 
a=0.9 a=1.0 
| @ | ® | © | @ , 6) | © 
0 | 1.0541 | 1.0907 | 0.036 | 1.0000 | 1.0000 | 0 
0.1 | 1.1799 | 1.2054 | 0.0255 | 1.1052 | 1.1052 | 0 
0.2 | 1.3122 | 1.3321 | 0.0199 | 1.2214 | 1.2214 | 0 
0.3 | 1.4561 | 1.4723 | 0.0162 | 1.3499 | 1.3499 | 0 
0.4 | 1.6138 | 1.6271 | 0. 0133 | 1.4918 | 1.4918 | 0 
0.5 | 1.7871 | 1.7982 | 0.0111 | 1.6487 | 1.6487 | 0 
1.0 | 2.9597 | 2.9647 | 0. 0050 | 2.7183 | 2.7183 | 0 
1.5 | 4.8855 | 4.8880 | 0.0025 | 4.4817 | 4.4817 0 
2.0 | 8.0577 | 8.0590 | 0.0013 | 7.3891 | 7.3891 | 0 
3.0 |21. 


9063 |21.9067 0.0004 |20.0855 |20.0855 | 0 


= 
a 

~ 
ll 

a 
to 


(b) = Ws = sinh v(s + 20) 


(b) — (a). 


— 


TABLE IV 
DENSITY NORMALIZED TO UNIT DENSITY AT THE BOUNDARY 














a=0 a= 0. 1 a=0.2 
(a) (b) (c) (a) (b) | (c) iy: a) (b) (c) 
es AN seme | a 
0 | 1.0000 | 1.2304 | 0.2304 | ,. 0000 | 1.2239 | 0.2239 | 1.0000 | 1.2159 | 0.2159 
0.1 | 1.2611 | 1.4037 | 0.1426 | 1.2451 | 1.3893 | 0.1442 | 1.2314 | 1.3735 | 0.1421 
0.2 | 1.4714 | 1.5769 | 0.1055 | 1.4510 | 1.5588 | 0.1078 | 1.4306 | 1.5382 | 0.1076 
0.3 | 1.6685 | 1.7501 | 0.0816 | 1.6485 | 1.7325 | 0.0840 | 1.6259 | 1.7106 | 0.0847 
0.4°| 1.8585 | 1.9233 | 0.0648 | 1.8440 | 1.9111 | 0.0671 | 1.8235 | 1.8916 | 0.0681 
0.5 | 2.0442 | 2.0965 | 0.0523 | 2.0401 | 2.0946 | 0.0545 | 2.0264 | 2.0821 | 0.0557 
1.0 | 2.9419 | 2.9625 | 0.0206 | 3.0891 | 3.1110 | 0.0219 | 3.1910 | 3.2139 | 0.0229 
1.5 | 3.8193 | 3.8285 | 0.0092 | 4.3334 | 4.3432 | 0.0098 | 4.7449 | 4.7554 | 0.0105 
2.0 | 4.6902 | 4.6946 | 0.0044 | 5.8724 | 5.8769 | 0.0045 | 6.8981 | 6.9031 | 0.0050 
3.0 | 6.4256 | 6.4266 | 0.0010 110. 3021 [10.3032 | 0.0011 |13.7121 |13.7133 | 0.0012 
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TABLE IV—Concluded 
DENSITY NORMALIZED TO UNIT DENSITY AT THE BOUNDARY—concluded 
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) = —* 
GC Ce 





- p(o) 


= 0.3 = 0.4 a = 0.5 

. (a) )) (c (a) (bh) (c) (b) (c) 

0 | 1.0000 2059 | 0.2059 | 1.0000 _ 1.1930 | 0.1930 | 1.1756 0.1756 
0.1 | 1.2174 3557 | 0.1383 | 1.2023 | 1.3345 | 0.1322 | 1.3089 | 0.1225 
0.2 | 1.4085 .5145 | 0.1060 | 1.3843 | 1.4869 | 0.1026 | | 1.4543 | 0.0962 
0.3 | 1.5996 .6839 | 0.0843 | 1.5692 | 1.6516 | 0.0824 | 1.6130 | 0.0780 
0.4 | 1.7963 8648 | 0.0685 | 1.7623 | 1.8299 | 0.0676 | | 1.7866 0.0645 
0.5 | 2.0022 0586 | 0.0564 | 1.9675 | 2.0235 | 0.0560 | | 1.9764 0.0539 
1.0 | 3.2454 .2695 | 0.0241 | 3.2526 | 3.2772 | 0.0246 | 3.2392 | 0.0244 
1.5 | 5.0384 0498 | 0.0114 | 5.2054 | 5.2173 | 0.0119 | 5.2585 | 0.0120 
2.0 | 7.7037 7093 | 0.0056 | 8.2436 | 8.2497 | 0.0061 ; | 8.5063 | 0.0063 
3.0 {17.7512 |17.7527 | 0.0015 .4908 20.4925 | 0.0017 |22.1848 |22.1866 | 0.0018 

= 0.6 = 0.7 = 0.8 
(a) (b) (c) (a) (b) (c) (b) (c) 

0 1.0000 | 1.1518 | 0.1518 | 1.0000 | 1.1190 | 0.1190 | 1 | 1.0756 | 0.0756 

.1 | 1.1702 | 1.2776 | 0.1074 | 1.1533 | 1.2380 | 0.0847 | 1. 1.1888 | 0.0532 

.2 | 1.3306 | 1.4157 | 0.0851 | 1.3020 | 1.3692 | 0.0672 | 1 | 1.31389 0.0421 

.3 | 1.4980 | 1.5676 | 0.0696 | 1.4590 | 1.5142 | 0.0552 1 | 1.4521 | 0.0344 

t | 1.6769 | 1.7348 | 0.0579 | 1.6281 | 1.6741 | 0.0460 | 1 | 1.6049 | 0.0286 

5 | 1.8701 | 1.9188 | 0.0487 | 1.8119 | 1.8508 | 0.0389 | 1. 1.7738 | 0.0241 

.0 3.1386 | 3.1613 | 0.0227 | 3.0340 | 3.0523 | 0.0183 | 2 2.9245 | 0.0112 

.5 | 5.1757 | 5.1870 | 0.0113 | 5.0195 | 5.0289 | 0.0094 | 4 4.8216 | 0.0056 

.0 | 8.4924 | 8.4983 | 0.0059 | 8.2775 8.2824 | 0.0049 | 7 | 7.9494 | 0.0031 

.O |22.7799 |22.7817 | 0.0018 |22.4567 22.4582 | 0.0015 (21.6058 |21.6067 | 0.0009 

a = 0.9 = 1.0 
(a) (b) (c) (a) (b) (c 

0 | 1.0000 | 1.0347 | 0.0347 | 1.0000 | 1.0000 0 

.1 | 1.1193 | 1.1435 | 0.0242 | 1.1052 | 1.1052 0 

.2 | 1.2449 | 1.2637 | 0.0188 | 1.2214 | 1.2214 0 

3 1.3814 1.3967 | 0.0153 1.3499 | 1.3499 0 

.4 |} 1.5310 | 1.5436 | 0.0126 | 1.4918 | 1.4918 0 

.5 | 1.6954 | 1.7059 | 0.0105 | 1.6487 | 1.6487 0 

0 | 2.8078 | 2.8125 | 0.0047 | 2.7183 | 2.7183 0 

.5 | 4.6348 | 4.6371 | 0.0023 | 4.4817 | 4.4817 0 

0 | 7.6442 | 7.6454 | 0.0012 | 7.3891 | 7.3891 0 

0 20.7820 20.7824 | 0.0004 20.0855 |20.0855 0 

(a) = f(z) 

(b) = (= 4) sinh V(s + Zo) 
4 ae inh + Zo 

(c) = (b) —(a). 


This normalization would be equivalent to unit density at the boundary if 
¢(u) were the correct angular distribution. Since 





i 
j 
| 
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it is evident from the values of ¢(0) given in Table VI that the accuracy of 
1 


the approximation to ¢(u) is such that | $(u)duis practically equal to its true 
value. 


Fig. 3a gives a plot of ¢(u) for a = 0, 0.1, 0.2, 0.3, and Fig. 3b gives $(u) 
for a = 0, 0.3, 0.4, 0.5, 0.6. In the limit as a > 1, @(u) is zero for wu ¥ 1, and 
for » = 1, becomes infinite. 


The accuracy of the approximation to the angular distribution is indicated 
in Tables VII and VIII. Table VII shows, for » = 1, and wu = 0,anda = 0.1, 
0.2, 0.4, 0.6, the values of ¢(u) obtained by the method of variations, and the 


g (py) 








0 0-2 0-4 0-6 08 0 
pL 


Fics. 3a AND 3b. Emerging angular distribution o(u), normalized to unit density: 


f,eWdu = 1. 
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TABLE V 





DEVIATION FROM ASYMPTOTIC SOLUTION 
NORMALIZED TO UNIT DEVIATION AT THE BOUNDARY 


AS A FUNCTION OF DISTANCE FROM THE BOUNDARY 





VOL. 28, SEC. A. 






values obtained by numerical integration (2). Table VIII shows the values of 
¢(u) obtained by the method of variations and those obtained by integration 
(2), for a = 0.4 and p = O, 0.4, 0.8, 1. 

















| 
| 
| 
| 
| 





ANGULAR DISTRIBUTION OF NEUTRON EMERGING AT THE BOUNDARY 
NORMALIZED TO UNIT 


a=0.1 |a=0.2 





a=0.3 


}a=0.4 


TABLE VI 





.8213 
. 9146 
.0064 
.0971 
. 1870 
. 2764 
. 3653 
. 4008 
.4185 
. 4362 
. 4539 





0. 
0. 
0. 
0. 
0.: 
0.3 
0. 
0. 
06 
0. 
0.8 
0.¢ 
0.9 
0. 
0. 
ai 





0.4502 





.4642 | 0 
-5077 | 0 
.5562 | 0 
-6490 | 0 


-9497 
. 0686 
. 2027 
. 3569 
. 5378 
6196 | 
6629 | 
7079 | 
7548 | 


0 

0 

0 
0.6490 | 
0.7430 | 0. 
0.8423 | 
0 

1 

1 

13560 | 

1 | 
:- } 

1 
1 
1 


NKR eee ROO 


.4001 
4117 
4483 
.4900 


.5725 


6609 
7604 


.8774 
.0172 | 
. 1933 
.4234 


7402 
9037 


.9963 


0973 


. 2081 


| 0.3501 


0.3595 | 


| 0.3896 
0.4245 

| 0.4959 
| 0.5757 
0.6703 
0.7875 
.9398 
1.1487 
4568 
.9621 
. 2692 


NNwNNRK 


6840 
9514 


4599 | 


DENSITY AT THE 


}a=0.5 


0.2500 | 0.: 


ROUNDARY 


}a=0.6 |a=0.7 





0.3000 2000 | 0 
0.3074 | 0.2556 | 0.2040 | 0.1526 
0.3315 | 0.2741 | 0.2175 | 0.1618 | 
0.3599 | 0.2964 | 0.2337 | 0.1733 
0.4195 | 0.3442 0.2705 | 0.1989 
0.4887 | 0.4012 | 0.3146 | 0.2304 
0.5738 | 0.4732 | 0.3715 | 0.2716 
0.6841 | 0.5696 | 0.4492 | 0.3284 
0.8354 | 0.7076 | 0.5633 | 0.4129 | 
1.0598 | 0.9239 | 0.7488 | 0.5527 | 
1.4297 | 1.3151 | 1.1058 | 0.8296 
2.1622 | 2.2456 | 2/0840 | 1.6455 
2.7073 | 3.1160 | 3.2091 2.7013 
3.0946 | 3.8605 | 4.3896 | 3.9728 | 
3.6087 | 5.0673 | 6.9361 | 7.4954 | 
4.3237 7.3616 |16.464 65.355 








.1500 | 








.1000 | 
1023 

-1073 
(1143 | 
'1299 | 


1498 


.1759 
2124 
. 2666 
3565 
5364 
0751 
7925 
2.6879 
5.3659 
186.5 


a=0 |ja=01j)/a=02);a =03)/a = 04/0 = 0.5 
0 | 1.0000 1.0000 | 1.0000 | 1.0000 | 1.0000 1.0000 
0.1 | 0.6286 | 0.6435 | 0.6576 | 0.6712 | 0.6852 | 0.6979 
0.2 | 0.4635 0.4810 | 0.4977 | 0.5150 | 0.5317 0.5476 
0.3 0.3573 0.3749 | 0.3920 | 0.4091 0.4272 | 0.4441 
0.4 | 0.2825 | 0.2992 | 0.3158 | 0.3326 | 0.3501 0.3668 
0.5 0.2271 | 0.2426 | 0.2577 0.2739 | 0.2904 0.3068 
0.7 0.1520 | 0.1645 | 0.1774 0.1911 | 0.2055 | 0.2199 
1.0 | 0.0880 | 0.0971 0.1063 | 0.1166 | 0.1276 | 0.1389 
1.5 | 0.0386 | 0.0437 0.0486 | 0.0549 | 0.0615 | 0.0686 
1.8 0.0244 0.0279 | 0.0316 | 0.0359 0.0408 | 0.0460 
2.0 | 0.0181 0.0209 | 0.0237 0.0273 | 0.0312 | 0.0355 
3.0 | 0.0045 | 0.0053 | 0.0062 | 0.0074 | 0.0087 | 0.0102 
| | | 
ja =06,a=07);a =08);a = 09} a =1 
0 | 1.0000 | 1.0000 1.0000 1.0000 1.0000 
0.1 | 0.7077 0.7114 0.7044 | 0.6970 | 0.6805 
0.2 | 0.5604 0.5655 0.5570 | 0.5454 | 0.5262 
0.3. | 0.4580 | 0.4637 0.4550 | 0.4416 | 0.4217 
0.4 0.3809 | 0.3870 | 0.3783 | 0.3645 | 0.3448 
0.5 0.3204 | 0.3267 0.3185 | 0.3048 | 0.2858 
0.7 | 0.2323 | 0.2383 | 0.2313 | 0.2186 | 0.2016 
1.0 | 0.1489 | 0.1541 0.1487 | 0.1384 0.1247 
1.5 | 0.0752 0.0788 | 0.0757 0.0689 | 0.0599 
1.8 | 0.0509 0.0538 | 0.0516 | 0.0464 | 0.0396 
2.0 | 0.0396 0.0418 0.0402 0.0359 | 0.0303 
3.0 0.0117 0.0128 | 0.0122 | 0.0106 | 0.0084 






j}a=0.8 |a=0.9 


0500 
0506 


.0530 
0562 
. 0637 
.0731 
.0855 
. 1029 
. 1289 
.1722 
. 2586 


5178 


. 8635 
. 2955 
.5915 


. 2578.10? 
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TABLE VII 
ANGULAR DISTRIBUTION, # = 0 AND p = 1 























A= 0 = l 
. Value obtained True Value obtained | Value obtained 
by variation value by variation by integration 
0.1 0.4502 | 0.4500 1.7548 1.7542 
0.2 0.4001 0.4000 2.2081 2.2075 
0.4 0.3000 | 0.3000 4.3237 4.3232 
0.6 0.2000 0.2000 16.4640 | 16.4629 
TABLE VIII 
ANGULAR DISTRIBUTION, a = 0.4 





bl Value obtained by Value obtained 
method of variation by integration 
0 0.3000 0.3000 
0.4 0.5738 0.57407 
0.8 1.4297 1.4301 
1 4.3237 4.3232 


Milne Problem with Constant Production 


It has been shown by Davison and Placzek (4) that the density and angular 
| distribution in the Milne problem with constant production are simply related 
to the density and angular distribution in the Milne problem with the same 
| capture and no production. We use the notation of that paper, namely: f(z) 
and y,(u) denote the density and outgoing angular distribution in the Milne 
| problem with no production normalized to unit outgoing current; p,(z) and 
W,(z) denote the density and angular distribution in the Milne problem with 
production normalized to unit production; g(z) denotes the fictitious source 


density, g(z) = (1 — a) p(s) + 1. (39) 


From Davison and Placzek, we have 


g(z) = g(0) V4 FH) - | sooae, (40) 


v 


where the value of g(o) is arbitrary and depends on the conditions at infinity. 
In particular, if g(z) is bounded at infinity, 


] 
oo 41 
g(o) > (41) 
and abi ate | “f(t)dt. (42) 
: v 0 


Also, under the same assumption, 


dln) = oer (i ed (43) 


(1 — a)p 
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Using for f(z) the formula (37) obtained by the variational method, 


a 2 — v) ena Lea ete) 
we WG LEP 5 Raa AOS 
g(z) 1 4/20”) : es {vA E:(z) 


a(v? — my 
+ (A — B)E,(z) — (A — B).~” Ex((1 — »))} (44) 


where v, Zo, A, and B are given in Table II. 


TABLE IX 
FICTITIOUS SOURCE DENSITY, MILNE PROBLEM WITH CONSTANT PRODUCTION 
NORMALIZED TO UNIT PRODUCTION 


0.1 a = 0.2 





(c) 


(b) (c) | (2) 


0112 | 0.8489 2. 8560 | 0.6199 
| 
| 


0. 4281 
0.3313 
0. 2650 
0.2160 
0.1783 


.3178 | 0.5679 | 3.0031 
6086 | 0.4321 | 2. | 3.1400 
8846 | 0.3412 .0026 | 3.2676 
1464 | 0.2750 | 3. | 3.3864 
3949 | 0.2249 | 3. 3.4970 
4588 | 0.0926 | 3. 3.9464 | 0.0762 
2770 | 0.0424 | 4.2256 | 4.2614 | 0.0358 
9061 | 0.0205 | 4.46 4.4822 | 0.0177 
.7619 | 0.0054 | 4. | 4.7455 | 0.0047 
0.0000 | 5. | 5.0000 | 0.0000 





WHE MSSSOS_ 
COMCURwWNH | 
See 
I “J 
SOND EE ee | 


_ 


8 








a = 0.4 








(c) 


2.0102 0. 4291 
2.0526 0.3111 
2.0915 0.2477 
2. 1269 0. 2025 
2.1593 0. 1693 
2. 1888 0. 1413 
2.3023 0. 0644 
2.3744 0.0318 
2. 4202 0.0164 
2. 4678 0.0047 
2. 5000 0.0000 


o 


4129 | 0.3596 
.4861 | 0.2828 
.5535 | 0.2285 


2. 
0533 | 2. 
2 
2 | 
2.6155 | 0.1883 
2 
2 
3 


| 3334 | 0.5077 
2033 
3250 | 
4272 | 
5158 
8275 | 


6726 | 0.1568 
8967 | 0.0692 | 
0115 | 3.0448 | 0.0333 
1259 | 3.1427 | 0.0168 
2454 | 3.2501 | 0.0047 
3.3333 | 0. 


oo ee re eo 
COMNCURWNHH 





G9 29 OO SS KORO NNN | 
zt I 


8 











| 
| 
i 
| 








(c) 


0.2919 
0.2190 
0.1779 
0.1481 
0.1249 
0. 1063 
0.0513 
0.0265 
0.0141 
0.0043 
0.0000 








. 6980 : 

7349 . 8607 
. 8546 . 9137 
9167 . 9465 
. 9512 . 9669 
.9827 | . 9873 
.0000 | 


es SSeS 
cCoonocourwn- 
De tt et et et et et et 


| 
BND a reek et et pee pe pet fet et et | 
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NEUTRON DENSITY 


TABLE IX—Concluded 


FICTITIOUS SOURCE DENSITY, MILNE PROBLEM WITH CONSTANT PRODUCTION 
NORMALIZED TO UNIT PRODUCTION—concluded 


























a = 0.7 a = 0.8 

(a) (b) (c) (a) (db) | (c) 
0 1. 1952 1. 4100 0.2148 1.1180 1.2491 | 0.1311 
0.1 | 1.2500 1.4118 0. 1618 1.1514 1.2491 0.0977 
0:2. | 1.9818 1.4133 0. 1321 1.1698 | 1.2492 0.0794 
0.3 1.3045 1.4148 ). 1103 1.1832 | 1.2493 | 0.0661 
0.4 1. 3228 1.4161 0.0933 1.1985 | 1.2494 | 0.0559 
0.5 | 1.3376 1.4173 0.0797 1.2019 | 1.2494 | 0.0475 
1.0 1.3827 1.4217 0.0390 1. 2265 1.2497 0.0232 
1.5 1. 4040 1.4244 0.0204 1. 2378 1. 2498 0.0120 
2.0 1.4150 1.4261 0.0110 1.2434 | 1.2499 | 0.0065 
3.0 1.4242 1.4276 0.0034 1.2479 | 1.2500 | 0.0020 
oo 1.4286 | 1.4286 0.0000 1.2500 | 1.2500 | 0.0000 

a = 0.9 a=] 

(a) (b) (c) (a) (b) (c) 
0 1.0541 1.1111 0.0570 1 1 0 
0.1 | 1.0689 1.1111 0. 0422 1 1 0 
O32 | era 1.1111 0.0340 1 1 0 
0.3 1.0830 1.1111 0.0281 1 1 0 
0.4 1.0876 1.1111 0. 0235 1 1 0 
0.5 1.0912 1.1111 0.0199 1 1 0 
1.0 | 2.1017 1.1111 0.0094 1 1 0 
1.5 | 1.1063 1.1111 0.0048 1 1 0 
2.0 | 1.1086 1.1111 0.0025 1 1 0 
oo.) ie 1.1111 0. 0008 1 1 | 0 
0 1.1111 1.1111 0. 0000 1 1 0 


(a) = g(z), (b) = gas(z), (c) = (b)—(a). 


The values of g(z) determined from (44) are tabulated in Table IX, column (a). 
Column (0) of Table IX daa 


as(2) = Lae ee (45) 


a(v? — ae 





and column (c) gives the deviation of g(z) from ga;(z). For a = 0, g(z) is 
infinite for every z, and for a = 1, g(z)= 1. Fig. 4 is a plot of g(z) 


The density p,(z) can be found easily from (39), which is equivalent to 


(z) — 1 
Pq(2) = — : 
p_(z) is tabulated in Table X column (a). Column (0) gives 


1 1 BOB Pa ein So 
reste fn 
a l—a a(v? — a) 


and column (c) gives the deviation from the asymptotic, Pags(2) — Pq(z). For 
a = 0, p,(z) is infinite, and for a = 1 





(46) 
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ss 
p4(z) =1 — 2 F(z). (47) 


Formula (47) is the limit as a > 1 Toe with g(z) given by (44). It can 
= @ 


easily be shown to be the exact value of p,(z) fora = 1. Fig. 5isa plot of p,(z). 


TABLE X 
DENSITY, MILNE PROBLEM WITH CONSTANT PRODUCTION NORMALIZED 
TO UNIT PRODUCTION 


(c) 


3. 0.9433. | 1.5451 | 2.3200 
3.0554 a 0.6310 | 1.9687 | 2.5038 
3.5294 4. 0096 0.4802 | 2.2609 | 2.6750 
3.9372 | 4.3162 | 0.3790 | 2.5032 | 2.8344 
4.6 
4. 
6. 
6. 
7 





o 


4.3016 0.3055 | 2.7130 | 2.9830 

4. 6333 0. 2499 2.8984 | 3.1213 

5.9624 0. 1030 3.5878 | 3.6830 

6. 9273 0.0471 4.0320 4.0767 

7.6506 | 0.0228 | 4.3306 4.3527 

| 8.6184 8. 0.0059 | 4.6760 | 4.6819 
10. 0000 10. 0000 0.0000 | 5 


Sit Hao S 


oovourwn-e 
osessssosse | 


a = 0.3 a 

bi (a) | (b) 

o | 1.1796 9049 | 0.7253 | 0.9686 | 1.6836 
1.5047 0185 | 0.5138 | 1.2358 | 1.7544 
1.7191 1231 ‘4040 | 1.4064 "8191 
1. 8929 2193 "3264 | 1.5406 | 1.8781 


(c) aye <4 (b) 


2.0389 3079 .2690 | 1.6500 9321 
| 1.7459 | 1.9813 
7096 .0989 | 2.0631 1705 
9212 | 0.0477 | 2.2377 | 2.2907 
0610 | 0.0241 | 2.3397 . 3670 
2144 | 0.0066 | 2.4386 | 2.4463 
3333 | 0.0000 | 2.5000 | 2.5000 


2.6107 
2.8735 | 
3.0369 | 
3.2078 
3.3333 


WNKrROoSSSS 


CcCoucutwnre 
Seessosc—es9 | 


: 


i. 
2. 
2. 
2. 
2. 
2. 1654 2.3894 . 2240 
2. 
2. 
3. 
3. 
3. 


8 


a = 0.5 


(a) (b) (c) (a) (c) 


l 
| 
| a 
| 


1.5503 | 0.7219 | 0.7275 | 1. | 0.7297 
1.5913 | 0.5338 | 0.9295 | 1.4769 | 0.5474 
1.6286 | 0.4297 | 1.0499 | 1 | 0.4448 
1.3051 | 1.6600 | 0.3549 | 1.1407 | 1 | 0.3701 
1.3961 | 1.6934 | 0.2973 | 1.2133 | 1.5255 | 0.3122 
1.4699 1.7214 | 0.2515 | 1.2730 | 1.5387 | 0.2657 
1.7092 | 1.8274 | 0.1182 | 1.4603 | 1.5885 | 0.1282 
1.8333 | 1.8930 | 0.0597 | 1.5528 | 1.6189 | 0.0661 

1. 1. 1 | 

1. 1. 1 

2. 1. 1 | 


| 0.8284 
1.0575 
1. 1989 


1rOoUFWNe 


. 6375 
. 6558 
. 6667 


Nee Ooooo 


1.9024 6022 | . 0353 
6451 | 
6667 


9337 0.0313 
9746 0.0092 
0000 0. 0000 


eoog 


- 0107 
0.0000 


1. 9654 
2.0000 


w 
ae 


8° 





DENSIT 


CANADIAN JOURNAL OF RESEARCH. 


VOL. 28, SEC. 


TABLE X—Concluded 


y, MILNE PROBLEM WITH CONSTANT PRODUCTION NORMALIZED 
7O CHET ERE TOR — sere 


A, 





= 0.7 


a 








0. 6508 


( 


rOoFWNe 


oS 


Gee SSeS 
oo 


o 


). 8334 
. 9378 
.0151 
.0759 
. 1252 
. 2756 
. 3467 
. 3835 
. 4141 
. 4286 


Pree fe et fh fret feet et et 


Ls 
1. 
. 3778 
- 3826 
. 3870 
. 3910 
. 4057 
. 4147 
. 4202 
. 4255 
. 4286 


3666 
3725 


(c) 





| 0.7159 | 0. 
0.5392 | 0. 
0.4400 | 0. 

5 0. 
0. 


5902 
7568 
8489 
9158 
9677 


. 0093 
. 1327 
. 1888 
. 2169 
. 2397 
. 2500 


0. 6551 
0. 4889 
0.3972 
0. 3307 
0.2791 
0. 2378 
0. 1156 
0.0601 
0.0325 
0.0100 
0. 0000 








| 

i 

| 
Lp 


Qo 


0. 
0. 
0. 
0. 
0. 
be 
A 
2. 
3. 


oonourwnd- 


8 


DENsITY, MILNE 


Vere ere > 


1 
.2784 | 1. 
4715 | 1 
6686 | 1 
8587 | 1. 
0443 | 2. 
.9421 | 2 
.8197 | 3 
.6906 | 4 
4261 | 6. 

@ 


> Hm C9 DD DD et ee ee | 


Het SSSOSS | 


(a) 


. 5409 
6885 
7712 
. 8304 

759 
9120 
0169 
. 0631 
. 0856 
. 1034 


(c) 


8 


(c) 





Reider) 
1.1111 
1 33i2 
1.1111 
1.1111 
1. 
1 
1 
1 
1 
1 


1111 


sath 
ae 
-1111 
it 
head 


0.5702 
0. 4226 
0.3399 
0. 2807 
0. 2352 
0.1991 
0.0942 
0.0480 
0.0255 
0.0077 
0.0000 


fet pet fre peed fet feet feet fet fet fed fet 





PROBLEM 


2306 
4038 


0967 


2. 9628 
3.8289 
. 6950 


4272 


co 


(b) = Pq,,(2)» (c) = (6) — (@). 


TABLE XI 


AT THE BOUNDARY 
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DENSITY, MILNE PROBLEM WITH CONSTANT PRODUCTION NORMALIZED TO UNIT DENSITY 
AT THE BOUNDARY—concluded 
a = 0.3 a = 0.4 a= 0.5 
, a) | () (c) (a) (b) () @) | ©@® (c) 
- | c ane ——_ | 
0 1.0000 | 1.6148 | 0.6148 1.0000 | 1.7382 | 0.7382 | 1.0000 | 1.8713 | 0.8713 
0.1 | 1.2756 | 1.7111 | 0.4355 | 1.2758 | 1.8114 | 0.5355 | 1.2766 | 1.9209 | 0.6443 
0.2 | 1.4573 | 1.7998 | 0.3425 1.4520 | 1.8781 | 0.4261 | 1.4472 | 1.9659 | 0.5188 
0.3 | 1. 6047 1.8814 0. 2767 | 1.5906 1.9391 | 0.3485 | 1.5754 | 2.0038 | 0.4284 
0.4 | 1.7284 | 1.9565 | 0.2281 | 1.7036 | 1.9948 | 0.2912 | 1.6852 2.0441 | 0. 3589 
0.5 | 1.8357 | 2.0256 | 0.1899 | 1.8026 | 2.0457 | 0.2431 | 1.7743 | 2.0779 | 0.3036 
1.0 | 2.2131 2.2970 0.0839 | 2.1301 | 2.2409 | 0.1108 | 2.0631 | 2.2058 | 0. 1427 
1.5 | 2.4859 2.4764 | 0.0405 | 2.3103 2.3650 | 0.0547 | 2.2130 | 2.2851 | 0.0721 
2.0 | 2.5745 | 2.5949 | 0.0204 | 2.4157 | 2.4438 | 0.0281 | 2.2964 | 2.3342 | | 0.0378 
3.0 | 2.7193 | 2.7249 | 0.0056 | 2.5177 | 2.5257 | 0.0079 | 2.3724 | 2.3835 | 0.0111 
co | 2.8258 | 2.8258 | 0.0000 | 2.5810 | 2.5810 0.0000 | 2.4143 | 2.4143 | 0.0000 
| 
a = 0.6 a = 0.7 a = 0.8 
~ | @ (6) (c) (a) >) (c) (a) | () (c) 
— | a - = _ | ~aeene as | eomesseemnesi | aoe — 
0 | 1.0000 | 2.0031 | 1.0031 | 1.0000 | 2.1000 | 1.1000 | 1.0000 | 2.1100 | 1.1100 
0.1 | 1.2776 2.0301 | 0.7525 | 1.2806 | 2.1091 | 0.8285 | 1.2823 | 2.1108 | 0.8285 
0.2 | 1.44382 | 2.0546 | 0.6114 | 1.4411 | 2.1173 | 0.6762 | 1.4883 | 2.1115 | 10. 6732 
0.3 | 1.5680 2.0768 | 0.5088 | 1.5598 | 2.1247 | 0.5649 | 1.5517 | 2.1121 | 0 . 5604 
0.4 | 1.6677 | 2.0969 | 0.4292 | 1.6532 2.1314 | 0.4782 | 1.6397 | 2.1127 0. 4730 
0.5 1.7498 | 2.1151 | 0.3653 | 1.7291 | 2.1374 | 0.4083 | 1.7102 | 2.1132 | 0.4030 
1.0 | 2.0074 | 2.1836 | 0.1762 | 1.9601 2.1601 | 0.2000 | 1.9193 | 2.1151 | 0.1958 
1.5 | 2.1345 | 2.2254 | 0.0909 | 2.0694 | 2.1739 | 0.1045 | 2.0143 | 2.1163 | 0.1020 
2.0 | 2.2024 | 2.2509 | 0.0485 | 2.1259 | 2.1823 0.0564 | 2.0620 | 2.1170 | 0.0550 
3.0 | 2.2613 | 2.2760 | 0,0147 | 2.1730 | 2.1905 0.0175 | 2.1006 | 2.1176 | 0.0170 
co | 2.2910 | 2.2910 | 0.0000 | 2.1951 | 2.1951 0.0000 | 2.1179 | 2.1179 | 0.0000 
a = 0.9 a=] 
8 Ge (b) (c) @ | ® | © 
0 | 1.0000 | 2.0541 | 1.0541 | 1.0000 | 2 | 1.0000 
0.1 1.2728 | 2.0541 | 0.7813 | 1.2775 | 2 | 0.7225 
0.2 | 1.4257 | 2.0541 | 0.6284 1. 4258 | 2 | 0.5742 
0.3 i 1.5351 2.0541 0.5190 | 1.53809 2 | 0. 4691 
0.4 1.6192 | 2.0541 0.43849 | 1.6106 2 0. 3894 
0.5 1.6860 | 2.0541 0.3681 | 1.6734 | 2 0.3266 
1.0 1.8799 2.0541 0. 1742 1.8515 | 2 0. 1485 
1:5 1.9652 | 2.0541 0. 0889 1.9269 | 2 | 0.0731 
2.0 2.0070 | 2.0541 0.0471 1.9625 2 0.0375 
3.0 2.0399 2.0541 0.0142 1. 9894 2 | 0.0106 
co | 2.0542 2.0542 0.0000 | 2.0000 2 | 0.0000 
| | 
(a) = Pg(z), (b) = Pag,(2), (c) = (b) — (a). 
It follows also from (46) and (41) that 
] 
alle eee (48) 
. vV ‘a (1 asad 


° ke ° ,. e e = 
If, instead of normalizing to unit production we normalize to unit density at 


the boundary, 
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- 7, _ Ba(2) _Vag(s) -Va ‘ 
pz) = 4 = (49) 
p(0) i< Va 
The values of p,(2), b,,,(2), and p,, (2) — pg (z) are tabulated in columns (a), (0), 
(c) of Table XI, with g(z) given by (44). For a = 0, p,(z) exists and is equal 
1 
V3 
The normalization to unit density is convenient in discussing the angular 
distribution. We have 


to f(z). Fig. 6 is a plot of p,(z). 


1 — py o 
,(u) > es (x), (50) 
l1—vVa 
where ¢(u) is the angular distribution for the Milne problem with no pro- 
duction, normalized to unit density (Table VI). ¢, is tabulated in Table XII 


and plotted in Fig. 7. 
TABLE XII 


EMERGING ANGULAR DISTRIBUTION MILNE PROBLEM WITH CONSTANT PRODUCTION 
NORMALIZED TO UNIT DENSITY AT THE BOUNDARY 


m 
M 
S 


| 
| 


7741 | 0.816 3 0.8873 .9183 | 0.9472 | 0.9745 
7883 | 0.8289 0.8643 | 0.8962 .9251 | 0.9590 | 0.9762 
8257 | 0.8611 0.9174 .9412 | 0.9653 | 0.9830 
S608 0.8905 0.9346 9551 .9740 | 0.9867 
.9147 | 0.9342 0.9635 9748 9844 | 0.9930 
.9565 | 0.9677 0.9830 .9885 .9933 | 0.9971 
9908 | 0.9946 9984 9994 | 0.9997 | 0.9997 
0198 0170 0108 0080 .0059 .0028 
0446 0356 a 1.0212 0151 | 1.0101 | 1.0048 
0666 0521 0300 0212 .0133 | 1.0067 
O858 0663 0376 .0263 | 1.0165 | 1.0079 
. 1029 0792 0443 | 1.0308 .0192 0093 
.1092 0837 | 0465 .0325 0202 | 1.0098 
1124 0859 } 1.0478 | 1.0333 .0206 | 1.0100 
1152 0880 .0487 | 1.0340 .0210 | 1.0102 


1183 0904 .0501 | 1.0347 0215 | 1.0104 
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| 0. 
| 0. 
| 0 
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0 
| 0. 
i; 
oe 
| 1 
| 1 
1 
| 1 
1 
1 
1 


Acknowledgments 


The research in this report was undertaken at the suggestion of Dr. George 
Placzek who offered valuable help and encouragement throughout. The com- 
putations and graphs included were prepared under the direction of Mr. Bengt 
Carlson. The author wishes to express her gratitude to Dr. Placzek and to 
Mr. Carlson and his staff. 

References 


ADLER, F. T. and Mark, C. Phys. Rev. To be published. 
Davison, B. Can. J. Research, A, 28: 303. 1950. 

Davison, B. Phys. Rev. 71: 694. 1947. 

Davison, B. and PLaczEK, G. Phys. Rev. To be published. 
. LECAINE, J. Phys. Rev. 72: 564. 1947. 

LecatnE, J. Montreal Report, MT-119, 1945. 

Mark, C. Phys. Rev. 72: 558. 1947. 

. MarRsHAK, R. E. Phys. Rev. 71: 688. 1947. 


GOVE? Oe GO BO 


(The Appendix follows on the next page) 





CANADIAN JOURNAL OF RESEARCH. VOL. 28, SEC. A. 


APPENDIX 


The explicit expressions for the constants L,, Lo, Ki, Ke, K3, are listed below 
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A GEOMETRIC TREATMENT OF “DIMENSIONS” IN PHYSICS! 


By PARRY Moon AND DoMINA EBERLE SPENCER 


Abstract 


“Dimensional analysis”’ is an established and venerable subject, associated with 
such names as Fourier, Maxwell, and Lord Rayleigh. But only a genius can 
obtain dependable results with it, a state of affairs that is caused by the vague, 
intuitionalistic way in which the subject has always been presented. The pur- 
pose of this paper is to formulate a precise theory. 

An idon or complete physical quantity is an entity with two attributes: 
quality and quantity. Formal rules are given for the manipulation of idons, 
and it is shown that they can be pictured as vectors in a k-dimensional space. 
This geometrization is not necessary, of course, but it often gives added grasp 
of the subject; just as vector analysis is never necessary in physics, though it is 
often convenient. There are six practical applications of the theory, all of which 
are handled in a precise and dependable way by the geometric method. 


Introduction 


The ‘‘dimensional”’ theory of physical quantities seems to have been first 
used by Fourier in his celebrated treatise (9), Théorie analytique de la chaleur 
(1822). It was further developed by Maxwell (11, 12) (1863); and Lord 
Rayleigh (15) was one of its ablest supporters. The subject has been extended 
by Vaschy (20), Riabouchinsky (16, 17), Buckingham (3), Bridgman (1), 
Ehrenfest-Afanassjewa (7, 8), Tolman (19), Campbell (4, 5), Hill (10), and 


others. 


4c 


Despite the effort that has been expended on ‘‘dimensions,”’ however, the 
theory has remained in an unsatisfactory state, as evidenced by the endless 
and largely meaningless disputes that have taken place. (2, 4, 5, 6, 14, 17, 21). 
The difficulty is primarily semantic; for in reading the voluminous and con- 
flicting literature on the subject, one can never be sure of the exact meaning 
attached to such words as ‘‘dimension’”’, “unit’’, ‘‘magnitude’’, ‘‘quantity”’. 
And this difficulty is not solely an academic question. The complaint is often 
made that it takes the physical intuition of a Lord Rayleigh to apply ‘‘dimen- 
sional analysis’ with any certainty, a situation that is caused at least partially 
by the lack of precision in the traditional treatment of the subject. 

Within recent years, there has been progress toward greater exactness in 
“dimension” theory. The present paper attempts to contribute toward that 
end by showing how the subject can be treated as a branch of geometry. 


Fundamentals 


Maxwell stressed the use of ‘“‘dimensions”’ in changing from one set of units 
to another. Thus he associated “‘dimensions”’ with units; and this conception 
of the subject has persisted down to the present day, despite the fact that 

1 Manuscript received in original form September 23, 1949, and as revised, January 26, 1950. 


Contribution from the Massachusetts Institute of Technology, Cambridge, Mass., and 
Brown University, Providence, R.I. 
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change of units is now one of the least important applications of ‘‘dimension” 
theory. 


The modern approach associates ‘‘dimensions’’ with physical concepts rather 
than with units, and one of the principal applications of the subject is the classi- 
fication and correlation of concepts.* Concepts may be classed as primary 
and secondary. We prefer these names to fundamental and derived because of 
the idea of uniqueness that is often associated with the latter terminology. 
In mechanics, length, mass, and time are usually chosen as primary concepts, 
but there is nothing unique about this choice. It is now generally realized 
that the number and selection of primary concepts is largely arbitrary: we 
choose a set of primary concepts that is convenient for the problem in hand. 


Secondary concepts are specified in terms of primary concepts by use of the 
defining equations for the secondary concepts. For example, if /, m, ¢ are 
chosen as primary concepts,** then velocity may be designated as /f~ and force 
as Imt.~* An alternative specification is obtained by writing only the expo- 
nents, as was done by Fourier. Then a concept is designated by writing the 
exponents of /, m, and ¢t in order. The concept of velocity is designated as 
{1, 0, — 1], while force is [1, 1,— 2]. These bracketed quantities are called 
designation symbols and provide a convenient way of specifying concepts. 
Physical concepts should be so chosen that, for a given set of primary concepts, 
there is a 1:1 relation between concept and designation symbol. 


Consider now a physical quantity or physical magnitude. It differs from a 
concept inasmuch as it represents a specific amount, as 4.5 kgm. Evidently a 
physical quantity can be designated in either of two ways: 

(a) By specifying the number of units and the kind of unit (as 4.5 kgm.); 
or, if m.k.s. units are presupposed, 

(b) By specifying the number of units and the concept (as a mass of 4.5). 


For most purposes, (a) is the more convenient since it leaves no doubt 
about the system of units employed. Maxwell and his followers used (a) also 
for ‘dimensional’ specification. The modern treatment of ‘‘dimensions,”’ 
however, employs (4) so that the emphasis throughout can be on concepts, not 
on units. We indicate a specific physical quantity, expressed as in (b), by a 
letter enclosed in curly brackets: 

{A} = Ala’, a’,...a*). (1) 


This quantity is called an idon (From Gk eétéos = class, species) and the 
study of these quantities is called idon theory or idonics. In Equation (1), the 
magnitude of the physical quantity {A} is denoted by A (a pure number). 
The bracketed quantity is the designation symbol, which specifies what con- 


*The operational definition of concepts is treated by P. W. Bridgman in ‘‘The logic of mod- 
ern physics’ (The Macmillan Co., New York, 1927). 

**The advisability of considering two fundamental lengths (il, and 11) instead of one is 
considered elsewhere: 

Moon and Spencer, ‘Modern terminology in physics,” Amer. J. Phys., 16:100. 1948; 
“A modern approach to ‘dimensions’, J. Franklin Institute, 248, 495. 1949. 
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cept is being considered. For instance, a mass of 4.5 kgm. is written {m} = 
4.5 (0, 1, O}. 


A serious attempt was made to employ the traditional notation and names, 
but the semantic difficulties were found to be almost insuperable. How can 
one distinguish, for instance, between the ‘‘dimensions” (length, width, and 
height) of a model and the ‘‘dimensions”’ (length, mass, time), which may be 
referred to in the same sentence? And do the latter ‘“‘dimensions’’ refer to 
the concept, the unit, or the complete physical quantity? In the interest of 
precision, therefore, we were forced to use a new word to represent the one 
particular meaning of ‘‘dimensions’’ that we wished to present. 


Geometry 


The foregoing notation suggests a geometric interpretation. A vector in 
3-space is often specified by writing a number-triple or ordered set [x, y, 2], 
where the three numbers represent the co-ordinates of the vector.* Similarly, 
a vector in k-space is written [a', a?,...a*], or in tensor notation (17), 


a‘ = fa', a?,...a*], Ge A Dee Re 


Comparing with Equation (1), one sees that the idon {A} may be regarded 
as a vector in k-space. The designation symbol [a', a?,... a*] is a base vector 
which fixes the direction of {A}, while the length of the vector {A} is given 
by the scalar quantity A. The numbers a’, a?,...a* are the co-ordinates of 
the base vector, The fictitious space in which the idon vectors lie is called 


idon space. The dimensionality k of this space is equal to the number of pri- 
mary concepts employed in the problem under consideration. 


In kinematics, for example, /ength and time are usually taken as primary 
concepts. Then idon space is two-dimensional with co-ordinate axes / 
and ¢t. Any specific distance is represented by a vector** {/}: 


{7} = 7[1, 0}, 


which is in the direction of the /-axis. A definite value of time is another 
vector, 
ee eR 
while a specific velocity is 
{vy} = vl, — 1). 


Evidently each concept is associated with a definite direction in idon space, 
this direction being specified by the base vector a‘. An increase in the mag- 
nitude of the physical quantity merely stretches the idon vector, whose 
direction remains that of the base vector a‘. 


*See, for instance, L. E. Dickson, Linear algebras (Cambridge Univ. Press, 1930) or 
Algebras and their arithmetics (G. E. Stechert, New York, 1938). 

** Needless to say, the vector nature of the physical quantity in idon space is distinct from 
its possible vector nature in the euclidean space of the phenomena. We assume without loss of 
generality that the physical quantity in euclidean space is a scalar. If it is actually a tensor, we 
deal in idon theory with only the magnitude. 
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Note that with a suitable choice of primary concepts, each base vector has a 
different direction in idon space, and thus each distinct concept is distinguished 
geometrically from all others. Of course this desirable condition will not 
obtain unless the investigator has exercised due care in the selection of his 
primary concepts. The subject is discussed in detail in another paper (J. 
Franklin Institute, Dec. 1949). In the present paper, we assume that such 
an unambiguous selection of primary concepts has been made. 


Idon space is a fictitious space; moreover it does not have all the properties 
of ordinary euclidean space. For example, the ‘length’ of the base vector 
a‘ cannot be obtained by taking the square root of the sum of the squares of 
its co-ordinates a', a?,..a*. By definition, 

{A} = Aa‘ 
and the magnitude of any idon vector {A} is A, irrespective of what numbers 
may appear in the designation symbol [a', a?,...a*] = a‘. Thus a‘ always 
acts as a unit vector. 

This principle may be stated in another way: in general, lengths in different 
directions in idon space cannot be compared. It is evident also that angles 
in idon space have no physical significance. These characteristics lead one to 
suspect that idon space is similar in many respects to affine space,* a fact that 
will be elaborated later. 


Algebra 


The following operations will be postulated: 
(1) Addition. Given two idons {A} and {B}: 
{A} = Aa’, {B} = Bd’. 
Then {A} + {B} is defined only if the two vectors are collinear.** If this 
condition is satisfied, then the sum is defined as 
{A} + {B} = (A+ B)a’. (2) 
Subtraction is defined by 
{A} — {B} = (A — B)a’. 
If A — B = 0, the result of the subtraction signifies zero amount of the con- 


cept denoted by a‘. Addition and subtraction are commutative and associ- 
ative. Like addition, subtraction is defined only for vectors in the same 


direction. 


(2) Multiplication by a pure number (idonic scalar). Given an idon {A} and 
a pure number , then the product is 


k{A} =kAa'. (3) 


*A very readable account of spaces is given in Felix Klein, Elementarmathematik vom 
hoheren Standpunkte aus (Springer, Berlin, 1925, Vol. II); English translation by E. R. Hedrick 
and C. A. Noble (The Macmillan Co., New York, 1939). 

**This corresponds to the well known principle of ‘‘dimensional” homogeneity. 
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Multiplication by a scalar is commutative: 
k{A} = {A} k. 


(3) Multiplication. The multiplication of two idons is effected by multiply- 
ing their magnitudes to obtain the magnitude of the product, and adding the 
co-ordinates of their base vectors to obtain the co-ordinates of the new base 
vector: 


{C} = {A} {B} = AB[a‘ + 5°). (4) 


Multiplication is commutative and associative, and is distributive with respect 
to addition as defined in (1). 


~ 


(4) Division. The division of two idons is effected by taking the quotient 
of the two magnitudes to obtain the new magnitude, and taking the difference 
of the co-ordinates of the base vectors to obtain the co-ordinates of the new 
base vector: 

; hs 


(B}  B {a* — d*). (5) 


These four are the only operations in idon theory. They are called allowabie 
operations. An idon can be raised to an integer power by repeated use of 
Equation (4), and this idea can be extended to a noninteger power y by the 
followin finition: 
ollowing definition (A}” = AT [ya]. (6) 


Of particular interest is the fact that if aa‘ + Bb' + yc' +... + on‘ = 0, the 
product of vectors in idon space is a scalar: 
{A}* {B}*§...{N}” = ASB®.. N’[0,0,...0] 
= A*B®.. N"(0]. 
Geometrically, this condition results when the sum of the base vectors a’, 


b*...n", multiplied by suitable scalars a, 8,..v, is zero. Such quantities, 
1 = {A}* {B}®.. {N}’, are called ‘dimensionless products’ or idonic scalars. 


The determination of idonic scalars is one of the important applications 
of idon theory. For a given physical problem, however, the idonic scalars 
are not unique but may exist in an infinite number of forms. A set of idonic 
scalars is said to be equivalent to any other set that can be obtained from the 
original set by allowable operations (assuming, of course, that the original set 
can again be obtained by allowable operations from the final set). 


Any equation of physics may be written in idon form by merely replac- 
ing each numerical magnitude by a vector in idon space. For example, 
Newton’s relation, 
F = ma 
may be written {F} = {m} {a} 
and since {F} = Fa‘, {m} = mb’, {a} = ac’, 


Fa’ = ma (bt + c’‘). 
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This breaks into the original equation which deals with magnitudes only, 
and the relation among base vectors: 
a’ =-b' +c’, 


The latter equation gives the designation symbol for force. 


Differentiation and Integration 


It is convenient to have rules also for the differentiation and integration 
of idons. By the customary definition of differentiation with respect to an 


arbitrary variable é, d{A} _ , (ey 

d{g} — as—o\ (a8) 
But {A} = A a‘; and since {AA} is the difference of two idon vectors, which 
has significance only if a‘ is the same for both, {AA} = A a‘. Similarly, since 
{&} = £b*, {AE} = AE-D*. 


Therefore {44} a AA 


{Ag} AE 


wi la® — 89 tenn (44) 
ato \ AE 

dA | 

dé 


[a* a, b‘), 


and 


= [a' - 4 
In other words, to differentiate an idon we differentiate the scalar magnitude, 
and subtract the designation symbols. 
Similarly, integration is a process of addition and multiplication. Thus 
Ji{A} até} = [ai + BUS A ae. (8) 
The integration of one idon with respect to another is the integral of the cor- 


responding scalar magnitudes but the resulting designation symbol is the sum 
of the original designation symbols. 


As an example, consider the equation 


d 


In idon form, it becomes 


d 
{F} = dit} ({m} {v}), 
{m} = ma’, 
{vo} = vb* 
{t} tc’, and {F} = Fd’. 


Substitution gives Fd‘ = [a'+ dt - “le (mv). 
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Thus the base vector for force is equal to the sum of the base vectors for 
mass and velocity minus the base vector for time: 

d=a'+b'—c' 

With /, m, t as axes, 

a‘ = (0,1, 0], 6’ = [1,0, — 1], c* = (0,0, 1]. 
Therefore, d* = [1, 1, — 2], 
which is the designation symbol for the concept of force. 


As a second illustration, take the potential difference V between two 


points, x = aandx = bd: 
b 
V -| nie. 
oQ 


The corresponding idon equation is 
ean 
(Vv) -| SF pte. 
2 {Q} 
where {V} = Va', {F} = Fb’, {Q} = Qc' and {x} = xd‘. 
By the rules enunciated above, 


Va' = [b' —c' +d‘ ¥ ls 
a uo 


so a=b'—c+d. 
If the primary concepts are chosen as /, m, t, Q, then 
bi = [1, 1, — 2, 0], c: = [0, 0, 0, 1], d' = [1, 0, 0, 0] 
and therefore the designation symbol for potential difference must be 
a’ = [2,1, — 2, — 1]. 


Of course these same results would have been obtained intuitively from 
theory, but it seems safer to base them on 


” 


the customary ‘‘dimension 
an exact procedure. 


The I] Theorem 


What is the most general relation among idons? Evidently it is the equa- 
tion that is obtained by employing the four allowable operations in all possible 


ways, or 


- K{A}*:{By*;.. {NN}: = 0, (9) 


where K; are scalar coefficients, which may themselves be arbitrary functions 
of idonic scalars. 


Examination of the equations of physics shows that all of them can be written 
in this form. There are a few cases where the foregoing principle seems to be 
violated, but further examination shows that such equations can be rewritten 
in the form of Equation (9). If transcendental functions occur, their arguments 
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are always idonic scalars (‘‘dimensionless’’ quantities). With idons that 
represent physical quantities having ‘‘dimensions,’’ the only function that 
can occur is a power function. Equation (9) therefore represents the most general 
equation of physics. 

The most general equation of physics will now be written in an alternative 
form which often has advantages. Divide Equation (9) by one of its terms. 
The result is of the form, 

> kifA}?. {By'r... {N}%r +1 = 0. (10) 

But addition is only allowable for idons having the same base vectors. Since 

1 is a scalar in idon space, all other terms must also be scalars. Therefore any 
equation of physics can be expressed in terms of idonic scalars: 

Th; = {A}; {By?;... (NY: (11) 

Since the II’s are scalars, any function of them is allowable, and we have as 

an expression of the most general physical equation, 
¢(IIi, Te, ... Um) = 0. (12) 


This is the II theorem, stated in 1896 by Vaschy (20), in 1911 by Riabouchinsky 
(16), and in 1914 by Buckingham (3). 


Derivation of New Equations 


The II theorem is utilized in the derivation of equations where there is in- 
sufficient information to allow a derivation in the ordinary manner. The 
result is usually less complete and less satisfactory than an equation derived 
from a thorough analysis, but it often gives a broad view of the subject that 
can hardly be obtained in am —aier way. 
In using the II theorem, one decides (on physical grounds) on the 7 variables 
that enter the problem: {A}, {B},...{N}. A choice is also made of & pri- 
mary concepts (/, m, t, for instance). Then an idonic scalar may be written 
Tl = {A}* {B}*... {Ny 

= A°B’.. N” [aa' + Bobi +.. + vn'], (13) 

where i = 1,2,..k%. Since II is a scalar, the bracketed quantity must be 
equal to zero, giving the vector equation: 


aa‘ + Bb'+...m' =0. (14) 


Equation (14) is equivalent to a set of k homogeneous algebraic equations 
which the exponents a, 8, ..v must satisfy. The matrix is written immedi- 
ately from the known co-ordinates of the base vectors, a’ = {a',a’,..a*], 
b* = [b', b?,.. b*), etc. and is 
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The rank of this matrix will be denoted* by r._ The only case of physical 
interest is where ry < n; in which case (m — r) of the unknowns (a, B,... v) 
may be given arbitrary values, and the remaining unknowns are uniquely 
determined. Equation (14) is rewritten 

aa’ + Bb'+..+«k* = N'+..4 rn! 


Se a —————" + 


r terms (n — r) terms 
The simplest procedure is to set one of the unknowns on the right (say, \) 
equal to unity, with the remainder of them zero. The solution of Equation 
(14) then evaluates a, 6,. . vy, and Equation (13) gives an idonic scalar, say I. 


Now set vy = 1 and X,..u = 0 to obtain IIx. Evidently there are (” — r) 
of these independent II's, and the final solution is 


pits, Is; 2s. 2 H,>) ==Q: (12a) 


If other arbitrary values had been assigned to X, .. », different idonic scalars 
would usually have been obtained. These II’s, however, would have been 
equivalent to those of Equation (12a). As a matter of fact, the latter set of 
II’s may not be the simplest and most convenient; so it is often advisable to 
transform them by use of the allowable operations and select the set best 


suited to the problem. 


As an example, consider thermionic emission from a plane cathode. On 
the usual basis, using /, m, t, and Q as primary concepts, one writes down the 
pertinent quantities: 


Symbol ‘‘Dimensions”’ 

Current density © rr19Q 
Potential difference 

between electrodes V Pmt?Q"} 
Charge of electrons Q Q 
Mass of electron m m 
Distance between electrodes l l 
Permittivity . 3m PQ? 


Then, according to the II-theorem, 
T= J°V8O'm' ite” 
= (F7671Q)* (PmtQ7)? Orme)’. 
Thus we have four equations, 
| —2a+26+y4 —3r =0 
B+é6 — v=0 
—a — 26 +2» =0 
a—-B+y +27=0. 
Since there are six unknowns, two may be assigned arbitrarily and (if the rank 
of the coefficient matrix is four) there are two II’s. 





*If the n concepts are distinct, each base vector being in a different direction in idon space, 
thenr = k. 
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Let a = 1,» =0. Then 8 = — 3,7 = — 3/2, 6 = 1/2, uw = 3 and 
Im 473 
Re al 


Now let a = 0,y = 1. Wenow find 8 = 1,y = — 1,6 = 0,4 = land 
- Vie | 


Therefore the usual ‘‘dimensional”’ solution is 


J= fie Q*/? 2 r( Ue) 
Q 
Vie 
where F 0 is an unknown function. 


Now consider the same problem, but in idonic notation with /,, l;, m, t, Q 
as primary concepts. The idons are written (See Table IV of Moon and Spencer, 
“‘A modern approach to ‘dimensions’’’, J. Franklin Institute, Dec. 1949.) 

{J} =J[-1, — 1,0, —1,]] 
{V} = V0, 2,1, — 2, — 1] 
{Q} = Q[0, 0, 0, 0, 1) 
{m} = m[(0, 0, 1, 0, 0] 
{2} I[0, 1, 0, 0, 0] 
{fe} =e{—1, — 2, — 1,2, 2]. 
Let {J} = const {V}* {Q}* {m}7 {1}* {e}*. 
Substitution of the above evaluation of the idons leads to an equation for the 
base vectors: 
{[— 1, — 1,0, — 1,1] = [0, 2a, a, — 2a, — a] 
+ [0, 0, 0, 0, 8] 
+ [0, 0, y, 0, 0] 
+ [0, 6, 0, 0, 0) 
+[— 4, — 2u, — w, 2p, 2] 
from which we obtain the five simultaneous equations, 


—-l=-4u 
[-1- 2045-2 
O=aty-u 
| == - 20+ 
1= —a+ 6+ 2u. 
Thus, a = 3/2,8 = 1/2,y = —1/2,6 = —2,n=1 


and 73/2 
= J = const 4/ 2 pL 


which is the usual formulation of Child’s equation. 
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Change of Units 


Another application of idon theory occurs when physical units are changed 
from one system to another (British to c.g.s. units, for instance). Such a 
transformation is equivalent to a change of scales along the co-ordinate axes. 
Consider a k-space with vectors {A}, {B},.. {K} in the directions of the axes: 

{A} = A[1,0,0,..0] 
‘B} = BiO, 1,0,.. 0] 


{K} = K[0,0,0,.. 1]. 
Any other vector is expressible in terms of the foregoing, or* 
{L} = Of{A}* (B}*... {K}* 
= @A°B®.. K*{aa' + Bb'+.. + xk’, 
Cm. tk (15) 


Now suppose that the scales along the axes are changed in the ratios “1, uo, .. 


Uy: : 
, Sal old unit 
cei et te a e 
new unit 
Then in terms of the new units, 
{A’} = um A[1, 0, 0, . . 0] 


{Bi} == 490,12, 0,0] 


{K'} = u,K[0, 0, 0,.. 1]. 
Since the form of Equation (15) must be invariant with respect to a change of 


units, (L’} = 0'{A’}* {B’}8... {Ky 


Qe’ 
“ee. .  ag— {L}> (16) 
8 
Equation (16) allows units to be changed in a systematic manner. Many of 
the older writers on the subject considered this to be the primary use of ‘‘dimen- 
sions.’’ The modern view, however, is that change of units is a rather trivial 
application, particularly since any change of units is unproductive, a complete 
waste of time which can be eliminated by the universal use of the m.k.s. system 


(12). 
Transformation of Co-ordinate Axes 


The old idea that /, m, and t, are sacrosanct has been abandoned (1). Numer- 
ous choices are possible. Select a set of k primary concepts, which corresponds 
to a choice of rectilinear co-ordinate axes in idon k-space. Any idon is then 


represented by the vector : 
: {A} = Aa‘. 

*Here © is a scalar constant, which allows the use of peculiar units for} LL‘. For example, 

if axes are l, m, t, and } L{ represents energy, then m.k.s. unit of } L{ is obtained directly from the 

three fundamental units and is m*kg sec * = joule. In this case, Q9=1. But if L is expressed 


in calories, 8 = 1/4.182. 
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Now make a different choice of primary concepts, keeping the number equal 
to k as before. Then the above idon is written 
{A’} = Aa”. 

The use of the same base letters (A and a) indicates that the idon vector is 
unaltered (18), while the primes show that the co-ordinate axes have changed. 
Geometrically, we have a given vector {A} which remains fixed as the co- 
ordinate axes are rotated to new positions. The co-ordinates of the base 
vector a” are different from those of a', but both are unit vectors. The length 4 
A is identical in both cases and need not be primed. 











It is well known (18) that the transformation from one set of rectilinear 
axes to another is always accomplished by a linear homogeneous transforma- 










tion, so 2 uw 
(a = Bia’ + B, a +...+B, a® 
9 * 1 of 2 or | 
a” = Bi’ a' + BE a’ +... + Bi at | 
a” = BY a' + BY a? +... + Brat. | 
In index notation (18), this set of equations is written | 
a” = Bi a‘, G,#’ = 1,2,...), (17) 
which is the equation of an affine transformation of co-ordinate axes. Conse- 
quently, idon space has affine properties. 


For instance, consider the transformation from /, m, t axes to l, P, t axes. 
In terms of /, m, t, we have idons of the form 





{A} = Aa’, or specifically, 


{1} = 1[1,0, 0], {P} = P[2,1, — 3], {t} = ¢[0,0, 1]. 
In terms of /, P, f, (4") = Ae”, or 


(1) = I[1,0, 0], {P’} = PlO, 1, 0}, {e’} = 40, 0, 1). 
Substitution in Equation (14) evaluates the co-ordinates of Ai, which are 
l —2 0 
0 1 0 
0 ee 
Any quantity expressed in terms of /, m, t can then be transformed to 
the /, P, t, system by means of Equation (14). 

Energy, for example, has the base vector a‘ = [2, 1, —2] in the /, m, t 
system. Substitution into Equation (14) leads to the new specification, 


a” = (0,1, 1]. 
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COMPARATIVE HEAT LOSSES FROM BARE AND ANODIZED 
ALUMINUM SURFACES TO STILL AIR! 


By A. C. Burk AND RoBerT H. Hay 


Abstract 


The effect of anodic films on the total heat losses to still air from the surfaces 
of cylindrical aluminum containers maintained at temperatures between 23°C. 
and 100°C. has been studied by recording the power required to maintain the 
containers at various temperatures above ambient. Anodic coats of the order 
of 20u thick required increases in power of 43°% to 60% at 30°C. and 68% to 
79% at 90°C., the variations at any one temperature being due to variations in 
thickness and color of the coating. The emissivity coefficient for the clear 
anodic coats used is estimated from the data to be about 0.80. 


Introduction and Theory 


Recent consideration of the heat loss problems involved in the design of 
certain aluminum tanks revealed the lack of reliable heat-loss information in 
the literature for aluminum surfaces at temperatures between 20°C. and 100°C.., 
particularly if the aluminum has been anodized or anodized and dyed black. 
It was therefore decided to perform simple experiments which might lead to 
reasonably reliable information on such heat losses. 


A body in air and at a temperature higher than that of its surroundings will 


lose heat to them by conduction, convection, and radiation. In the first two 
cases, particularly for the conditions of this experiment, the heat loss will be 
at a rate approximately proportional to the temperature difference between 
the body and surrounding air. In the last case, however, rate of heat loss is 
proportional to the fourth power of the temperature, i.e., 


Q « «(T, — T3), (1) 


where ¢ is the emissivity of the surface and 7, and T, are respectively the 
absolute temperatures of the radiating body and of its surroundings. Heat 
transfer by radiation is, of course, unaffected by the presence or absence of air. 
As the quantity proportional to the amount of heat lost in unit time by the 
surface by all three of the above mentioned modes of heat transfer, the writers 
chose the electrical power used by an immersion heater to hold the water in 
an aluminum tank at a given temperature. Different surface treatments can 
be applied to any one tank and the different power inputs to the heater will 
be a measure of the ability of the various surfaces to dissipate heat. 


Some measure of the emissivity of anodized surfaces can be gained by com- 
paring the power input for an anodic coat with that for a similar coat dyed 
black on the same tank. Since the nature of the surfaces in the two cases would 
be almost exactly similar (i.e., conduction plus convection would be equal for the 


1 Manuscript received January 20, 1950. 
Contribution from the Physics Division, Aluminium Laboratories, Limited, Kingston, Ont. 
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two cases) the difference in power would be due to a difference in radiation alone. 
Using the subscript a for the anodized tank and 6 for the anodized tank dyed 
black and somewhat arbitrarily taking the emissivity for the dyed tank to be 
90% (the emissivity of a “black body” is 100%), the total heat loss per second 
in watts for each of the two cases would be, 
Q. ecA(T: a. T?) + q; (2) 
and Q, = 0.9cA(T; — T?) +4, (3) 
where ¢ = 5.74 X 107” watts cm.~ deg.~! (Stefan’s constant), 
and = coefficient of emissivity (€ = 1 for a perfect black body), 


radiating area in square centimeters (taken as surface area of the tank), 


€ 

A 

T, = temperature of radiating surface in degrees absolute, 
Ts 


= temperature of surroundings in degrees absolute, 
g, = aterm to represent heat loss by conduction plus convection. 


g, can be calculated from Equation (3) and substituted in (2). The value for e¢ 
can thus be calculated from observed heat losses and temperature differences. 


Experimental 


Two cylindrical tanks were made up, each nominally 4 in. in diameter and 
9 in. high. One tank was made from Alcan 57S* alloy sheet 0.040 in. thick and 
the other of Alcan 2S* alloy sheet 0.062 in. thick. The top end of each cylinder 
was closed by a flat plate which carried three diametrically located holes over 
which were brazed “chimneys”, or 3-in. long cylinders, through which a knife- 
type immersion heater, a stirring mechanism, and a thermometer and thermo- 
couple lead wires were inserted. The amount of water surface exposed to the 
air was thus kept at a minimum. All anodizing was done in 15% (by weight) 
sulphuric acid at 21°C. with a current density of 18 amp. per sq. ft. and sealing 
by immersion in boiling water. Coatings were stripped after each test by im- 
mersion for 10 to 15 min. in a phosphoric— chromic acid bath** at 95-100°C. 
Black dyeing was accomplished by immersing the coated but unsealed tanks 
for 10 min. in a solution of Nigrosene Black—12525 (National Aniline Co.) 
made up to a 12 gm. per liter concentration and a pH of 6.1 and maintained 
at 65°C. For test, a tank was filled with distilled water to within 1/2 in. to 
3/4 in. of the top of the chimney tubes. The knife-type immersion heater was 
fed by a variac transformer which in turn was fed from a Sola constant voltage 
transformer on the 60-cycle main. Power in watts was calculated from voltage 
and current readings on Westinghouse type Py-5 one-half per cent meters. 
The center opening admitted a mercury thermometer to obtain readings of 
water temperature and a thermocouple connected to a sensitive potentiometer 
to read small temperature differences when approaching the equilibrium tem- 
perature of the water. A second mercury thermometer was used to read the 

* Aluminum Company of Canada designations. Corresponding Alcoa designations are 52S 


and 2S respectively. 
** ASTM Designation B 137-45. 
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temperature of the surrounding air. A propeller type stirrer was used to agitate 
the water. An initial test was run to see if there was any heating effect due to 
the agitation of the stirrer, but no rise in water temperature could be noted 
from this cause in several hours’ operation. 


The procedure of the experiment was to obtain the equilibrium temperature 
for a given amount of power input to the heater. Power steps were adjusted to 
give appropriate temperature intervals. The time necessary to reach equili- 
brium for any one steady power input was from two to four hours, but the 
experiment needed only occasional adjustment and readings taken. When the 
water temperature reached equilibrium for any one power this temperature 
was reduced or increased for plotting by the amount that the room temperature 
was above or below 23.0°C. respectively, as 23.0°C. was the average room 
temperature. The amount of this correction was never more than 2° and since 
most of the heat transfer is by conduction and convection (which are roughly 
proportional to the first power of the temperature difference) the resulting 
error would be small. 

Results 


Fig. 1 gives the power necessary to hold the water in a 57S tank at various 
temperatures for different surface finishes while Fig. 2 gives the same data for 
a 2S tank. All figures are for a room temperature of 23.0°C. The two tanks used 
in these tests are of the same shape and of approximately the same size but 
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Fic. 1. Heat loss from an aluminum water tank to its surroundings as measured by power 
necessary to maintain equilibrium temperature above ambient; Alloy 57S bare and anodized; 
ambient temperature, 23°C. 
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they are not nearly enough identical to permit comparison to be made of losses 
from surface finishes on different alloys. 


In Fig. 3 is shown the increase in power input in per cent for tanks anodized 
on the outside over the power input for the same tanks bare and polished on 
the outside. The power figures for Fig. 3 are taken from curves in Figs. 1 and 2. 
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Fic. 2. Heat loss from an aluminum water tank to its surroundings as measured by power 
necessary to maintain equilibrium temperature above ambient; Alloy 2S bare, anodized, anodized 
and dyed black; ambient temperature, 23°C. 





On the assumption that an emissivity value of 90% can be given to the 2S 
tank anodized 40 min. and dyed black, then values for the 40 min. anodized 
undyed finish were calculated and are given in Table I. The average e of these 











TABLE I 


EMISSIVIVITY IN ©) FOR 2S ANODIZED 40 MIN. WHEN 2S ANODIZED 40 MIN. 
PLUS BLACK DYE IS TAKEN AS 90% 














Heat loss, watts % Emissivity 



















Water ee es Bi eels So ay 

temp., ; ; aa : ean) , 
°c. | Anodized | Anodized 40 min. | Anodized 40 min. | Anodized 

40 min. plus black dye | plus black dye 40 min. 

90.0 |; 88.9 93.3 90.0 81.1 
80.0 72.5 74.7 90.0 84.5 
70.0 56.5 57.9 90.0 85.5 
60.0 42.0 43.2 90.0 84.9 





50.0 28. 29.6 90.0 85.1 
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values is 84.2%. A 40 min. anodic coat on 2S with conditions as given above 
would be of the order of 15 to 20 in thickness, so that a figure of 84% for 
emissivity is in fair agreement with other published figures for anodic films (3), 
provided the published curves are extrapolated slightly. Additional valuable 
reference material is contained in papers by Draper (1) and Lux (2). 


.G - 87S ANODIZED 90 MINUTES 
© - 2S ANODIZED 40 MIN.& DYED BLACK 
4&~ 57S ANODIZED 40 MINUTES 
¥-2S ANODIZED 40 MINUTES 





Fic. 3. Curves for anodized tanks showing extra power input required to maintain given water 
temperature as percentage increase over power necessary to maintain same temperature in bare tank. 


In the formulas given above for calculating emissivity, 7, is the temperature 
of the radiating surface. This temperature was not measured, the water temper- 
ature being used instead. It would be slightly higher than the actual temper- 
ature of the radiating surface. However, a 5-degree error in T, at 90°C. has 
approximately 1% effect on the value of ¢ while at 50°C. a 2.5-degree change 
in 7, alters « by about 0.7%. Since calculations showed that the difference 
between the water temperature and the temperature of the radiating surface 
would be of the order of 0.5°C. and not more than 1°C., the error arising from 
the assumption that 7, equals water temperature has been neglected. 


Discussion of Results 


The curves of Figs. 1 and 2 give all the data obtained. From them have been 
constructed the curves of Fig. 3. These show clearly the most interesting con- 
clusions to be drawn from the experiment. 
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The agreement between values obtained for the two runs with the bare 
surfaces is very good, the second run being made after anodic coatings had 
been stripped off and the bare aluminum repolished. From Fig. 1 (or Fig. 2) 
we see that the heat losses for the various thicknesses of anodic coating (thick- 
ness being taken, to a first-order approximation, to vary directly with anod- 
izing time) are of the same order of magnitude and that they all differ very 
much from those for a bare, polished aluminum surface. Even as thin a coat 
as is obtained with five minutes anodizing time has a marked effect on heat loss. 
It is clear, also, that though in general heat loss depends on some factor related 
to the time of anodizing and also is: greater for black than for clear surfaces, 
the differences we have observed due to these factors are small. 


From Fig. 3, referring to the 57S alloy curves, it is seen that the heat loss is 
distinctly higher for a 90 min. coating than for the thinner 40 min. coating, 
and that the two curves show a greater difference in heat loss at the lower 
equilibrium temperatures than at the higher, even though the loss difference 
is not very great. At first glance there appear to be consistent differences 
between the heat loss for the 57S alloy tank and for the 2S alloy tank. It is 
well known that each alloy has its own anodizing characteristics, and it may 
be that a difference in heat loss could be attributed to this difference in anodic 
film character. However, the two tanks used are sufficiently different in size 
to rule out drawing any such conclusion from the data presented here. 





Applications may occur for which the temperatures are well above 100°C. 
It would be unwise to extrapolate the curves given here for more than 10° or 
20°C., and further experimental work would have to be done in the region of 
interest if it lay outside such an extrapolation range. 
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A 60 db. NONLINEAR AMPLIFIER' 
By J. A. CARRUTHERS 


Abstract 


A nonlinear amplifier is described which has been used to drive an Esterline 
Angus 0-1 ma. recording galvanometer as part of the automatic recording 
apparatus set up for taking microwave antenna patterns. The amplifier output 
is approximately proportional to the Jogarithm of the input voltage, and permits 
plotting a 60 db. variation in input signal. The speed of response is limited by 
the galvanometer movement but is adequate for obtaining a pattern through 
100 degrees of azimuth in about two minutes. The plotting accuracy is better 
than 1 db. if care is taken to check calibration frequently. 


Introduction 


In the investigation of fundamental properties of microwave diffraction, as 
well as in the routine measurement of antenna directivity characteristics, it is 
necessary to measure the power transmitted between the system under test 
and a distant point, as a function of the angle of rotation of the test apparatus. 
In order to eliminate much of the time required for a point-by-point plot of 
such patterns an inexpensive automatic recording system has been built which 
permits obtaining the complete information in a few minutes. 


The essential parts of the apparatus are an Esterline Angus 0-1 ma. recording 
galvanometer and a distorting amplifier which is described here. The amplifier 
must be nonlinear since it is desirable to plot up to a 30 db. range of radio- 
frequency power picked up at the receiver; because a square detector (bolo- 
meter or crystal) is used at the receiver, it follows that a 60 db. variation in 
the signal from the detector must be recorded. To cover this plotting range, 
an approximately logarithmic characteristic has been obtained by the use of 
automatic volume control (A.V.C.) bias on supercontrol tubes whose gain is 
thereby decreased as the signal level increases. Amplification takes place at 
600 cycles, which is the frequency used to amplitude modulate the R.F. carrier 
and which therefore appears at the output of the detector. The 600 cycle 
voltage is converted to d.c. to actuate the galvanometer. 


Since an amplifier cannot be truly logarithmic starting from zero output, 
it is designed to have a characteristic which is lin-log, i.e., linear for small 
signals up to a crossover point beyond which the output increases in proportion 
to the logarithm of the input. This is expressed by the following relations: 


Ey =k E; E; <E, (1) 
zs) Bi> i (2) 


+ 


“Cc 


Eo =k E.(1 + log, 


1 Manuscript received December 23, 1949. 
Contribution from the Department of Physics, McGill University, Montreal, Que. 
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where Eo= output voltage, 
E;= input voltage, 
E.= crossover voltage which is a constant of the design of the am- 
plifier. 
The constant, k, is related to the voltage amplification of the amplifier 
as follows:— 


From (1) V.A. = =e k E;< E, (3) 


Ly 


4% 


ee E E Ss Seite 
(2) VA. =— = k= fe Eo &, (4) 
OF E; E 
If it is considered that E, is the lower extreme of the plotting range, then for 
a 60 db. increase in signal E;(max.) = 1000 E,. Substitution in Equation (4) 
shows that the A.V.C. bias must act to decrease the gain of the amplifier by a 


factor of 130 in order to cover the desired plotting range. 


Description of the Amplifier 
Fig. 1 shows the circuit diagram of the amplifier which is described below 
in detail. 


(1) The Input Stage 
The cathode following V; was found necessary to isolate the gain controls 


from the main amplifier. Any change in the impedance of the grid circuit of 
V2 brought about a change in calibration, an effect which was eliminated by 
the cathode follower. The controls permit a 12 db. variation in gain in 3 db. 
steps as well as a 3 db. continuous adjustment. The amplifier has been used 
either as a termination on a 500 ohm line or in parallel with a terminated line. 
For this reason the input impedance can be altered from 500 ohms to high 


impedance. 


(2) The Main Amplifier 

The two amplifier stages V2 and V; have A.V.C. bias applied to the grids 
in order to decrease the gain with increasing signal. V; is an A.V.C. amplifier 
which is isolated from the A.V.C. diode (Ve_1) by a cathode follower stage 
( V5—~2). The signal to be measured is coupled to the detector Vs_2 through V5_1, 
which is a triode amplifier whose gain has been adjusted to give the most 
desirable shape to the cramping characteristic of the amplifier, without af- 
fecting the feed-back loop. This is particularly helpful since it allows a control 
on the crossover point without danger of setting the amplifier into oscillation. 

The time constants have been adjusted to allow as rapid a response as 
possible while maintaining stability. It requires approximately one second for 
the amplifier to reach a sensibly steady state for a 60 db. change in the input 
signal. The 2.2 megohm resistor, and the 0.05 and 0.01yu/. condensers between 
V; and V4 were found to decrease the tendency to oscillate by reducing the 
low frequency gain without introducing large phase shifts at the low frequencies. 
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(3) The Metering Circuit 

The remaining stages make up the detector and balanced voltmeter circuits 
which drive the d-c. meter. The output impedance of the voltmeter stage has 
been adjusted to give optimum speed of response of the galvanometer without 
overshooting. The 5K potentiometer between the plates of V7; and V7_2 is 
for setting the meter to zero at no signal. 


Operation and Performance 


Fig. 2 is a plot of relative meter deflection versus change in input voltage, 
expressed in decibels. Allowing for the square law detector ahead of the ampli- 
fier it is seen that a 30 db. plotting range is feasible. Full scale deflection 
requires an input of 0.05 v. 
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Fic. 2. Over-all calibration curve for amplifier and recording galvanometer. 


The calibration has been found to vary somewhat with tube temperature 
and age, but a half-hour warm-up period is sufficient to ensure a precision of 
1 db. after a thorough initial break-in of all tubes. For maximum precision, a 
500 ohm attenuator box has been used in the line feeding the amplifier so that 
calibration can be checked very rapidly: after completing a pattern the antenna 
is rotated to the angle of maximum signal and the attenuation is varied in 
10 db. steps to give six points on the calibration curve. 


The response time for large changes in signal is less than two seconds, which 
allows taking a normal horn pattern through 100° of azimuth in about two 
minutes. The accuracy in following small changes in intensity is limited by 
the friction between the pen and the chart. This friction was decreased con- 
siderably by rebalancing the galvanometer needle so as to exert minimum 
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pressure on the paper, although this requires keeping the inkwell filled to the 
same level at all times in order to maintain constant buoyancy effect from the 
part of the needle in the ink. 


Fig. 3 shows a comparison of the curves obtained for a horn pattern when 
replotted from the recorder chart, and when obtained by a point-by-point plot 
(hand pattern). A plotting accuracy of about 1 db. is indicated. Fig. 4 shows 


E-PLANE PATTERN AN 18 —— PLOTTED FROM 
HORN SLANT LENGTH 99.5CM. / RECORDER PATTERN 
FLARE ANGLE 18°28’ i ——---POINT- BY- POINT PLOT 
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Fic. 3. Comparison of replotted recorder trace and hand pattern. 
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Fic. 4. Section of recorder chart showing the actual trace. 
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the trace of the recorder chart from which the graph in Fig. 3 was obtained. 
Although this trace is distorted by the curved arc of the pen and the departure 
from logarithmic response of the amplifier, it is often necessary only to observe 
the effects of varying some part of the antenna system, and replotting the 
trace is not required. 
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DEFLECTIONS OF AN INFINITE PLATE! 


By Max WyMAN 


Abstract 


Problems associated with the thickness of ice sheets on Canadian lakes led the 
author to investigate mathematically the deflection of a loaded plate resting 
upon an elastic foundation. Certain simplifying assumptions, which appear to be 
not unreasonable, permit the solutions obtained to be applied to the strength of 
a floating ice sheet. Relations for maximum stress and maximum deflections in 
terms of known functions are derived, both for a concentrated load and a uni- 
form load distributed over a circular area. The resulting expressions should be 
tested experimentally to determine their domain of validity. Owing to similarity 
of conditions, these results may also be applied to the design of concrete road- 
ways and airport runways. 


Introduction 


During the last war considerable activity took place in the northern part 
of Canada. The problem of transporting heavy equipment over frozen lakes 
became important and Prof. I. F. Morrison of the University of Alberta 
received several requests for a formula which would give the thickness of ice 
required to support loads of various proportions. Engineers interested in the 
problem had no accurate formulas and were using rough charts whose validity 
seemed doubtful. In addition to the requests made to Professor Morrison 
several other members of the Faculty of Engineering were asked to consider 
the problem of landing aircraft on frozen lakes. Here again the problem was 
to determine the thickness of ice in order to make a safe landing. 


A rigorous mathematical solution of the problems enunciated would be very 
complicated. Some of the factors which would require consideration are 
that ice is a nonisotropic elastic medium and that its elastic properties show 
wide variability under different conditions. In addition the problems of trans- 
portation involve moving loads which are of course much more difficult to 
treat than the corresponding problem involving only stationary loads. 


In order to make a start toward the solution of the above problem we shall 
make several simplifying assumptions which will enable us to carry out a 
mathematical analysis. The final formulas which we shall obtain should be 
tested experimentally in order to obtain their domain of validity. 


As the ice deflects into the water there is a restoring force due to the water 
pressure which is proportional to the deflection. For this reason the mathe- 
matical problem to be solved resolves itself into determining the deflection 
of a loaded plate resting on an elastic foundation. Since we are considering a 
frozen lake there is no serious error in considering the ice plate to be of infinite 
extent. More serious is the assumption that we shall make that ice is an iso- 


1 Manuscript received in original form August 13, 1949, and, as revised, January 10, 1950. 
Contribution from the University of Alberta, Edmonton, Alberta. 
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tropic elastic medium. In order to overcome some of the error introduced by 
this assumption we shall be conservative in choosing the values of our elastic 
constants and in this way we hope that the numerical values which we shall 
obtain for the safe thickness of ice will be on the conservative side. Finally 
we shall deal only with stationary loads and shall assume that these loads are 
big enough to make the weight of the ice itself negligible. 


On the basis of the above assumptions the problem has been considered to 
some extent by others (3, pp. 275-282). In most cases the final form of solution 
has been left as an infinite series which does not lend itself to easy numerical 
calculations. For this reason we shall treat the problem from the beginning 
and shall be able to show that the final form of solution can be expressed in 
terms of known functions for which tables of values exist. 


The fundamental differential equations, upon which our mathematical prob- 
lem is based, are all readily available in Timoshenko’s book on Plates and 
Shells. For this reason we shall try, as far as possible, to restrict our notation 
to that used in this book. 


Section 1—Concentrated Load 


In Fig. 1 we consider an infinite isotropic plate of thickness / resting on a 
liquid of density p. If we let w be the deflection of the middle plane of the 
plate then it is known (3, p. 58) that w must satisfy the differential equation. 


324 es ee dw 
Las apia (dyn a ey 
r dr \ drbr adr dr D 
where g is the intensity of the external load and D is the flexural rigidity of 
the plate. For a concentrated load P the intensity g at every point, except the 


origin, is entirely due to the pressure of the liquid on the plate. Its magnitude 
is equal to — pgw. Thus the differential equation to be solved is 


oe ae i of 2) )} pew 
5 ak ee mel eee "= 0. 
r dr \ dr Lr dr dr sr D 


Pp 


w 


Fic. 1. 


If we let r = Ix, where / = (D/pg)”', we find that (1.2) can be put into the 
form 
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(z ne <) (fe + 1 oe) 4 ~ (1.3) 
act «dx dx? x dx 


In order to find the general solution of (1.3) we note that any solution of 





te 1 dw 
OO SOU can Os (1.4) 
dx? x dx 
will satisfy (1.3) providing a = + ft aed, However, for these values 
of a it is well known that the Bessel Functions 


Jo(i*/2x), Jo(( — i)**x), Ko(i'/2x), Ko(( — i)”*x) 
are all solutions of (1.4). Since 
J (23 2) = ber x +i bei x 
Jo((— 7)*2x) = ber x — i bei x 
Ko(t! 2y) = ker x +ikeix 
K.((— 1)"2x) = ker x — i kei x 
we see that the modified Bessel Functions ber x, bei x, ker x, kei x provide 


four linearly independent solutions of (1.3) and hence the general solution of 
(1.3) is r , 
w = A, berx + Az beix + Az kerx + Ag keix, (1.5) 
where Ai, Ao, A3, A4 are arbitrary constants. This gives the general solution 
of (1.2) to be 


w = Ay ber (r/l) + Az bei (r/l) + Ag ker (r/l) + Agkei (r/l). (1.6) 


The modified Bessel Functions in (1.6) are discussed in great detail in 
McLachlan’s book, Bessel Functions for Engineers (2). Many of the formulas 
involving these functions are given on pages 168-172 of this book and numeri- 
cal tables are given on pages 177-178. Since we shall be using these formulas 
and tables throughout the rest of the paper we shall refer to it by the abbrevi- 
ation M in the remaining part of this paper. 


From M, page 169, formulas 185, 186 we see that both ber x and bei x 
increase without limitasx— ©. Similarly formula 203 shows that ker x — © 
as x +0. Since our displacement w must be finite at r = 0 and r = © we 
see that A; = Az = A; = 0. Hence 


w= A, kei (r/l). (1.7) 


From the fact that our plate is in equilibrium we can conclude that the total 
restoring force of the water must equal the concentrated load P. Hence 


| pgw 2xr dr = P. (1.8) 


0 


Thus, 
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2rpes| r kei (r/l) dr = P. (1.9) 
From M, p. 172, formula (235), this integral can be evaluated to give 
— 2rpgAgl E ker(r/t) | = FP. (1.10) 
Since M, p. 172, formula (232) implies lim [r ker’(r/l)] = 0 and M,p. 170, 
formula (203) gives lim [r ker’ (r//)] oe we find that 
fie Ay = — P/2mpgi’. (1.11) 


Thus the deflection, at any point, due to a concentrated load P is given by 
w= — P kei (r/l)/2xpgi?. (1.12) 


This immediately implies that the maximum deflection is at the origin and is 


riven by faa 
' Wmaz = P/8 pgl’. 


Section 2—Arbitrary Loading 


In Fig. 2 we are considering the top of the plate and shall assume that the 
plate has been loaded in an arbitrary manner over the region enclosed by a 
closed curve. We desire to find the deflection at an arbitrary point Q due 
to this load. It should be noted that although Q is shown as an external point 
to the area of loading, the following derivation also applies if it is an interior 


pe int. 
¥ 
oa 
Zek 7 
Fic. 2. 


Let us take polar co-ordinates with respect to the line OQ with O as the polar 
origin. Let S be any interior point with polar co-ordinates (r’, 0). Further 
suppose that SQ = R and that the density of loading at S is given by some 
function g(r’, @). Thus the concentrated load at S is given by q(r’,0)r’ dr’ dé 
and hence the deflection dw at Q due to this load is given by (1.12) to be 


dw = — q(r’, 0) kei (R/I) r'dr'd0/2mpgi*, (2.1) 
where R? = r’2 + r? — 2rr’ cos 6. Thus the total deflection can be obtained 
by 
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w= | aw 6)r’ kei (R/1)dr'd0/2rpgi?, (2.2) 
A 


where A is the region of loading. 


Although (2.2) gives the general expression for the deflection due to any 
type of loading, it will be difficult to evaluate the integrals involved in most 
cases. For complicated boundaries it will be necessary to use numerical 
integration. However, in a few cases it is possible to carry out the 
integration in terms of known functions. We shall show in a later section 
that (2.2) can be used to find an exact expression for the deflection due to a 
uniform loading over a circular region. Before doing this however it is neces- 
sary to obtain the so-called addition theorem for the function kei x and an 
identity satisfied by the modified Bessel Functions. 


Section 3—The Addition Theorem for kei « and an Identity 


It is quite possible that the addition theorem for kei x and the identity, 
involving the modified Bessel Functions, which we shall prove, are already 
known. However we have not been able to find them in print and, since the 
proofs are short, we felt it advisable to include the proofs in this paper. 


Addition Theorem 
It is well known (4, p. 361) that the function Ko(x) satisfies the addition 


theorem 


Kola) = > KelZ) Tals) cosimd, (3.1) 


m= — © 
where x? = Z? + 2? — 2Zz cos 6. Moreover the expression (3.1) is valid as 
long as |z|< |Z|. Replacing x, Z, z in the above expression by ity, i#Z, ide 
respectively we have 
CO 


Kf(#*x) = > Knit 2Z) Im (i'/?z) cos mé, (3.2) 


m=— 


where again, x? = Z*? + 2° — 2Zz cos #6. From M, formulas 201, 205 and 163 


we have 


Ko(i'’*x) = kerx +i keix (3.3) 
K,(1'/?Z) = i™(kerm Z + 1 kein Z) (3.4) 
In(i'/?z) = i-™(berm z + 1 beim 2). (3.5) 
Thus 
kerx+tikeix= > (ker,Z +72keinZ) (berms +72 beinz) cosmé (3.6) 
n=— Cc 


If we assume x, Z, z are all real, we can, by equating imaginary parts, find that 


eo 


keix = > (keinZ bernz + kernZ beinz) cos mé. (3.7) 


m=— © 


Expression (3.7) is the so-called addition theorem for kei x and will hold as 
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Be |g! * . 
long as |z|< |Z}. From the symmetry of the expression x? = Z? + 2*— 2Zz 


cos 8, it is obvious that 
Cc 


keix = >) (keimz bermZ + kermz beimZ) cos m0 " (38) 


m=— 0 


0 


Since | cos modd0@=0,m=2+1,+2,+3,... 


=e,  i-—/0 
we have that 


| kei xd 0 = 2m (kei Z ber z + ker Z beiz), if |z|< [Z| and — (3.9) 
0 


“ 


| ke xd 0 = Qn (kei z ber Z + ker z bei Z), if |Z|< |e/. (3.10) 


~~ 
0 


We shall make use of these relationships in evaluating the integral (2.2). 


An Identity 


From M, page 156, problem 60, we have 


. 1 cae, 
=| Kole =a Io(z) — Io(z) ame Kote) = 1. 
dz dz 


. “1/2 
Replacing z by i'/*z we have 


=| Ko(i"s) “ Io(i'?2) — Io(t’/*2) ‘ Kui") | = 1. 
Zz az 


From Relations (3.3) and (3.5) of the present paper we find 


| (ker z+ikeiz) (2 bers +72 & bei :) 


az az 


(ber 2 + 7 bei 2 (= ker 2 -+12 zs kei :) = 1. 


az az 


Assuming z is real, we can equate real and imaginary parts and find that 


d ~ 2a ; d conde ; 
Zz ker . ber z — kei z = bei z — ber z— ker z + bei z—keiz]| = 1 
dz dz dz dz 


and (3.12) 


se i ‘ itl d : 
kei g — ber z + ker s— bei z — bei g — ker zs — ber z — keiz 


= (). 
dz dz dz dz 


(3.13) 
Denoting differentiation by the prime notation these identities may be written 


z (ker z ber’s — kei zg bei’s — ber g ker’z + bei g kei’s] = 1 


(3.14) 
(3.15) 


kei z ber’ s + ker zg bei’ zs — bei z ker’ s — berz kei’ z = 0 


Here again we shall make use of these identities in the next section. 
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Section 4—Uniform Circular Loading 


We shall assume that the loading is uniform over a circular area of radius a. 
Since the loading is uniform g(r’, #) is a constant which we abbreviate to q. 
Thus (2.2) gives us the deflection to be 


a f2r 
ok ay ee (| |r’ kei (R/I)dr’ is). (4.1) 
2rpgi* oJo 









In evaluating this integral we distinguish two different cases. 
g g § 










External Point — r' <r. 





Case | 


For this case (3.9) applies and 








a 

w= — a_| r’ [kei (r /l) ber (r’ /l) + ker (r/1) bei (r’ /1)] dr’. 
pgl? 0 

Evaluating the integrals by means of M, p. 170, formulas 189, 190 we have 






sg wi bisa? (hy er le 13) hae) Week (4.2) 
pg 









l. 





where b = a 


Thus (4.2) gives the deflection at any external point. 






Case 2. Internal Point — OS rX& a. 






As before, the deflection at any interior point is given by 






—1_f {kei (R/I) dr’ dé. 
Qmpgl?* a” 


w= 





This can be split up into two integrals as 





| 


2 
| 
| 


q * | | 2 yr’ kei (R 1)dr’ dé +| [ yr’ kei (R syardo (4.4) 


2rpgi* 0/0 0 





In the first integral (3.9) applies and for the second (3.10) may be used. Thus 





fr 
w= — 2 4] r’ {kei (r /1) ber (r’ /1) + ker (r/1) bei (r’ /1)] dr’ 
pgl® Jo 





ote | r’ [kei (r /l) ber (7/1) + ker (r’ /1) bei (r /1)] wt. 





These integrals may be evaluated as before and the terms rearranged to give 





w= a ker (r /l) ber’ (r /l) — bei’ (7/1) kei (r /l) + kei’ (r /l) bei (r /L) 
pg 
— ker’ (r/l) ber (r/l)] + Ci! (b)ber (r/l) — kei’ (b)bei (r /1)]. 
pg 





Using (3.14) we find 
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b 3 ‘ : Be 
welt +f [ker’ (b)ber (r /l) — kei’ (b)bei (r /1)] (4.5) 
PE = pk 
and that (4.5) will give the deflection at any interior point. The maximum 
deflection occurs at r = 0 and is given by 


Wmax = g(1 + b ker’ (b)) /pg. (4.6) 


Section 5—Application to Safe Thickness of Ice 


Once the deflection function is known then it is possible to determine the 
stress distribution in the plate. From these it should be possible, from the 
rupture theories available, to determine the thickness at which the ice plate 
will break. We shall assume that the ice plate will break when the tensile 
stress reaches a certain maximum. 


Using Timoshenko’s notation we let 17, and W/; denote the bending moments 
per unit length and co, the tensile stress. The expressions given by him are 
dw v dw 
M,=- D( + — “*) 


dr? r dr 


M, rai yee) 


dr? 


(M,) 


(Gr) max = — 6-—— max. 
h? 
It is not difficult from (4.5) to show that (M,)max occurs at r = 0 and at this 


point M, = M;. Calculating from (4.5) we find 
3qbD (1 + v) kei’ b 
(Or)max = ee . (5.1) 
Since the load is uniform the total load P over the circular area is given by 
P =rna*q. Thus 
oO = 3PbD (I +») kei’ b . 
\“ ry} max ra? pgh?l? 
Since D/pg = I, and b = a/l 


ve eo (5.2) 
3(1 + v) kei’d 


Formula (5.2) thus relates the maximum load P to the thickness h of ice. 
It should be pointed out that the flexural rigidity D is given by 
D = Eh®/12(1 — v*) 


and hence b will also depend on h. 


NUMERICAL ANALYSIS 
It was pointed out that the elastic properties of ice show wide variability 
under different conditions. For this reason the values we choose for the elastic 
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constants are definitely on the conservative side. From work done by Hess 
(1, pp. 49, 214) and others we take 

E = 320,000 Ib. per in’. 

(o+)max = 140 lb. per in’. 

v=1/3 

p = 0.0361 lb. per in*® 
For ice 1.5 in. thick we find the formula 

D = Eh? /12(1 — v?) 


that 
D 1.01 * 10° lb. in. 


Hence 1 = (D/pg)'* = 40.8 in. 
For a circular loading area of diameter lft. we find 
b=a/l = 0.147. 
From M page 180 
kei’ b = (0.6159 — log, b) (6/2) + rb? /644+.. 
This gives for 6 = 0.147 that 
kei’ b = 0.186. 
From Formula (5.2) we thus find 
P = 195 lb. 
This result seems to be in agreement with the statement made by Barnes that 
ice 1.5 in. thick will hold a man. 


From M, page 180, we have 


1 — 
ker’ b= — ni + inv — 4 (1.3659 — log.b)vi+..., 


where v = 0/2. Thus 

1+ bd ker’ b = 3xv* — (1.3659 — log.b)vt +.... 
For b = 0.147 we find 

1 +b ker’ b = $9(0.0735)? 
is a satisfactory approximation. Thus wax as given by (4.6) is 
Wmax = 7q(0.0735)?/2 pg = P(0.0735)* /2a*pg 
Hence Wmax = 0.4 in. 
Conclusion 


In the introduction of our paper we gave a brief outline of the motivating 
problems for the present paper. Throughout the paper we have tried to 
point out the simplifying assumptions that we have made in order to obtain a 
mathematical analysis. The present work should be considered at best a 
start toward the solution of the motivating problems and that a satisfactory 
solution will be obtained only by taking into account some of the factors 
which we have neglected. 
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It is hoped that formulas (4.6) and (5.2) will soon be tested experimentally 
and that their domain of validity will be established. When this is done these 
formulas may prove of some value in determining some of the elastic properties 


of ice. 

In our present work we have dealt at some length with the cases of a single 
concentrated load and also a single loading over a circular area. From these 
solutions it is easy to determine the deflection due to any number of concen- 
trated forces or any number of loadings of circular areas.’ These solutions 
will of course be required in determining the thickness of ice for vehicles of 
different types. 


Finally we would like to point out that our solutions may have applications 
in a different field. Our solutions will of course be valid for any isotropic 
plate resting on an elastic foundation. The constant pg is simply replaced by 
the modulus of the foundation. As such the solutions which we have 
determined for the deflection of the plate may prove of value in the design 
of concrete roads or airport runways. 
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ANGULAR DISTRIBUTION OF NEUTRONS AT THE INTERFACE 
OF TWO ADJOINING MEDIA' 


By B. Davison? 


Abstract 


Two media separated by a plane interface and each filling an infinite half- 
space are considered. Neutrons come in infinite number from infinity in one 
medium and the neutron density vanishes at infinity in the other medium. 
Isotropic scattering in the laboratory system of co-ordinates in both media is 
assumed and all neutrons are assumed to have the same speed. The angular 
distribution of neutrons emerging from either medium into the other is obtained 
in terms of the angular distribution of neutrons for the Milne problem, that is, 
the angular distribution of neutrons emerging from a scattering medium into a 
vacuum. The latter angular distribution is tabulated in the preceding paper 
by LeCaine. 


1. The preceding paper (1) (hereafter referred to as A) contains the tabu- 
lation of the angular distribution of neutrons emerging into a vacuum from a 
half-space filled with a scattering material. One may also be interested in 
the angular distribution of neutrons emerging from one into another scatter- 
ing medium. The purpose of the present paper is to show that this can also 
easily be obtained from the tabulations of A. 


With regard to each of the two media we make the same assumptions as 
are made in A with regard to the single medium considered there: we assume 
that one-velocity theory is applicable to both media (i.e., all neutrons have 
the same velocity and are scattered without loss of energy), that scattering is 
isotropic in the laboratory system of co-ordinates in both media and that we 
are dealing with plane geometry, each medium occupying an entire half- 
space, with the neutron density in it depending only on the distance from the 
interface. 


The case of the neutron density’s tending in one medium to zero at large 
distance from the interface would seem to be of greater interest than the more 
general case when the neutron density does not tend to zero at large distances 
in either medium. The former case is also simpler in the sense that the solu- 
tion is unique, apart from a multiplicative constant, while in the latter case 
the solution can be looked upon as a sum of two linearly independent solutions 
with arbitrary coefficients. In view of this we limit ourselves in the present 


note to the former case. 


Since we assume plane geometry, it is immaterial whether the total mean 
free paths in the two media are the same or not. For in the one-velocity 
theory in plane geometry we may use different units of length in the two 
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media without changing the form of the equations governing the distribu- 
tion of neutrons. Thus we take the total mean free path to be unity in each 
medium. The media then differ only in their absorptive or multiplicative 
properties. 

Let c be the probable number of secondary neutrons per collision in any 
medium and a = 1 — ¢; so that a> 0 for media in which capture predomin- 
ates over multiplication and a < 0 for media in which multiplication pre- 
dominates over capture. If capture and scattering are the only processes 


possible then a 
lotal mean free path 





* ™ Mean free path for capture ° 

2. Take the interface as the plane x = 0 and characterize the spaces x > 0 
and x <0 by subscripts 1 and 2 respectively. Let neutrons come in infinite 
number from infinity in the first medium and let their density vanish at 
infinity in the second medium. This implies that az> 0. We impose no 
restrictions upon the sign of a1, or upon the comparative values of a; and a2* 
We require the angular distribution of the particles entering the second 
medium. 


Let f(x) be the density in our problem, n(x), be the density which would 
exist in the first medium if the second medium were replaced by a vacuum, and 
t(x) be the density which would exist if the second medium occupied the half- 
space x > 0 with an infinite number of particles comimg from infinity and 
with vacuum in the half-space x <0. Let yy(x, uw), wa(x, uw) and W(x, uw) be 
the angular distributions corresponding to f(x), n(x), and ¢(x) respectively, 
u being the cosine of the angle between the velocity of a neutron and the 
normal to the interface drawn into the first medium. Also let v; and ve be the 
roots of the equations 

1+ ] 1+ 


= , log : log ——— - (1) 


1 — a, 2r1 l1-—- vy 1 — az 2v2 1 — vw 





Since az is positive, the roots of the second of these are real. We shall 
understand by v2 the positive root. 


The densities f(x), n(x), and ¢(x) satisfy the integral equations 


1 — a a 


2 1—a. |? 
= -—— | f(y) E(\x — yi)dy + — | f(y) E(x — y\)dy 
= 0 > —~o 


- - 


ia | n(y) E(\x — yl)dy 


0 


; | t(y) E(\x — yl)dy 
0 


co 
in which E(x) = | e** dt/t. 
1 
*Except that the description ‘‘neutrons come in infinite number from infinity in the first 
medium" is applicable only for a, > 0. For a; < 0 it should be replaced by ‘‘neutron density 
be a bounded oscillating function of position in the first medium.” 
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Taking the Fourier transforms, and denoting, following Titchmarsh (2), 
these Fourier transforms by the corresponding capital letters we have 


F*(u) li so tase tan™ a4) 4 rw) Ug as tan” *|- a 
u u 


tan”! u 


N*¥(u) E — (1 — a) | + N-(u) =0 , (3) 


tan 


8] + 7-(u) =0 | 


é 


T* (u) E — (1 — a2) 
u 
in which, according to the definitions, F*(u), N+(u), and T*(u) are regular in 
the half-planes, I(u) > R(v1), I(u) > R(v1), and I(u) > v2 respectively, while 
F-(u), N~(u), and T~(u) are regular in the half-plane J(u) <0. Also, all 
these functions are quadratically summable over the lines parallel to the real 
axis situated within the half-planes in which they are regular. Notice next 
that we can put 
Pod a) hee ‘ 
u u2+ 1 py, (u) 
(¢ = 1, 2) 
in which p* (uw) is regular in the half-plane J(u) > — 1, pj; (uw) is regular in 
the half-plane J(u) < 1, and both are zero-free in the whole plane. Equation 
(3) can now be rewritten as 


F*(u) (u? + v1?) Pr aay = — F-(u) (u? + v:*) 5 a) 
u) Pr, 
N*(u) (w? + vi) ph (u) = — No(u) (w+ 1) oe (u) (5) 


T+(u) (u? + vo®) pi (u) = — T~(u) (uv? +1) #7, (u). 


One easily sees that the left sides of (5) are regular in the whole plane. Also, 
for u tending to infinity, they behave as u? times a quadratically summable 
function, i.e., they can only be linear functions of u. Next, since F~(u), N~(u), 
and 7~(u) are regular in the lower half-plane, the linear functions in question 
vanish at u = —ivo, u = — i, and u = —i respectively. Thus (5) gives 


Ft(u) (u? + 1°) Ph na = Aj(u + iv) 
Uu 


). . 
Nt(u) (u? + v:*) pk (wu) = An(u + 4) (6) 
T*(u) (u? + v2) px, (u) = A(u + 4) 


in which Ay, An, and A; are some constants. Eliminating pt (u) and pe (u) 
_ between the equations (6) gives 


C = Nt(u) 


Fr(u) = 
() u — tivo T*(u) 


(7) 


in which € is some constant. 
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The formula (7) can easily be reinterpreted in terms of the angular distri- 
bution at the interface. For we have, for » > 0, 


i= ai - , —x/ 
— S(x)e-" "dx 
Qu 0 


-a-ay4/5 tw?) (8) 


and similar expressions for ¥,(0, — ») and y,(0, — uw). Using these expressions 
we rewrite (7) as 


¥;(0, — pw) = ie | f(\u|R)e-*dR = 
0 


Cc y,(0, = ML) 
Og he ree 9) 
A ae TS eee 


the constant C depending, of course, upon the normalizations used. 


Suppose, for the sake of definiteness, that f(x), (x), and ¢(x) are normalized 
so as to be associated each with unit current at the interface, i.e., 


1 1 1 , 
| pw (0, udp -| _ben(0, udp = | el, udu = —1, (10) 
1 


and determine, for this normalization, the values of T*(u), Nt(u), and Ft(u) 
in the vicinity of « = 0. Consider first T*+(0). T'+(u) is defined, for I(u) > ve, 


by 


1 PR ae 
T* = oe oe dx 11) 
(u) Vin I, e (x)dx ( 


and for I(u) < v2 it is the analytical continuation of (11). Now it is known 
that for large values of x the asymptotic behaviour of ¢(x) is given by 


t(x) = Be’* + O(€~*) = tas(x) + teuppi(*) (12) 


say, in which B is some numerical coefficient. Thus (11) becomes 


Tt (u) = = es J ae eaaeae a | ee \ . 


fin (— 5 — in 0 


which gives at once the analytical continuation of the definition (11) into the 
region J(u) > — vs, and in particular 
Pode oe _ 1 
T*(0) = ——4 --— t x)dx 7. le 
( ) V/2 Vo ais I, suppl (xe | ( 3) 


Tv 


Consider now the current associated with ¢(x). Following (12) we can also 
split the current into two parts, one associated with ¢,,(x) and the other with 
tsupp\(x). From the transport equation we have at once 


Wilx) — Wilxo) = a: | WYo(x’)dx’, (14) 


in which y; denotes the first moment, i.e., the current, and y is the density. 
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For fsuppi(x) the current vanishes at x = ©, and consequently the current 
corresponding to fsuppi(x) at x = 0 is given by 


co 
a2 | tsuppi(*) dx. 
0 


For ¢,,(x) the current would vanish at x = — ©, and consequently the current 
corresponding to ¢,,(x) at x = 0 will be given by 
—O ro 
a | tas(x)dx = a2 B e dx = 

0 0 
Combining these expressions and remembering that according to the normal- 
ization conditions (10) the total current is equal to — 1, we have 

a.B foo 


= a as | tsuppl (x)dx = — 1. (15) 
0 


Vo z 


Comparing this with (13) we obtain 


T+(0) = — z (16) 
ao/2r 


If a; is positive the same argument can be applied to N*+(u) and F*(u) giving 


Ft(0) = Nt(0) = sites : (17) 


ay V/2r 
We can now rewrite (6) more explicitly as 


pi (u) —1 _, pt (0) ( 
Fae ve 1— 


Fr(u) (u? + v,°) - 


Pi, (uw) 1/20 . pi, (0) 


N*(u) (u? + vi) pf (u) = eas vr pr, (0) (1 — au) 
aiV2r 


1-/ ° 9 \ 7 I 9 ° 
T*(u) (u? + v2?) py, (u) = ——= v2? px, (0) (1 — mm). 
ao 2r 

The formula (7) thus becomes 
hf) 2 N*(u) 
P(g): ee (19) 

vo t+iu Tr(u) 
and (9) 


¥,(0, Wel 
1 — pre ¥.(0, — pw) 


As in A let g(a, u) denote the angular distribution of neutrons emerging 
into a vacuum from a medium whose capture is characterized by a, normal- 
ized to unit density at the surface. uy is defined as in A, i.e., having a sign 
opposite to ours. Then using (6) of A we easily see that our ¥,(0, — u) and 
¥:(0, — wu) are connected with g(a, u) by the relationships 


(20) 
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¥(0, —p) = ee y(ai, M) 


(0, — wu) —. ¢(as, w) 
V a2 


and that (20) becomes 
1 — a3 Vi ¢y(a1, h) (21) 


v;(0, _ bh) = a , 
2Vaiaz 1 — wre (ae, p) 


in which g(a, u) is the function tabulated in Table VI of A. 

Turn now to the case of a: < 0. Since the coefficient appearing in (20) is 
independent of a; for a; > 0, one might expect it to be independent of a; 
altogether. We verify this surmise first for a, = 0. In this case, in the half- 
space x > 0, the solutions f(x) and n(x) will be associated with constant cur- 
rent and, remembering the normalization condition (10) and the fact that in 
absence of capture the current is equal to one third of the asymptotic value of 
the density gradient taken with the opposite sign, we have, for} x — ©, 


f(x) = 3x + a bounded function (22) 
n(x) = 3x + a bounded function. 


Hence, in the vicinity of u = 0, 


Fr(u) = —- : — +0 (1) 
u?./2n u 


3 i\ [7 
Nt =— — O{— 
(u) u2»/2r * (=) 


so that the first two of (12) are replaced by 


- + 
u?F+ (4) Po" eons ee © (1 _ i*)) 


pz, (1) V2n px (0) va (24) 


Ae 3 
2Nt(u) pt (wu) = ——= pf (0) (1 — iu), 
u u) pg (u) ag Pe (0) ( ) 


while (19), (20), and (21) obviously remain unaltered. Consider finally the 
case of a; < 0 (i.e., the medium filling the half-space x > 0 is a multiplying 
medium.) Now f(x) and n(x) for x > 0 are bounded oscillating functions, 
and the integrals 


| f(x)dx and | n(x)dx 
0 


0 
exist in the Cesaro-summability sense. However it is known that by applying 
the Cesaro summability we obtain the correct values of the Fourier transforms, 
so that we have 


RT aie or ie -+_|" \ de 
F+(0) = ae i, f(x) dx; N*(0) el n(x)dx, (25) 
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the integrals being understood in the Cesaro-summability sense. On the other 
hand the first moment at large distances from the interface will behave, both 
for f(x) and for n(x), as A sin {|»:|(x — xo)}, (A and xo being, of course, differ- 
ent for n(x) and f(x)) so that on the average there will be no current at large 
distances, either for f(x) or for n(x), (though at any particular cross section 
there will be some current). Then the current at x = 0 associated with f(x) 
will be, on account of (14), 








0 


a; { average with respect to xo ot | ; f(X) dX}. 





But this is exactly the definition of the value of a divergent integral in the 
sense of Cesaro summability. Thus we obtain from the normalization con- 
ditions (10), exactly as in the case of a; > 0, 









F+(0) = N+(0) = — r 
ai\/2nr 













which is identical with (17), so that (18), (19), (20), and (21) are applicable 
to the case a; < 0 without any alterations. 










So far we have discussed the angular distribution of neutrons entering medi- 
um 2. The angular distribution of neutrons entering medium 1 from medium 
2 can be also obtained without difficulty. For this purpose notice first that, 
pz, (u) 
yi (U) 





according to (4), is an even function of u. Then, using the definitions 





vi 





of p,, (uw) and p,, (u), one easily sees that 
1 

+ . 

Oia) ste enor’, 26 

‘ Py ( u) i 





Combining this with (5) and (6) one derives 


F*(u)F-(— u) (uv? + ») = Af. (27) 






On the other hand, similarly to (16) and (17), one obtains 







and (27) becomes, more explicitly, 


Ft(u)F-(— u) (u? + v;’) 










Further, similarly to (8), we have 
St nk w(t > &,) Vif (=) 
and combining this with (8) and (28) gives 









_ -_ 2 
¥;(0, ry) v,(0, _ h) = j (1 — ai) (1 — a2) ee 





aoai 1 — yey ? 
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Thus from (20) 


(0, 4) = 2 ‘a! (29) 


or, from (21) 
1—ai vi(l — pre) (aa, pL) (30) 
2Vaiaz l—pry? (ai, pz) 





¥7(0, uw) = 


Thus the angular distributions of neutrons entering medium 2 from 1, 
v;(0, — w), and of neutrons entering medium 1 from 2, ¥;(0, «) can both be found 
from the angular distributions of neutrons emerging from mediums 1 and 2 


into a vacuum. 


APPENDIX A 


Elementary Derivation of the Formula Connecting Neutron Density 
and Neutron Current at the Interface 


In passing from (20) to (21) we have used the relationship, quoted in A, 
f(0) = —= j(0) (A.1) 
J Ja ’ . 


which gives the relation between the neutron density and the neutron cur- 
rent at the free surface of a half-space. However the formula (A.1) is not 
derived there, but reference is made to a paper which has not, so far, been 
published. In view of this it appears desirable to append the derivation of 
(A.1) to the present report. 


Actually we shall prove a slightly more general result, namely that at the 
interface between two isotropically scattering media, each occupying a half- 
space, the neutron density and the neutron current are connected by the 
relationship 

. Vi . 
f(0) = j(0), (A.2) 
Vaiaz 
provided the neutron density vanishes at infinity in the medium (2). Here 
f(x) is the neutron density at the optical distance x from the interface, (0) 
is the current at the interface, the positive direction of the current being 
taken from the medium (1) into the medium (2), a: and a2 characterize the 
absorption in the two media respectively, and »; is the solution of the equation 
1 1 v 1 
7 (A.3) 


na = 
2v1 ak a Ps 1 — a, 


(a2 will be assumed positive, but a; can have either sign or be zero). 


A free surface can be treated as the surface of a medium absorbing all 
neutrons, in which case a2 = 1, and the formula (A.2) goes over into (A.1). 
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The function f(x) is known to satisfy the integral equation 


eo 
0 


fle) = a | F()E(\x — y\)dy + t f(y) E(|x — y\)dy 


(the medium (1) occupies the half-space x > 0, and the medium (2) the half- 
space x <0). Putting, for brevity, 


1 “= Gi 43 7} ly 1 — @e ° ‘| 
Ag(x) = ae I P(Y)E(\x — yl)dy + 5 c. o(y)E(|x — y|)dy (A.A) 


we can rewrite the above equation as 
(1 — A) f(x) = 0. (A.5) 
The expression for j(0) is given by 


0 
j(0) = a | f(x)dx (A.6) 


—c 

since the current through the interface should be equal to the number of 
neutrons captured in the medium (2) in unit time. On the grounds of (A.5), 
the expression for 7(0) can also be written as 

lag 

j(0) = ag Af(x)dx. 

J —0o 
On changing the order of integration and slightly rearranging we get a further 
expression for j(0), namely 


aos = i 0 
7j(0) = i a f(x) E2(x)dx — . 3 | f(x) E2(|x])dx. (A.7) 
0 


Gp 
(In the above 
00 ett Oo .—<t 
E(x) -| — dt; E2(x) -| dt. ) 
tJ 1 
In order to prove (A.2) it is convenient to consider an auxiliary function 
Q(x) defined by 
p x 24 0 
Q(x) = f'(x) —vyr \; f(y)dy —_— vj (9) re (A.8) 


ai 


Let us construct the integral equation satisfied by Q(x). If g(x) is any differ- 
entiable function such that Ag(x) exists (i.e., the integrals involved in (A.4) 
converge) then differentiating (A.4) and eliminating the derivatives of the 
E functions by integration by parts gives 


d 


dx 


a@2—~ ai 


Ag(x) = g(0)E(\x|) + Ag’(x). (A.9) 


° z 
Taking here g(x) = | f(y)dy and remembering (A.5) we see at once that 
0 


d - 
# | — J ) = 
ax | ( A) E fta)dy 0, 
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x 


so that (1 — A) | fovay is a constant. To determine this constant take 
0 


xs = 0. Then 
(l—A)] f(y)dy oe E(y)dy | f(x)dx 
0 z=0 \ 2 0 


-1= 8" kdbdey [° soeyax \. 


By changing the order of integration and using (A.7) we see at once that 
this is equal to — 7(0). Thus 


(1 — A) | sonay = — j(0). 


On the other hand, differentiating (A.5), and using (A.9) again, but with 
g(x) = f(x), gives 


(1 — A) f’(x) = SS MOE (x1) 


while A(1) can be easily worked out directly from (A.4). 
Combining these results gives 
(1 — A) Q(x) = A(x), say, 
in which 


2S (0) E(x) — v2 j(0) 2— E(x) for x > 0 
2 2a, 


H(x) = (A.11) 


<2 —_S* f(0)E(\x|) — 270) 2 —"* {2 — Ey(\x|)} for x < 0. 
2 2a, 


Let us examine now the behavior of Q(x) for large positive x. It is known 
that, considering say the half-space x > 0, we can present f(x) in the form 


f(x) = fas(x) + f(x), x>0 (A.12) 
in which f,,(x) is a solution of 
f''as(%) — v1? fas(x) = 0 
while f(x) is a small quantity of the order of 0(e~*). Hence 
f(x) — we f(x) = f(x) — ve fx) = O(e*). 


Integration of the last relationship and comparison with the definition of 
Q(x) shows that 


Q(x) = const + 0(e77) = Q(~) + O(€7), (> 0). (A.13) 
On the other hand, once Q(x) tends toa definite limit for x + «, (A.10) gives 
aiQ(o) = H(@),. 
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But according to (A.11), for positive x, H(x) = O0(e77). Hence H(o) = 0 
and consequently Q() = 0. Thus (A.13) reduces to 


Q(x) = O(e-*). (A.14) 


Notice now that H(x) is the free term of an inhomogeneous equation pos- 
sessing a solution, the corresponding homogeneous equation also having a 
solution. If we were dealing with bounded functions and the range of integra- 
tion were finite, this would immediately imply that H(x) is orthogonal to the 
solution of the transposed homogeneous equation. In our case this would mean 


1 — a 


ae 


@ “ae 0 
| F(x )f (x \dx + =" | A(x \f(x)dx == (A. ] 5) 
0 —oo 


since the solution of the transposed homogeneous equation in our case is 


(1 — ax) f (x) forx>0 
(1 — ae) f (x) forx <0. 


Actually in our case the range of integration is infinite. For an infinite 
range of integration the orthogonality of a possible free term to the solution 
of the transposed homogeneous equation does not always hold. But (A.15) 
can easily be proved directly, using (A.14) and the fact that f(x) is assumed 
to tend to zero for x > — @. 


For, similarly to (A.12), for negative x we have quite generally 
f(x) = Ae?* + Be? * + O(e!#!) (x < 0), 


in which v2 is the positive root of the equation similar to (A.3) (v2 is neces- 
sarily real since az is assumed to be positive.) If for x + — © f(x) tends to 
zero, the coefficient accompanying e”'”! should vanish, and the above formula 
gives 


f(x) = O(e~”"*!), (x < 0) (A.16) 


Consider now the left side of (A.15). Substitute for H(x) its expression 
in terms of Q(x), i.e., rewrite the left side of (A.15) in the form 


a om te 


| f(x) (1 — A) Q(x)dx + = : | 


0 


f(x) —A)Q(x)dx. (A.17) 
=i 
In the double integrals which appear if (A.17) is written out explicitly, one 
can interchange the order of integration, as one can easily verify directly 
from (A.14) and (A.16). Interchanging also the role of the two variables of 
integration we see that (A.17) is identical with 
1 =~ Gs 0 


| Q(x)(1 — A) f(x)dx + 3 | Q(x) (1 — A) f(x)dx, 
0 


ae 


l1—a, 


2 
which vanishes in view of (A.5). This proves (A.15). 


On the other hand, substitution into (A.15) for H(x) from (A.11) gives 
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= —*! £0) | = al f(x) E(x)dx + [== ’ f(x) E(\x|)dx 
Z 7 0 2 J —oo 


_ 


Er — a, {” Se 1 — ae |? 
— v°j(0) —— a f(x) E2(x)dx — = f(x) E2(\x|)dx 
0 —oc 


— »°j(0) — 2 ‘|! foaar. (A.18) 


Now the content of the first square bracket in (A.18) is, from (A.5), equal 
to f(0), of the second square bracket is, from (A.7), equal to j(0); finally the 
integral appearing in the last term is also expressible in terms of j(0) by 


means of (A.6). 


Thus we get 


Qo — a1. a a2 — a (a2 — ai) (1 — az) 
— f?(0) — v;°77(0) | + | = 0 
2 F ‘4 2a, 2aia2 





from which (A.2) follows at once. 
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PARTICLE IDENTIFICATION IN PHOTOGRAPHIC EMULSIONS 
BY THE DELTA-RAY METHOD! 


By L. Voyvopic 


Abstract 


A comparison is made in this study between the 5-ray densities expected theo- 
retically and those observed experimentally on charged particle tracks in Ilford 
G5 emulsions. The energies of 6-rays are determined here by means of grain- 
counting measurements and energy-loss considerations. Using this method of 
energy determination and the theoretical formula for 6-ray density, good agree- 
ment has been found with the experimental distributions on proton, u-meson, 
deuteron, and a-particle tracks. The agreement is satisfactory both in absolute 
value and in the shape of the 6-ray density distribution. These results are con- 
sidered to show that the 6-ray technique is suitable as an absolute method of 
charged particle identification in photographic emulsions. 


Introduction 


The study of 6-rays, or ‘‘knock-on”’ electrons, in photographic emulsions has 
proved to be a valuable technique in identifying energetic heavy nuclei in the 
primary cosmic radiation (2, 4), and the energetic fragments (8) occasionally 
emitted in nuclear disintegrations or ‘‘stars’’. The tracks of these heavy par- 
ticles appear under the microscope as solid filaments of continuous developed 
grains, accompanied by large numbers of short spurs due to ‘‘knock-on” 


electrons. 

Heretofore, the theoretical distribution ef 6-rays along the tracks of charged 
particle has been assumed as valid for photographic emulsions. For com- 
parison with the 6-rays observed on experimental tracks, theoretical curves 
have been normalized by calibration with a-particles of known energy (2, 4), 
or with recognizable protons and mesons (8). In the present study a more 
detailed comparison is made between the calculated and observed 6-rays on 
recognizable tracks in Ilford G5 electron sensitive emulsions. In particular, 
the absolute values of the 6-ray densities are used in this comparison in order 
to test the validity of the 6-ray technique as an absolute method of charged 
particle identification. 


Theory 


The theoretical production of 6-rays may be derived from the differential 
collision probability x(£, E’)dE’dx, which is the probability for a particle of 
charge Z, energy E, and velocity 8c, traversing a thickness dx, to transfer an 
amount of energy between E’ and E’+ dE’ to a free electron. In the classical 
treatment x(E, E’) is given by the Rutherford formula (7) as 

2u.pZ*Cdx dE’ 


x (E, E’) dE'dx = ————_—_ 


aif (1) 
e (EP 


1 Manuscript received February 2, 1950. 
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(2’) 


where yu, is the electron mass, p is the density, and C = 0.150 * +— for the 
material traversed. A) 


The 6-ray density , or the number of 6-rays per unit length, involving 
energy transfers 2 E’»jn, is obtained by integrating x(E, E’)dE’ over the 
2u.8? 

PHB Mev.). 


energy region EF’ »:n to E’ max (where E’ maz = - } Using the values 


/ 


ae te J ; ; 
6 = 4.0 gm. per cm.*, — = 0.45 for Ilford emulsions (3), one obtains the 





formula 3 27 X ee crs 1 — 8 
B° Seitie 1.0287 


The quantum mechanical treatment of x(£, E’) given by Bhabha (7) leads 
to a small correction term in the above formula for n. In the case of E’» in = 10 
kev. this correction is not more than 3%, and it is ignored in the present study. 


| 6-rays/100 microns. (2) 


From Formula (2) above it is seen that the 6-ray density is dependent on the 
velocity Bc of the primary particle and on the minimum energy E’ mi, chosen 
for the 6-rays. For comparison with experiment, the velocity may readily be 
related to the observed residual range of the particle by means of the experi- 
mental proton range-energy data (1, 5) for photographic emulsions. However, 
no satisfactory range—-energy results for slow electrons are yet available for 
determining F’,,;,. Recent measurements by Ross and Zajac (6) with a 6-ray 
spectrograph cover the electron energy region 30 to 100 kev., but, as will be 
seen in the following sections, the 6-ray energies concerned in this study lie 
below 30 kev. To overcome this difficulty, a method has been developed for 
relating the number of grains in a 6-ray track with the electron energy. This 
method makes use of grain-density measurements on protons and mesons, 
together with theoretical energy-loss calculations for charged particles in 
photographic emulsions. 

Experimental Method 


Illustrations of 6-rays are seen in the photomicrographs of Fig. 1. Fig. la 
shows a section of a proton track and its 6-rays at about 1000 microns from 
the end of the proton range. Some of the 6-rays are difficult to identify in the 
photograph, their ranges being only a few microns and Coulomb scattering 
very pronounced. Several 6-rays lie in the same direction as the proton track. 
Under the microscope one can follow each 6-ray by focusing on the individual 
grains, and in fact grain-counting is much more convenient and also more 
reliable than range measurements on these very short electron tracks. The 
long 6-rays shown in Figs. 1b and lc were observed along the track of a long- 
range proton, at 5.4 and 3.4 mm. from the end of its range. Unfortunately, 
clear 5-rays of this type are very rare and so have little statistical significance 
in the study of the 6-ray density distribution. Fig. 1d shows a section of an 
a-particle track with its more numerous 6-rays. Here again range measurements 
on the recoil electrons are very difficult whereas grain-counting is more feasible. 
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From experience it has been found that an optimum criterion exists for 
counting 6-rays. This empirical criterion in the present Ilford G5 emulsions 
requires that the electron track contains four or more developed grains, as- 
suming that one of the 6-ray grains lies inside the track of the primary particle. 
Under these conditions enough 6-rays are counted for satisfactory statistics, 
and at the same time the 6-rays are counted with good reproducibility after 
one has gained a reasonable amount of experience. Similar procedures for the 
selection of 6-rays have been used by other authors (2, 8). 


25k 


l= ig. ld 


Fic. 1. Photomicrographs of delta-rays on charged particle tracks. 


Energy Determination of 5-rays 


In order to relate a chosen minimum number of grains in a 6-ray with its 
corresponding energy E’»in, the energy-loss per developed grain has been 
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investigated. Fig. 2 shows the results of grain-density observations on the 
charged particles whose 6-rays were also studied. For regions of track con- 
taining more than 160 grains per 100 microns, the method of counting gaps 
rather than grains was used in these observations. Fig. 3 shows the grain- 
density results plotted versus E/M rather than residual range, and also the 
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Fic. 2. Grain-density observations on long range tracks. 
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Fic. 3. Grain-density and specific energy-loss for singly charged particles plotted versus E/M. 
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specific energy-loss k, calculated for the emulsion. It is seen that the experi- 
mental grain-densities for the various types of particles all follow the same 
variation with E/M, which suggests that the grain-density of an electron track 
will also follow this same variation with E/M. (The 12 mm. meson in Fig. 2 
which left the emulsion just before the end of its range is included in Fig. 3 
as a u-meson. ) 

The specific energy-loss k, shown in Fig. 3 is for energy transfers of less than 
10 kev. obtained from the Bethe-Bloch formulas (7) and the composition of 
Ilford emulsions (3). Small energy transfers (k,) have been considered rather 
) the latter includes 
dx 


than the total energy loss by ionization and collision ( 


energy transfers to 6-rays whereas in the grain-counting results of Figs. 2 and 
3 the 6-rays were ignored. 


From the ratios of energy-loss and grain-density in Fig. 3 one may readily 
obtain the energy-loss per developed grain as a function of E/M and hence 
as a function of electron energy. By integrating over the electron energy one 
then obtains the total number of grains in the corresponding electron track. 
The results of this procedure for electrons with energies up to 30 kev. are shown 
in Table I. 

TABLE I 


PREDICTED RELATIONSHIP BETWEEN ELECTRON ENERGY AND 
TOTAL NUMBER OF GRAINS IN ELECTRON TRACK 


: k,/N, 

Electron energy, rN 
ae ve kev. per 
kev. 


grain 


Number of grains 
in electron track 


0.01 
0.015 
0.02 
0.03 
0.04 
0.05 
0.06 


meee DoD oo 


By interpolating the results of the last column in Table ! we find that a 


6-ray containing four developed grains is expected to require about 11 kev., 
and one containing 10 developed grains to require about 21 kev. 

Beyond E/M = 0.06, the energy per developed grain remains constant at 
about 1.2 kev. per grain. The minimum calculated energy-loss is 38 kev. per 
100 microns, which leads to an expected minimum in grain-density of 32 grains 
per 100 microns. It is interesting to note that the mean grain-density for 50 
random “‘thin’’ tracks (assumed to be due to relativistic u-mesons) in this 
batch of emulsions was found to be 33 grains per 100 microns, and for 15 elec- 
trons from u-meson decay was found to be 37 grains per 100 microns. These 
latter would be expected to have rather more than minimum ionization. 
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Observations and Results 


Fig. 4 shows the 6-ray densities for protons. The three smooth curves are 
those calculated for values of E’ in equal to 8, 10, and 12 kev. The theoretical 
densities have been converted from velocity to residual range dependance by 
means of the range-energy results of Lattes et. al. (5) and of Bradner et.al. (1). 
The block curve is that obtained experimentally as the mean of four protons 
of residual range not less than 5 mm., using the criterion that each 6-ray con- 
tains four or more developed grains. The experimental curve is seen to agree, 
both in absolute value and position of the peak, with an electron energy E’ min 
hetween 10 and 12 kev. Taking into account the normal statistical uncertainty 
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Fic. 4. Delta-rays on protons. 


due to the finite numbers of 6-rays, and a comparable uncertainty found in 
reproducibility of observations, the experimental results are found to be con- 
sistent with £’,,;,= 11 + 1 kev. 


The results for six u-mesons (accompanied by decay electrons) of residual 
range not less than 2 mm. are shown as the block curve in Fig. 5. Here the 
experimental 6-ray densities indicate a value of E’,:n between 11 and 7 kev. 
corresponding to four or more developed grains. Considering the uncertainties 
due to statistics and reproducibility of observations, the u-meson 6-rays are 
found to be consistent with EF’ »in= 9 + 1 kev. 


Delta-ray density measurements have been made on two long tracks of 
residual range 12 mm. and 1400 microns which were ejected from an inter- 
esting cosmic-ray star containing 32 prongs. These two tracks were both 
ejected in a backward direction with respect to the axis of a wide-angle shower 
containing eight charged relativistic particles. 
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Fig. 6 shows the observed 6-ray densities for the 12 mm. prong. This track 
was tentatively identified in Figs. 2 and 3 as that of a deuteron with initial 
energy about 85 Mev., although its identification as a 100 Mev. triton cannot 
be excluded on the basis of grain-counting. The three theoretical curves in 
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Fic. 5. Delta-rays on p-mesons. 
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Fic. 6. Delta-rays on a high energy deuteron. 


Fig.” 6 are for a triton, deuteron, and proton based on E’m;n= 10 kev., the 
mean energy found above for proton and yu-meson 6-rays. It is seen that the 
experimental curve for 6-rays having four or more developed grains on the 12 
mm. track is in better agreement with the theoretical curve for a deuteron 
than with those for a triton or a proton. The maximum of the observed 6-ray 
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density curve in Fig. 6 occurs at a residual range beyond 1000 microns whereas 
the maximum for protons (Fig. 4) occurs near 500 microns. 


The bottom curves in Fig. 6 show the observed densities for 6-rays having 
10 or more developed grains on the 12 mm. track (block curve), and the theo- 
retical curve for deuterons with EF’,,;,= 20 kev. (smooth curve). These two 
curves seem to be in quite good agreement. The value E’,,;,= 20 kev. was 
chosen in computing the theoretical curve because the results of 6-ray energy 
determination in the preceding section indicated that an electron would require 
about 21 kev. to produce 10 developed grains in a 6-ray track. 
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RESIDUAL RANGE IN MICRONS 
Fic. 7. Delta-rays on alpha particles. 


The observed 6-ray densities for the 1400 micron track accompanying the 
12 mm. deuteron are shown as the broken block curve in Fig. 7. The full block 
curve is that for a similar track found in an eight pronged star in another plate 
in the same batch of emulsions. These two tracks were both too dense for iden- 
tification by grain-counting, indicating an atomic number greater than one. 
Comparison of the experimental curves in Fig. 7 with that for protons in Fig. 4 
shows a ratio of nearly four in 6-ray densities for the same residual range. This 
is just the result to be expected from Formula (2) if the tracks were produced 
by a-particles. It will be seen by comparing the experimental curves for 6-rays 
having four or more developed grains with the theoretical curves for a-particles 
with E’»;, equal to 10 kev. and 12 kev., also shown in Fig. 7, that there is 
agreement with EL’ »;n= 11 + 1 kev. 


No examples were found in the present batch of emulsions of tracks showing 
6-ray densities sensibly higher than those shown in Fig. 7 for a-particles. 
These emulsions were flown in aircraft at an average altitude of 14,000 ft. 
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Discussion 


In the section on energy determination of 6-rays, an electron energy of about 
11 kev. was predicted for 5-rays having four developed grains. From the experi- 
mental results discussed in the preceding section it is seen that the numbers 
of 6-rays for all the particles studied are in agreement with theory if it is 
assumed that a 6-ray with four developed grains has an energy of 10 + 2 kev. 
Similarly, the numbers of 6-rays having 10 or more developed grains is in 
agreement with theory (taking E’,,;,= 20 kev.) whilst the predicted energy 
for an electron track containing 10 grains was about 21 kev. 


If we assume that the 6-ray energy determination method used here is valid 
then we may conclude that the theoretical formula for the production of 6-rays 
in photographic emulsions has been verified both in absolute value and in 
distribution. Thus, observations on the numbers of 6-rays and their distri- 
bution with residual range allows an absolute method for charged particle 
identification. 


It is difficult to see how the validity of the method used here for determining 
6-ray energies may be tested directly since the electron energies are so very low. 
For example, external electrons of 10 kev. energy would all be absorbed in the 
surface of the emulsion, and observations on grain counts or range under such 
conditions would be very poor. Although the grain counting method used here 
is considered to be reliable, a more direct test of the 6-ray theory may be made 
with higher energy electrons for which experimental range-energy data are 
available (6). This method would be feasible only for energetic heavy nuclei 
where the numbers of high energy 6-rays are much higher than is the case with 
the light particles available for the present study. 
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ON THE ORIGIN OF TEN CENTIMETER RADIATION 
FROM THE POLAR AURORA! 


By P. A. ForsytH?, W. PETRIE, AND B. W. CURRIE 


Abstract 


Short-lived bursts of 10 cm. radiation from auroral displays have been received 
by radar equipment. The sources of continuous radiation from a partially 
ionized medium are briefly discussed. From a knowledge of the constants of the 
equipment used it is deduced that the power density at the receiver is at least 
2 X 10-° watts per (meter)?, and it seems that the most likely source of this 
radiation is a plasma oscillation of the ionized volume associated with the auroral 
display. If this is so, the electron density in at least localized regions must be of 
the order 10" per cm.’ 


Introduction 


A radar set operating at a frequency of 3000 Mc. has detected short pulses 
of radiation from auroral displays (6). The pulses arrive in a random manner, 
occuring in bursts lasting only a small fraction of a second, and are more 
numerous with bright displays exhibiting a definite ray structure. In view of 
this result it is of interest to consider the mechanisms responsible for the 
emission of long wave length radiation from the polar and nonpolar auroras. 
The following discussion is concerned with the emission of radiation at a 
frequency of 3000 Mc. (wave length of 10 cm.). However, the calculations 
can readily be modified to give the emission at other wave lengths. The 
emission of long wave length radiation from a partially ionized medium may 
result from a number of processes which we will discuss in turn. 


1. Atomic and Molecular Transitions 


Transitions between atomic levels and also between molecular levels pro- 
duce radiation in the ultraviolet, visible, and infrared regions of the spectrum, 
but little or no radiation in the radio spectrum. The atoms and molecules 
in the earth’s upper atmosphere are excited to only the lower energy levels, 
and transitions between these levels do not produce radiation in the wave 
length region in which we are interested. 


2. Capture Processes 


(a) Electron Captures by Protons 

We have no definite information regarding the number of protons in the 
earth’s upper atmosphere, although cosmic ray studies indicate that these 
particles do exist at high altitudes. A study (21) of the intensities of the Balmer 
lines in the auroral spectrum indicates that the capture process is not the pri- 
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mary mechanism responsible for the production of radiation from hydrogen 
atoms. However, although the capture of an electron by a proton is probably 
an unimportant process in the upper atmosphere, the results of this capture 
process can be predicted quite accurately, and give us information on similar 
processes involving atomic and molecular ions which certainly exist in large 
numbers. 


A number of investigations (20) indicate that the excitation temperature 
during an auroral display is of the order of 5000° K. The number of ions 
present during an auroral display is unknown, whereas the number present 
in the normal night sky is reasonably well determined in the E and F layers. 
Hence we will consider the processes which might produce long wave length 
radiation from the nonpolar aurora, and attempt to estimate the differences 
which occur during an auroral display. It should be realized that the electron 
temperature of the free electrons in the ionosphere is probably much less than 
the excitation temperature which exists during an auroral display. 


Table I lists certain ionospheric data given by Grimminger (10), and Bates 
and Massey (2). 


* 


TABLE I 


Height Temperature r . 4 ’ 
ent, eelvin , No. particles/cm.’ | No. electrons/cm.* 


500 
1000 
2000 
2000 
2000 
2000 
2000 


n | 


~ 


104 
105 
108 
10° 
10° 
10° 
10° 


10" 
10" 
10° 
105 
10° 
107 
10? 


om 


XK XK KK KK 


x 


Tt em 


~ 


7 
6 


o 


i 
xX KKK KX | 


j 
| 
| 
| 
| 


The values of temperature and electron density above 300 km. are not given 
by the above mentioned authors but are assumed values. The mean kinetic 
energy of an electron at temperatures of 500°, 1000°, 2000° K. may be found from 
the classical relation 1/2 mv? = 3/2 kT. The mean kinetic energy values are 
10.4 X 10-", 2.1 K 107%, 4.1 K 107" erg. Electrons with these energy 
values, when captured on the level » = © of the hydrogen atom will emit 
radiation of wave lengths 191,000, 94,500, 48,400 A. However, a number of 
the electrons will have less kinetic energy than the mean values, and our 
problem is to calculate what fraction of the electrons have velocities such 
that when captured on the level n = ©, radiation of wave length 10° A is 
emitted. The velocity of a free electron, which when captured on the level 
n = © radiates at \ = 10° A, is 2.09 X 10° cm. per sec. If we assume that the 
free electrons are in thermal equilibrium, then the fraction of electrons having 
a velocity between v and (v + dz) is given by the Maxwell law, 


mv? 


3/2 
= - ae a2 ~2RT - 
f(v) = (57) v-e 
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Substituting v = 2.09 X 10° cm. per sec. and the three temperature values, 
we find that f(v) ~ 5 XK 10™, 1.9 X 107", 0.66 XK 10~™". Hence, very few 
electrons will emit radiation at \ = 10° A if the captures occur on level n = ©, 
The energy radiated by an electron of velocity 2.09 K 10° cm. per sec. when 
captured on the level n = © is 19.6 K 107% erg. Hence, the energy radiated 
per cubic centimeter per second at a wave length of 10° A due to electron cap- 
tures on the level 2 = © will be the product of the figure given above for the 
energy release per capture, and the number of captures per.cubic centimeter 
per second on the level = ~ of electrons with velocities of 2 X 10° cm. per 
sec. Actually very few captures will take place on the higher levels of the 
atom at the temperatures we are considering. The number of captures per 
cubic centimeter per second per unit velocity range may be expressed as 
N;N.f(v)ve, where N; is the number of ions per cubic centimeter, NV, the 
number of electrons per cubic centimeter, f(v) the Maxwell velocity distri- 
bution law, v the velocity of the free electron before capture, and o the target 
area of the ion for capture of the electron. If we substitute the value of ¢ as 
given by Menzel and Pekeris (18), the values of f(v) and v which are approxi- 
mately 10-" and*2 X 105 cm. per sec. respectively, and assume the product 
N;N, is 10", it is apparent that the number of captures on the higher levels 
is very small. For instance, the number of captures per cubic centimeter per 
second of electrons with velocities 2 X 10° cm. per sec. on the level ” equals 
1000, is around 10~. Hence the radiation emitted per cubic centimeter per 
second at 10° A due to this capture process will be of the order of 10-® erg or 


less. 


(b) Electron Captures by Atomic and Molecular Ions 

The positive ions which we might expect to find in the upper atmosphere in 
any abundance are Ot, O.t, Nt, Net. Although the results of the capture 
process for these ions cannot be predicted accurately as for the case of hydrogen, 
captures on the high levels of these ions will still be very few, and the radiation 
emitted per cubic centimeter per second at a wave length of 10°A will probably 
be comparable to the value given for the hydrogen problem. 


(c) Electron Captures by Neutral Hydrogen and Oxygen Atoms 


The properties of the H~ ion have been carefully investigated by Chand- 
rasekhar and Breen (3), and of the O~ ion by Bates and Massey (1). The 
electron affinity of H~ is around 0.7 ev. and that of O~ around 3 ev. Hence the 
capture of an electron by a neutral hydrogen or oxygen atom produces a contin- 
uous spectrum at much shorter wave lengths than the region in which weare in- 
terested. Bates and Massey (1) also discuss the properties of the O.~ ion. The 
energy radiated when a very slow electron attaches itself to an O2 molecule 
is of the order of 0.1 ev.; this energy is at shorter wave lengths than the region 
in which we are interested. 
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3. The Cascade Process After Capture 


After capture, an electron cascades to lower levels and certain transitions 
may produce radiation in the radio spectrum. We will consider the cascade 
process for hydrogen and use the Bohr expression to calculate what n-values 
must be involved in order to give radiation at \ = 10°A (¢ = 0.1 cm.~) ice., 
use 7 = "a , and find that m, = 121, mz = 120, gives} ~ 0.1 cm.7}. 
Transitions between other levels of large » number will also produce radiation 
of wave length 10°A. But we have pointed out that very few captures take 
place on a high level such as m,:= 121, and then only a fraction of these 
electrons drop to the next lowest level n.= 120. We may calculate the number 
of captures on the level m;= 121, and the number of downward transitions to 
the level m.= 120 from the work of Menzel (17). The calculations indicate 
that the radiation emitted from this cascade process is negligible compared 
with other radiation sources which will be discussed later. We cannot make 
similar calculations of the cascade process for other particles, but it seems safe 


to say that the same arguments apply to these particles. 


4. Free—Free Emission Processes 


(a) Electrons and positive tons 

Using the classical orbit theory we may consider a free electron approaching 
a nucleus as moving in a hyperbolic orbit, just as a bound electron moves in an 
elliptical orbit. A free electron may undergo a free—free transition in the field 
of the nucleus, and either absorb or emit energy. Free—free transitions produce 
continuous spectra. 

This problem was first investigated quantitatively by Kramers (14) in con- 
nection with the continuous X-ray spectrum. Applying the Bohr correspon- 
dence principle, he derived an expression giving the energy radiated as the 
result of a transition between a free electron and a positive ion. The emission 
in ergs per cubic centimeter per second in the frequency range v to (v + dy) is 

_ o2m Ze® dy 
[= NiNe an? 
3V3 em? v 
where ; = the number of ions per cubic centimeter, 
, = the number of electrons per cubic centimeter, 
= the effective charge of the ion, 
= the mean velocity of the free electrons. 


Eddington (5) has used this expression to calculate the emission of energy 
from interstellar space. Gaunt (8) has given a correction to the simple Kramer 
formula. Reber (24) has used the corrected expression to calculate the emission 
from interstellar space. We will first apply the familiar Kramer equation to 
compute the emission from the nonpolar aurora, and then show that a more 
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accurate expression yields comparable results. The positive ion density in the 
night sky is unknown but is not likely to be much different from the electron 
density. Assume the product VN, in the E and F layers is 10° and 6.25 & 10" 
respectively. Substitute these figures and temperature values of 500°, 1000°, 


2000° K., and v= ‘ee 
Table II. m 


1/2 
) . The value of J at various. hefghts is given in 


TABLE II 


Height, km. I, ergs cm.— sec.—! 


100 2x. 10-%° 
200 Lx 1 
300 o- 4O-8 
400 
500 
600 
700 


Menzel (17) has given more accurate expressions for the free—free emission 
based on quantum mechanics. The emission per cubic centimeter per second 
from a free electron moving in the field of a proton is ~ 

N:NeRK _pyer 

— 7) = é h RT dy, 

2RT*’’ 
R is the Rydberg constant in terms of the frequency, and K a numerical 
constant. The frequency equivalent to \ = 10 cm. is v = 3 X 10°. Substitute 
as before V;N.= 10°, and 6.25 & 10", and temperature values of 500°, 1000°, 
2000° K. The values of J are given in Table III. 

TABLE III 


Height, km. I, ergs cm.— sec.—! 


100 3 X 10-*! 
200 LS. 10-** 
300 9 x<.10-* 
400 : <x 20-7 
500 1x 10-” 
600 Lx i 
700 LX 107% 


The results are close to the values obtained from the Kramer equation. 


(6) Electrons and atoms 

The free-free emission from an electron in the field of a neutral atom has 
been studied theoretically by Nedelsky (19). The radiation emitted per cubic 
centimeter per second from one electron encountering one atom is 


2 3 
* f(n) Va? (: _ 2) (2 - *)o, 
3c3 vo Vo 


where 1/2 mv?= 3/2 kT = hyo, and f(m) is a function of the parameters of the 
atomic model. This equation applies to slow electrons and may be used when 
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Va? < 1, where V is the energy of the incident electron in electron volts and a 
the radius of the screening electron shell in angstroms. a is of the order of 1A 
for the case we are considering. For temperature values of 500°, 1000°, 2000° K.., 
Va?= 0.06, 0.13, 0.26. The mean electron velocities at the above temperature 


values are 1.5 X 107, 2.1 X 10’, 3.0 X 107 cm. per sec. Nedelsky gives values 
for f(n) and we use the value f(m)= 7. If we substitute in the equation the 
total number of particles and electrons at the various heights, the J values 


are as given in Table IV. 
TABLE IV 


Height, km. I, ergs cm.— sec.—! 
1.3 X 10-* 
1.6 X 10-* 
1.8 X 10-77 
3.2 x 10-* 
x< 10-79 
x 10-29 
xX 10-30 


A more accurate expression for the free-free emission from an electron in 
the field of a neutral hydrogen atom has been given by Chandrasekhar and 
Breen (3). These authors have computed free—free absorption cross sections, 
and their values may be converted to emission transition probabilities through 
the use of the general relationships between emission and absorption. For 
discrete transitions the relations between emission and absorption are given 
by the well known Einstein A and B values. The relations in the case of free— 


free transitions have been given by Cillié (4). Chandrasekhar and Breen have 
pointed out that the earlier calculations of the free—free emission or absorption 
are inadequate, since these calculations are based on the Born approximation 
which does not accurately describe the motion of an electron in the field of an 
atom. Using the results of Chandrasekhar and Breen, and the relations between 
absorption and emission as given by Cillié, we find that the free—free emission 
values are about 10 times the values given in Table IV. 


The total radiation due to free—free transitions received from the nonpolar 
aurora, will include radiation which originates at different heights. The element 
of volume of a sphere is r°drdéd¢@. If we consider the radiation from a volume 
the surface of which subtends unit solid angle at the receiver, then dV = r’dr. 
SS te ; , ah ; . Irdrdvy — Idrdv 
rhe intensity of the radiation received is = ——— ergs per sec. per 

4rr* dor 

cm.? per unit solid angle. From the previous discussion, it appears that the 
mechanism of free-free emission gives the greatest radiation energy at a wave 
length of 10°A. Hence if we use the results obtained from the Chandrasekhar 
equation, the value of the radiation intensity received from the night sky is 
6.6 X 107° dy ergs per sec. per cm.” per unit solid angle. This is the radiation 
received from the E and F layers. The radiation received from the atmosphere 
above the F: layer is small compared with this figure. 
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These formulas for the free—-free emission apply to a mixture of electrons and 
one type of atom or ion, and the results must be modified when we are dealing 
with a mixture of gases such as are present in the earth’s upper atmosphere. 
In order to solve the problem exactly it is necessary to know the probabilities 
of all the possible reactions between the electrons and the ions, atoms, and 
molecules of the upper atmosphere. It is likely that the figures given represent 
the maximum free—free emission from the upper atmosphere, and the radiation 
from this process may actually be much less. 


5. Thermal Emission 


It is difficult to estimate the intensity of thermal emission from the night 
sky since the upper atmosphere does not radiate as a black body. The radiation 
density in an isothermal enclosure at long wave lengths in a wave-length 
interval dd is given by the Rayleigh—Jeans law, ¥,dA = 8rkTA~‘dX. Since the 
total radiation per cubic centimeter is a7‘, the fraction of the total radiation 
_ 8akTd~*dd 

a 
centimeter per second from the surface of a black body at wave length \ = 
0.450) *TdX. The surface area of the emitting surface if it subtends unit solid 
angle at the receiver = r*, Then the radiation per second at wave length A from 
this surface received by 1 cm.’ of the receiver = ta ll = 0.204 X 

iS 4rr° 
10~°\~*Tdx. In terms of frequency this expression is 6.8 K 10~"%Tdv. Although 
it is difficult to estimate the value of T which should appear in the equation 
for thermal noise, it appears that the free-free emission and thermal noise are 


at wave length A = 0.45\-*7—*dy. + The emission per square 


comparable in intensity. 


Comparison of the Calculated Emission and Sensitivity of 
the Radar Receiver 


The present investigation with radar equipment has shown that the intensity 
of radiation of wave length 10 cm. increases during an auroral display, and 
this increase is almost certainly due to the fact that the ionization density is 
greater in the region of the aurora. The radiation is received in bursts lasting 
only a small fraction of a second, and the intensity of the radiation depends 
upon the brightness of the aurora. A similar relationship exists between the 
appearance of a solar flare and an increase in the solar noise. 


A simple calculation will show that the intensity of 10 cm. auroral radiation 
must be considerably greater than the figure given for the energy due to free— 
free transitions, if this radiation is to be recorded by the radar equipment 
presently in use. If we use the value already calculated for the free-free emis- 
sion, i.e., 6.6. X 107! dv ergs per sec. per cm.” per unit solid angle, take the band 
width, cross-sectional area, and solid angle accepted by the radar receiver as 
1.2. X 10° c.p.s., 10 cm.?, and 107? steradians respectively, then the power 


received by the equipment is, 6.6. X 107° 1.2 K 10° 10‘ 107° 1077= 
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8 X 1078 watts. The sensitivity of the radar receiver was estimated by observing 
the echo received from a distant target of known range and dimensions. From 
the amplitude of the observed echo and the known constants of the transmitter, 
the detectable power density was calculated to be 2 K 107 watts per (meter)*. 
Since the value for the target gain used in this calculation was unity, it is 
certain that the resulting power density represents a minimum value. The 
actual detectable power density is probably somewhat greater than this value. 
Hence, if the radiation is the result of free-free transitions, the ion density 
during an auroral display must be increased by a factor of 2.5 & 107 over the 
normal night sky values, at least in localized regions. This calculation neglects 
any changes in the emission due to temperature changes during an auroral 
display ; temperature changes are of secondary importance as far as this problem 
is concerned. If the above figures are correct and the radiation is the result of 
free—free transitions, the ion density in localized regions during an auroral 
display must be of the order of 10" per cubic centimeter. It has been suggested 
that the value which we used for the electron density in the F layer is too low. 
However, even if a higher density of 2 X 10° per cm.’ is used, the conclusions 
which we have reached are unchanged. 


Discussion of Ion Density During an Auroral Display 


Little is known regarding the ion density during an auroral display. The 
interpretation of the few radar studies of the aurora is not entirely clear (12, 15). 
Since there is evidence that during at least certain auroral displays protons 
or hydrogen atoms enter the earth’s atmosphere, it is worthwhile showing that 
a diffuse stream of such particles may produce a large number of ions in the 
upper atmosphere. 

The excitation and ionization effects of particles passing through air have 
been summarized in an article by Francis and Jenkins (7). Positive ion col- 
lisions have been investigated theoretically by Massey and Smith (16). The 
cross section for excitation and ionization of N2 by proton collisions has been 
measured by Ramsauer, Kollath, and Lilienthal (23). Their results indicate 
that the cross section for a 10,000 ev. proton is of the order of three times the 
gas kinetic cross section. The gas kinetic cross section of N» is around 11.3 X 
10~'® cm.” (13, chap. 14). Hence the excitation and ionization cross section 
of No is 34 X 107! cm.*. If we assume that the cross section of all atmospheric 
particles is of the same order, the number of excitations and ionizations at the 
100 km. level as a proton travels 1 cm. is 34 X 107° X 8.5 K 10"°= 0.3. 


The ionizing power of protons in air at 1 mm. pressure has been measured 
by Gerthsen (9). Some of his results are given in Table V. 


TABLE V 


Proton energy, ev. No. of ion pairs/cm. 






23 
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The above figures indicate that a 10,000 ev. proton produces around 13 ion 
pairs per cm. Now the number of atoms per cubic centimeter at 1 mm. pressure 
and 0°C. = 3.56 & 10! Hence the number of ions per centimeter produced at 
different heights is given in Table VI. 


TABLE VI 


Height, km. No. ions/cm. 
100 0.03 
200 0.5 X 10-+ 
300 0.2 x 10 


The previous calculation indicated 0.3 excitations and ionizations per centi- 
meter are produced in the E layer. If these figures are reasonably accurate it 
means that the majority of collisions produce excitation rather than ionization. 
Since the energy loss per ionization is 36 e.u. and the ionization potential of Ne 
is 15.6e.u., the energy of an electron freed from N2is around 20e.u. This electron 
may itself produce one ionization per centimeter. Hence the total number of 
ions produced by incoming protons is given in Table VII. 


TABLE VII 


These calculations are admittedly approximate but will suffice to illustrate 
the point we wish to make, that a low density stream of protons will produce 
a large number of ions in the upper atmosphere. 


We have no information on the density of the solar material which reaches 
the earth’s atmosphere. However, we may estimate what may be the density 
of a proton stream which will produce any given ion density per second in the 
earth’s atmosphere. The equation relating these quantities is mv X 0.06 = mo, 
where 7; is the density of the incident stream, v the velocity of the incident 
stream, and m2 the number of ions formed per cubic centimeter per second in 
the atmosphere. If v is the velocity of a 10,000 ev. proton, and mz is 10", then 
nm; is 1.2 X 10°. 


This calculation indicates that the solar stream need not be very dense in 
order to produce a large number of ions in the upper atmosphere. Since the 
proton stream causes excitation as well as ionization, it is likely that a part of 
the radiation from excited molecules will be in the ultraviolet region of the 
spectrum and produce some secondary ionization. 




















































4r X 1014 


Radiation from a Plasma Oscillation 


" ae . ; dw 
The rate of radiation from the classical oscillator = i 
t 
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The short duration of the bursts of radiation suggests that the electron 
density decreases rapidly after the passage of the ionizing particles, and that 
these particles arrive in discrete bundles. 


(25). The frequency of oscillation of an electron plasma is given as 


3c3 






The rate of expansion of a cloud of electrons is readily determined. A simple 
calculation indicates that if the initial electron density of the cloud is 10", the 
radius of the surface of the cloud expands by a very large factor in 1 yu sec. 
The presence of positive ions will certainly modify the expansion process, but 
since the electrons produced by ionization will have high velocities relative to 
the positive ions, the electron density may decrease rapidly owing to the 
diffusion of these particles through the ionized volume. It seems likely that 
conditions in the ionized volume are such that this volume may execute a 
plasma oscillation, the energy being derived from the stream of particles passing 


through the ionized region. 


Whether or not the short duration of the bursts of radiation is due to a 
decrease in the electron density, and to the fact that the ionizing particles 
arrive in discrete bundles, we can estimate the energy radiated from a plasma 
oscillation by using the expression for the radiation from the classical oscillator. 
The theory of plasma oscillations was first investigated by Tonks and Langmuir 


en 
—=§ 


™m™ 


where e and m are the charge and mass of the electron, and m the electron 
density. If we consider that the electrons in the upper atmosphere are forced 
into oscillation by a proton stream, the electron density must be 1.1 & 10" 
per cm.’ in order that the plasma oscillate with a frequency of 3 X 10°. 

is 162'vig?A? 


, where 


q is the charge on the oscillator and A the amplitude of the oscillation. If the 
electron density is 1.1 X 10" per cm.*, the power radiated is 3.9 K 10~*A? watts. 
The power collected by the receiver from a radiating volume of 1 cm.’ at a 


3.9 X 10-6 X 10142 _ 


height of 100 km. is 3.1X10~'"A? watts. Arough estimate 


of the amplitude of the oscillation may be obtained as follows. The electrons 
produced by ionization have velocities equivalent to an energy of 20 ev., i.e. 
2.6 * 108 cm. per sec. If the frequency is 3 X 10° the amplitude of the oscil- 
lation will be of the order of 0.02 cm. The mean free path of a particle at the 
height of 100 km. is much greater than this figure, being of the order of 3 cm. 
Hence the power collected by the receiver will be 1.2X10~* watts per cm.*. 
Since the power must be at least 2X 10~"” watts to be detected by the receiver, 
the radiating volume will have a minimum value of 10" cm.*. It is well to note 
that this is a small volume, and a much larger figure would be necessary if the 
same amount of power were produced as a result of free-free transitions. The 
oscillation of the plasma will be damped by radiation and collisions, but these 
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effects are small for the case we are considering. The radiation emitted during 


sao ot. et OA ; j ; i 
one oscillation is eens 3 while the total energy associated with the oscil- 


, ; Nar e ae ; . 82vq? 
lation is 2x?mv?A®. The ratio of these factors is a 
mc 


7.4 X 10~". Hence the radiation loss per cycle is small. The collisional fre- 
quency at a height of 100 km. may be as high as 2.5 X 10’ if the electron 
density has the large value which we suggest, but this means that collisional 
effects are not too important when oscillations of frequency 3 X 10° are con- 
sidered. This lack of appreciable damping favors the diffusion of the electrons 
through the ionized volume, and the consequent reduction in electron density. 
Collisional damping will have more effect on the emission of longer wave length 


and has the value 


radiation. 


Concluding Remarks 


This discussion of the likely sources of long wave length radiation from the 
polar and nonpolar auroras has suggested that the emission of radiation of 
wave length 10 cm. from the polar aurora is a result of an ionized volume 
executing a short-lived plasma oscillation. It does not seem that other sources 
of radiation of this wave length, such as free—free transitions between electrons 
and ions or neutron particles, can supply sufficient energy to be detected by 
the equipment which has been in use. The ideas presented in this summary 
should be examined more carefully when further data have been accumulated 
on the intensity of the radiation at various wave lengths. It would, perhaps, 
be instructive to examine the radiation for possible polarization effects. The 
ideas presented by Pierce (22), Haeff (11), and others should be investigated 
for possible application to the auroral problem. 
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THE ANGULAR MOMENTUM OF THE ELECTRON 
IN CLASSICAL ELECTRODYNAMICS! 


By F. A. KAEMPFFER 


Abstract 


It is shown that not only the rest-mass of the electron but its spin as well is 
carried by the surrounding electromagnetic field. The ratio of magnetic moment 
and angular momentum (g-factor) of the electron is determined. The law of 
conservation of angular momentum is checked. 


It is a general belief that the angular momentum or ‘‘spin’’ of the electron 
cannot be properly described in terms of the classical theory. The use of the 
word “‘spin’’ indicates that one has in mind a rotating body of finite mass, 
whose angular momentum can be written down in terms of classical mechanics. 
In fact, the physicists of 1900 generally adopted the idea of H. A. Lorentz that 
the electron can be described as a sphere of mass m and surface charge e, 
having a radius ro = | which is defined by the assumption that the rest- 

2mc? 
mass m of the electron is represented by the rest-energy of the surrounding 


electric field E = r according to the equation* 


- [" EMV = me. (1) 


: , 5 en aad , 
If one tries to take into account the spin — x of the electron by assuming that 
T 


the sphere rotates with a certain angular velocity w, so that 


‘ is 
—mrjw =——, 
5 2 Qn 


one is confronted immediately with the paradoxical consequence that the 
outer parts of the electron have to rotate with a velocity 


h he a ; 
ron = en hes 137¢ (3) 
2rmro = 2re? 


9) = 
“max 


far beyond the velocity of light. 


1 Manuscript received in original form August 26, 1949, and, as revised, November 4, 1949. 
Contribution from the Department of Physics, University of British Columbia, Vancouver, 
British Columbia. 

* Actually one has to take into account the stress-energy which is required to hold the surface 
charge e together. This gives an additional factor 4/3 for the total rest-energy of the electron (see 
R. Becker, Elektronentheorie, Leipzig 1933, p. 353), and its radius becomes 

2 é 
= ae (J) 

3 me ) 
A more convenient way of deriving this value is to calculate the total energy of the magnetic field 
created by a slowly moving electron and by equating this energy to the kinetic energy of the electron, 
which leads to the same value (1). 
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It is, however, possible to take into account the angular momentum of the 
electron without assuming any rotational motion of the sphere. This can be 
done by providing the above-mentioned model with a magnetic dipole, which 
is located in its center and which has just the right magnetic moment | 7| = 


eh : - 5 a a : 
—— required by experimental evidence (Stern—Gerlach experiment), so that 


4amc ~ 
M 


: Meee Ne Bie is et 
one has besides the electric field E a magnetic field H = 3 (ru) agg ae 
According to classical electrodynamics this electromagnetic field carries an 
angular momentum 
eB Ree boas Se gece canara 
e=—), [r X[EXH]]dv = -—]| —[r XIrxau]] dV. . (4) 
4ac dre 4dacJdr r§ 
If the z-axis of the co-ordinate system is put into the direction of the magnetic 
dipole , one easily calculates 
ae ee ae eine 2en ig 
6, = 0; 0,= 0; 6, = — — (1 — cos?@) r’dr sin 6dédg@ = ——- (5) 
4c oJo r cro 


By using the exact value (1) for ro one obtains the equation 


To 


| 4 . 
zi = — || (6) 
mc 
which shows that the model under consideration has just the right g-factor 


required by experiment for the electron. In other words, if one puts into (6) 


‘ eh . - ; 
the experimental value —— for the magnetic moment, one gets with 
4amc 
a tk as 
e| = — (7) 
2 2r 


just the right value for the angular momentum of the electron. 


The above equation implies the statement that every dipole surrounded 
by an electric charge carries an angular momentum. In order to show that this 
is not in contradiction with the law of conservation of angular momentum, we 
f_ | which is carried along the radius-vector 

dar? 
r from infinity to 7» through the magnetic field of the dipole. By moving q with 


consider a charge-element g = 


the velocity v through H/ one has to act against the Lorentz force K = 115 HI. 
c 


so that at arrival at 79 one has imposed on gq the momentum 


Pa Kat - [44 x it | a ‘a -|" 2 Ald? = -| fF x ilar. 
J ao dt co Cr To cr 
(8) 


In other words: In order to bring a charge-element g from infinity to a distance 
ro towards a magnetic dipole of moment u one has to put into the system an an- 


gular momentum 
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~ oo q 
6,=-| <b x[r X ull dr. (9) 
ag fF 
; e : s 
Introduction of g = a5 and integration over all angles @ and ¢ leads then 
wre 
immediately to Equation (4), q.e.d. ' 


Thus the model has the advantage that all mechanical properties of the 
electron, including the spin, are represented by its surrounding electromagnetic 
field. It seems to have, however, on first sight the drawback that one departs 
from Ampére’s idea that all magnetism in matter is caused by electric circuits 
and that, if one finds a magnetic dipole in nature, one can always find a circuit, 
which is the ‘‘real’’ cause of it. Of course if one drops the idea that velocities 
beyond c are impossible, one can easily show that the rotating model leads to 
the right magnetic moment, if the surface charge is considered as the circuit 
which causes the magnetic moment according to Biot-Savart’s well known 
formula. In view of the importance of the general principles of relativity theory 
it seems, however, an advantage to drop the idea of Ampére and to postulate 
the existence of magnetic dipoles (which procedure is in fact not essentially 
different from the procedure by which the word “electric charge’ is defined) 
and so to avoid contradiction with first principles of relativity theory. 
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THE DURATION OF THE 16 VOLT LEVELS OF NEON! 









By FRANK A. GRANT 


Abstract 






Electronic apparatus has been developed to determine the rate of decay of the 
16 volt levels of neon after the interruption of an electrical discharge, by measure- 
ment of the absorption of a monochromatic light beam. The half life of the *P2 
metastable level has been determined as a function of pressure and discharge 
tube diameter. The *Po and #P; levels exhibit a pronounced departure from an 
exponential decay. The hypothesis of the imprisonment of resonance radiation 
is consistent with the measured duration of the 'P; and *P, levels. 














I. Introduction 







The most satisfactory method of determining the population of excited atoms 
in a discharge tube is the measurement of the absorption of a light beam of 
appropriate wave length. In neon gas the spectral lines most highly absorbed 
by the 16 volt levels are visible, and their absorption coefficients are known. 
The investigation of the anomalous dispersion of excited gases by Kopfermann 
and Ladenburg has provided data concerning the population of these levels in 
a continuous neon discharge. And because of the ease with which common 
impurities (except argon) may be removed, neon has been used for the investi- 
gation of the duration of excited levels. 












The light absorption method has been used by Dorgelo and Washington, 
Meissner and Graffunder, Anderson, and by Lippert to measure the half life 
of the *P, metastable level of neon, but it was not possible to formulate a 
simple hypothesis consistent with all their results. In particular, the depen- 
dence of half life on the discharge tube diameter was not established, and in 
the results of Meissner and Graffunder, and of Anderson, there appeared to 
be a dependence on the spectral line used to measure the half life. 












The use of the photomultiplier tube and cathode ray tube offers a method 
of obtaining experimental data much more quickly and accurately than was 
possible in previous investigations, and has been employed in the present 
research to attempt to find an explanation of these inconsistencies, and to 
investigate the effect of the imprisonment of resonance radiation. 









II. Absorption of Light—Theory 






Absorption as a Function of Number of Absorbing Atoms 

If a parallel beam of monochromatic radiation traverses a uniform layer of 
absorbing gas, and if a (defined below) is small compared with unity, the 
intensity of the transmitted light is given by the following equation :— 








1 Manuscript received December 18, 1949. 
Contribution from the McLennan Laboratory, University of Toronto, Toronto, Ont. 

2 Holder of National Research Council Fellowships 1946-7, 1947-8. Present address: 

Department of Physics, University of Maryland, College Park, Md. 
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= I,exp (—PWNl), 
length of absorbing path, 
number of absorbing atoms per cubic centimeter, 
absorption coefficient, 
PN, 
constant of proportionality, 
Breadth of emission line, 
Breadth of absorption line ’ 
Intensity of the radiation at a distance / from 
the origin, 
Intensity of the radiation at the origin, 
Intensity of the radiation at any / in the absence 
of absorbing atoms. 


In the case of principal interest, N is an exponential function of time, and the 
half life T is defined by the equation 


: : t ln 2 
N = No exp( - im). (2) 


Combining these two equations and taking logarithms twice, a linear function 


a... log loe(“*) 


of time is obtained :— 
T 
(3) 


Equation (3) is represented by the broken lines of Fig. 1. The slope is inversely 
proportional to the half life 7. The absorption A is defined by Equation (6). 


a) 
1-4 


i 


L0G LOG (7 
he 


T er 
TIME 


Fic. 1. log log I i) plotted as a function of time. Absorption calculated from Equation 


(8). Exponential decay of N assumed. 
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If the quantity a is not small compared with unity, then the absorption is 
given by the following expression :— 


spol oT Gen Rice ae 





(lI 


w is a parameter depending upon the breadth of the absorption line. Using 
the data of Mitchell and Zemansky (10), the corresponding curves for a = 1, 
1.5, and 2, have been plotted in Fig. 1. The departure of these curves from 
linearity is not great. 










Light Beam of Finite Cross Section, Nonuniform Medium 
The absorption A is defined in terms of light flux, so that it represents an 
average absorption across the light beam. 







es ° (5) 
Po 

ae Gi beam JdS “ 
oe beam IodS 


A 


Il 





(6) 





To study the effect of a nonuniform layer of absorbing gas, it is convenient to 
express Equation (4) in the following forms:— 


(Rol) (Rol)? 







(Rol)? 





_= ee (4) 
(1+ a?)* 2!(1 + 2a%)* 3!(1 + 3a?)’ 
r 9 T 2 f 
A, = (PNe) exp( _ tna S oa exp; — = +... (8) 





2 







If PNol varies with radius of the discharge tube, but is small, all terms except 
the first may be neglected, and no error is introduced when the absorption is 
averaged across the light beam. When higher terms cannot be neglected, how- 
ever, integration with respect to the radius results in an expression which is 
not a simple function of the half life 7. Moreover, any attempt to calculate T 
by the foregoing method results in a value that depends on the spectral line 
used, and on the orientation of the light beam with respect to the axis of the 








absorption tube (9). 





Since it is assumed that N = 0 at the wall, it is important to eliminate the 
effect of the light reflected from the wall of the discharge tube. The light beam 
must be restricted to the interior, where the variation of N with radius is 






assumed to be negligible. 
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III. The Decay of Excited Levels After the Cessation of the 
Excitation-Theory 


Transitions between the Triplet Levels 

Two of the triplet levels of neon are metastable (Fig. 2). Under the conditions 
of the present research, the probability that an atom will leave a metastable 
level in a collision is much greater than the probability of a radiative transition. 

The *P» and *P» levels are separated by an energy difference of 0.09 v., and 
kT at 300°K is 0.026. v. Consequently the probability of a transition from a 
metastable level to the *P; level in a collision with a normal atom is not negli- 
gible. The behavior of the *P; level is important in determining the rate of 
decay at the metastable levels. 


O 


Fic. 2. Energy diagram of the 16 volt levels of neon. 


The Imprisonment of Radiation 

Meissner and Graffunder (9) measured the duration of the nonmetastable 
’P, level. At a pressure of 2.15 mm. of mercury they obtained an apparent half 
life of 0.2 milliseconds. In the absence of imprisonment of resonance radiation 
the ‘‘natural”’ lifetime of this level is of the order of 1077 sec. In their work on 
the anomalous dispersion of excited gases (6) Kopfermann and Ladenburg 
found it necessary to assume the imprisonment of radiation to account for the 
concentration of the *P, level. 

Shortley (7) has calculated that the probability of the transition from the 
1P, level to the ground level is about 13 times that from the *P; level. To be 
consistent with the hypothesis of the imprisonment of resonance radiation a 
considerable increase in the lifetime of the singlet level is also to be expected. 


The Decay of the Triplet Levels 
In early studies of metastable levels, attention was directed to the mech- 
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anism by which a metastable atom could revert to the ground level. It was 
assumed that at pressures below 1 mm. of mercury the principal loss of the 
energy of excitation resulted from collisions of metastable atoms with the walls. 
At pressures greater than 2 mm. transitions to the *P; level became increas- 
ingly important. Following this transition, resonance radiation was emitted, 
and reversion to the ground level took place. 
If imprisonment of resonance radiation takes place, however, the emitted 
photon will be absorbed and re-emitted many times, until it 
(a) is absorbed by the wall, 
(b) is absorbed by a neon atom which is in collision with a foreign 
atom, 
(c) is absorbed by a neon atom whose co-ordinates and momenta 
with respect to a neighboring atom are such that a transition 
to another level (most probably *P:2) takes place (4). 


The metastable *P, atom in turn may 
(d) impart its energy directly to the wall, 
(e) be transferred to another level (most probably *P;) by collision 
with a normal neon atom, 
(f) be ionized by collision with another excited atom, 
or (g) excite or ionize a suitable foreign atom (Penning effect). 


Process (f) would cause a hyperbolic decay of the *P:2 level. Since no evidence 
has been obtained for a decay of this type, the process will be ignored. In pure 
neon processes (a), (¢), (d), and (e) determine the rate at which the energy of 
excitation diffuses and is imparted to the wall of the vessel. 


During the excitation process the relative concentrations of the triplet levels 
are determined by the electron temperature (7), and are approximately the 
same as the statistical weights of the levels. After the excitation is interrupted, 
this distribution is altered by the collision processes listed. 


The probability of a transition from one triplet level to another of greater 
energy is dependent on the fact that in only a fraction of the collisions is the 
kinetic energy of the colliding partners sufficiently great. No such restriction 
is effective in the case of the reverse transition. It is therefore reasonable to 
expect a period in which the number of atoms transferred from the *P» and *P; 
levels to the *P» level will greatly exceed the number transferred in the reverse 
direction. After this initial period the rates of decay of the three triplet levels 
should tend to the same value—the half life of the *P2 level. 


The Dependence of T on Pressure 
Zemansky (11) derived an expression for the half life of the *P2 level of the 


form: 1 é B 
Eve 


The constant B was expressed as a function of the atomic weight, temperature, 


+ CP (9) 
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dimensions of the vessel, and ¢? the cross section of the atom for the diffusion 
process. The CP term represented process (e). If imprisonment of radiation 
takes place, the effect of this process is reduced by process (c), which was 
neglected in the derivation of Equation (9). 

Equation (9) will be used empirically until it is found to be inconsistent with 
experimental results. 


Dependence of B and C on Discharge Tube Diameter 

In the investigations of Meissner and Graffunder (9), Anderson (2), and 
Lippert (8), and in preliminary experiments of the present research, Fig. 16, 
the maximum half life at 300°K. was reached at a pressure between 1 and 
2 mm. of mercury, regardless of the discharge tube diameter. This seemed to 
indicate that the constant C varied with diameter in the same way as the 
constant B. The contrary indication of the present research will be discussed 
in a later section. 

Theoretical treatments of the subject were made by Meissner and Graffunder 
(9), Anderson (1), and Zemansky (11). The constant B was found to be depen- 
dent on the tube geometry through the factor 


23.2 


(22. in) 
In a long discharge tube the derived relation indicated that at sufficiently low 
pressures T should be proportional to the square of the discharge tube diameter. 


The experimental results of previous research showed wide disagreement in 
the values of Ba? (Table IV). Much of the present research was undertaken 
to determine whether this disagreement was the result of the inaccuracy of the 
experimental data or of neglect of the imprisonment of resonance radiation. 


IV. The Technique of Measurement 


It has been shown that the absorption of monochromatic radiation is a 
function of a and N [Equation (8)]. If the quantity a is small or known, N (in 





Fic. 3. Diagram illustrating the apparatus used in the determination of the absorption of light 
as a function of time. 
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arbitrary units) can be determined from the absorption of a beam of light of 
the appropriate wave length. 


The emission and absorption discharge tubes, E and A, Fig. 3, both contain 
neon gas. After traversing A, the light from E is focused at the slit of the 
monochromator S. Monochromatic light is incident on the photomultiplier 
tube PMT and the voltage output displayed on the cathode ray tube CRT. 
After interruption of the excitation in A, the absorption is measured as a func- 
tion of time, and used to determine the number of excited atoms present in A. 


ros eee 
| 
i OUTPUT OF PGA 


‘. 


= Nuwax, 
s ‘, 


i N(BLEVEL) 


iv @ -TRANSMITTED LIGHT FLUX - N-Nexpect) 


G2#f#GéwOwO2 #640 
t- TIME MILLISECONDS 


Fic. 4. Diagram illustrating the pulse technique used in the determination of the absorption 
of light as a function of time. 


1. The first requirement of such a system is a voltage pulse (Fig. 41) with a 
rapid decay, to excite the atoms in A. N can then be determined during the 
interval between pulses, when the excitation is interrupted. If the duration of 
the pulse is sufficiently great, an equilibrium is established between the elec- 
trons and excited atoms (7), and conditions in the tube are independent of the 
pulse repetition frequency. For accuracy of measurement, however, it is desir- 
able to make the time interval between these pulses sufficiently great for the 
absorption to reach a negligible value. In this way a zero reference is obtained. 

Kopfermann and Ladenburg (6) have shown that the excitation curve N vs. 
current density exhibits a maximum. If the optimum current density is ex- 
ceeded, the effect of the ‘‘vernichtenden Elektronenstéssen”’ causes a reduction 
of NV. To obtain the maximum concentration of excited atoms regardless of the 
discharge tube diameter or of the pressure, the amplitude of these voltage 
pulses must be adjustable. 


2. The second requirement is a minimum disturbance of the gas in A by the 
light from E. Hence light pulses are used of duration short compared with the 
interval between them (Fig. 4iii). The light flux incident on the photomultiplier 
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tube must produce an output voltage large compared with the “noise’’ output. 
To enable this requirement to be met with the use of a source of low intensity, 
the design of the monochromator is important. Hence the need of a wide slit 
and prisms of high dispersion. 


3. The emission line should exhibit negligible self—reversal. 


4. All forms of broadening of the emission and absorption lines should be 
kept to a minimum. In particular, the value of a should not be large compared 
with unity. The use of pulse excitation reduces the temperature rise in both 
discharge tubes to a negligible value. 


V. Apparatus 
PGA 
Pulse Generator A (Fig. 5) was designed to produce an output voltage of 
up to 2500 v., of duration 0.5 to 1.0 milliseconds, and pulse repetition frequency 
50 to 200 c.p.s. To reduce the effect of stray capacitance, the standard practice 
of using low values of load resistance was followed. 


+ 320 V. VARIABLE TO 2-5 KV. 


DISCHARGE 
TUBE A 


6AG7 807 


Circuit diagram of Pulse Generator A. 


TABLE I 
COMPONENTS OF FIG. 5 


500 picofarads 9, 0.05 microfarads 17. 2 kilohms 
8 microfarads 10. 4.7 kilohms 18. 0.05 microfarads 


1 kilohm it. 5-6 19. 470 kilohms 
1.3 5 12. 0.1 microfarads 20. 10 os 
2. ‘s 13. 22 kilohms 21. 100 ohms 
5.6 - . 1000 picofarads 22. 10 kilohms 
500 ‘i 5. 800 kilohms 

9 ae 


8 microfarads | fae 


PGE 
Pulse Generator E (Fig. 6) was designed to produce an output voltage up 
to 1500 v., of duration 0.02 milliseconds, and a pulse repetition frequency 1000 


C.p.s. 
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UP TO /-SKK 


DISCHARGE 


807 807 
Fic. 6. Circuit diagram of Pulse Generator E. 
TABLE II 
COMPONENTS OF FIG. 6 





kilohms 
microfarads 
kilohm 
microf arads 


5 kilohms | 9. 15 kilohms 


; 4 
es . 100 - 
2 -( two 10k in 
picofarads 3.) parallel 
kilohms 4. 0.1 microfarads 
zp . 100 ohms 
picofarads . 470 kilohms 


© 
— 
“1S 0 St br 


kilohms 


QO MIR Oe Co br 





Discharge Tube A 

The absorption discharge tube A (Fig. 7) was circular in cross section with 
optically plane Pyrex windows at both ends. The electrodes were circular and 
made of 20 mil tungsten wire concentric with the walls of the discharge tube. 


Fic. 7. The absorption discharge tube A. 


Three tubes, of 30 mm., 46 mm., and 81 mm. inside diameter were used. The 
power dissipated in the tubes was so small that no increase in temperature 
could be detected by the hand. 


Discharge Tube E 

The light was emitted from a capillary tube 4 mm. in diameter and 4 cm. in 
length. The discharge tube was filled with neon at a pressure of 3 mm. The 
power dissipated in the emission tube was less than 0.5 w. 


Photomultiplier Tube 
A 931-A photomultiplier tube with accelerating voltages of 105 v. was used. 
The load resistor was 19,000 ohms. 
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Cathode Ray Tube 
This was a Dumont Type 208-B with a blue trace. 


Monochromator 

To enable a wide slit to be used, three carbon disulphide prisms were used in 
the monochromator. In this way a satisfactory signal—noise ratio in the output 
of the photomultiplier tube was obtained with the lowest possible intensity of 
the light source. 


Baffles 

Paper baffles, B,; and By (Fig. 3), were attached to both ends of the absorption 
tube. Circular apertures 25 mm. in diameter concentric with the axis of the 
tube determined the beam diameter. To eliminate the effect of light reflected 
from the wall of the absorption tube, an image of the first aperture was focused 
at the slit of the monochromator, and a third baffle was inserted to allow only 
the light forming this image to enter the slit. 


Gaseous Purity 

The neon was purified first by a discharge in a potassium pool discharge tube 
permanently connected to the absorption tube, and, second, by the heating to 
incandescence of the tungsten electrodes in A for several hours. The second 
method appeared to be the more effective in removing hydrogen, but the 
potassium was very effective in the removal of mercury vapor. The gas was 
found to contain a trace of argon, which of course was not removed by these 


methods. 


GROUND GLASS 


at JOINT 


TO MANOMETERS 
AND PUMPS 


8imm. 1D. 


LIQUID AIR 
TRAP 


POTASSIUM POOL 
DISCHARGE TUBE 


Fic. 8. Vacuum system. 


Since stopcocks caused a serious deterioration of the purity, all purification 
was accomplished with the system sealed off. To change the pressure, the 
discharge tube A was connected to the manometers and pump by breaking a 
glass seal S (Fig. 8). By the use of a magnet, it was possible to raise the steel 
pellet P and allow it to fall, breaking the seal S. To avoid contamination of the 


neon, the pellet was enclosed in glass. 
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Manometers 






The pressure was measured with a McLeod gauge and an oil manometer, 
using Octoil-S. 










Synchronism 

For convenience of time measurement, both pulse generators and cathode 
ray tube were synchronized with an audio frequency generater A FG, Fig. 3. 
In all cases the pulse repetition frequency of PGE was about 1000 c.p.s. so that 
the time interval between pulses was about 1.0 millisecond. In some cases 
both generators were synchronized with a 2000 c.p.s. voltage from AFG. The 
period of PGE was adjusted to be an even multiple (two) of 0.5 millisecond, 
and the period of PGA to an odd multiple, making a cathode ray tube pattern 
with an apparent interval between successive PGE pulses of 0.5 millisecond. 














Photography 

The noise level of the system was sufficiently high to make direct measure- 
ment of the cathode ray tube trace impracticable. This effect was eliminated 
by making a 30 sec. exposure when photographing the trace. 

The effect on the curvature of the cathode ray tube screen was reduced by 
limiting the amplitude of the image to 2} in. (one half the screen diameter), 
and by placing the reference line (zero absorption) within an inch of the center 
of the screen. In addition to this, absorption of less than 5% was ignored when 











making measurements. 







VI. Summary of the Phenomena Investigated 






Imprisonment of Radiation.—Because of the relatively small absorption, a 
long discharge tube (90 cm. long, 46 mm. inside diam.) was used to investigate 
the duration of the 'P; and #P, levels. 











The Decay of the Triplet Levels —Because of the convenient path length 
(90 cm.), the investigation of the decay of the *P levels was made with the 
same discharge tube. Reproducibility was checked by making two independent 
observations of the decay of the *P» and *P, levels; the results of these investi- 
gations are plotted in Figs. 12, 13, and 14. 







The *P. Level—Various Spectral Lines.—The dependence of the half life of 
the *P» level on the spectral line used was investigated at the same pressure 
and using the same discharge tube as in the foregoing investigations. The 
results are plotted in Fig. 15. 







8P. Level-The Constants B and C.—The dependence of the constants B and 
C of Equation (9) on the discharge tube diameter was determined using tubes 
of inside diameter 30, 46, and 81 mm., all 70 cm. in length. Two sets of obser- 
vations were made, and are plotted in Figs. 17 to 20. 







In all cases the interval between PGE pulses was 1 millisecond. 
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VII. The Treatment of Experimental Data 


In order to determine the half life of the *#P, metastable level, it is necessary 
to make an assumption of the value of a (see Fig. 1). Because of the negligible 
temperature rise in both tubes, the Doppler breadth of both emission and 
absorption lines was assumed to be the same, and a taken as unity. The dif- 
ference between the slopes of the curves a = 0 anda = 1, between the ordinates 
1.6 and 1.0 is 6%. Whenever possible the portion of the curve lying between 
these limits (Fig. 9 and Fig. 10) was used for the determination of the slope, 
and the value of the half life so obtained was reduced by 6%. 


8/ MM. DIAM. 
T = 34, MSECS, 


46 MM. DIAM, 
T= /-67 MSECS. 


\ 30 MM. DIAM, 
‘oe 0-96 MSECS. 


PRESSURE = 0-48 MM. 
L=70CM. 


3 4 Ss 
MILLISECONDS 


Fic. 9. The decay of the quantity log log (=) as a function of time. Wave length 5945A, 
8P» level. g 


It was then attempted to represent the half life of the *P2 level by the em- 
dirical relation :— 
' a 


ae ae GR: 9 
: a (9) 


Using the method of Mitchell and Zemansky (10), the quantity z was plot- 


ted against P? in Figs. 17 and 18. The value of B is the intercept made by the 
straight line on the vertical axis, and the value of C is equal to the slope. The 
curves of Figs. 19 and 20 are plotted using the values of B and C obtained from 
Figs. 17 and 18. 

Measurements of the absorption were made from the photographs of the 
cathode ray tube screen by placing the photographic plate on a piece of squared 
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paper. In these photographs, of which Fig. 11 is an example, zero absorption 
corresponds to the asymptote of the envelope of the figure as time increases 
without limit. 


= 
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Fic. 10. The decay of the quantity log log (; tS as a function of time. Wave length 59454, 


ae 


3P, level. 
VIII. Experimental Results and Discussion 


Imprisonment of Radiation 

To investigate the effect of imprisonment of resonance radiation in length- 
ening the life of the *P2 level, measurements of the duration of the 'P; and *P, 
were made. The half life of the 1P; level was found to be 0.8 milliseconds (Fig. 
11). The decay of the *P; level is not exponential (Fig. 12), but the absorption 
does not decrease to 10% until 3 milliseconds after the interruption of the 
excitation. 

Both these measurements are consistent with the hypothesis of imprisonment 
of resonance radiation. 


The Decay of the Triplet Levels 

In Figs 12 and 13, the first three milliseconds constitutes the period of tran- 
sition discussed previously. A pronounced departure from an exponential decay 
is exhibited by the *P» and *P; levels. For comparison, the rate of decay of the 
’P, level is represented by the broken line. Measurements of the rate of decay 
of the *P> and *P; levels were made by Dorgelo (3) and by Meissner and 
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Fic. 11. 'P, level. Photograph of the oscillogram showing the absorption as a function of time. 
Spectral lines 6678 and 67174, Discharge tube diameter, 46 mm.; length, 90 cm. Pressure, 1.3 
mm. Intervals, 1/2.06 millisecond. . 
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Fic. 12. The decay of the quantity log log (; I i) as a function of time. Wave length, 6096A. 
3P, level. 


Graffunder (9). The accuracy of measurement, however, was insufficient to 
indicate a departure from an exponential decay. 

In the decay of the *P2 level, the period of transition is less evident (Fig. 14). 
The *P, Level-Various Spectral Lines 


The results of Meissner and Graffunder (9) and of Anderson (2) indicated 
that the apparent half life of the *P, level depended on the spectral line used. 
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A possible cause of this effect is the variation of the quantity PN/ with radius 
equation (8)]. A similar effect would be produced if stray light from the 
emission tube were allowed to enter the photomultiplier tube. Precautions were 


2 
~ 


6266 (°P,) 
DIAMETER, 46 MM. 
LENGTH, | 90 CM, 
PRESSURE, /-3 MM, 


0-2 0-3. 040-506 07 


ABSORPTION 


0 
MILLISECONDS 


Fic. 13. The decay of the quantity log log (=) as a function of time. Wave length, 6266A. 


3Po level. 


taken to eliminate these effects by the introduction of baffles and by reducing 
to a negligible intensity the stray light entering the photomultiplier tube. In 
spite of these precautions, there still remained the following effect to be 
explained. 


6 8 /2 14 
MILLISECONDS 


Fic. 14. The decay of the quantity log log a as a function of time. Wave length 5O45A, 
a oa g 
3P » level. 


.) is plotted as a function of time for 


In Fig. 15 the quantity log loe(- 


the wave lengths 6402, 5945, 6143, 6334, and 5882A. The absorption of all 
these spectral lines is dependent on the *P» level. The plotted curves have been 
calculated in the same manner as those of Fig. 1. In every case a half life of 
2.35 milliseconds has been assumed. The origins of the time scales have been 
chosen to give the best fit to the experimental points. If the same half life is 
assumed for all curves, a good fit can only be obtained by assuming different 
values of a. 
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Fic. 15. *%P: level. Curves calculated as explained in the text. 


The *P, Level—The Constant C 

No satisfactory interpretation of the constant C has been made (5). Since 
the neon used in the present research was found to contain a trace of argon, it 
is possible that this impurity is responsible for the CP term. The slope of —1 
(Fig. 16) indicates that the effect is proportional to the first power of the 
pressure. 
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Fic. 16. *P2 level. log T plotted as a function of log P, illustrating the predominance of the 
CP term of Equation (9) at higher pressures. 


It was stated in a foregoing section that the results of previous investigations 
seemed to indicate that the maximum half life occurred at a pressure between 
1 and 2 mm. of mercury regardless of the diameter of the discharge tube. To 
investigate this point, it was important that discharge tubes of different 
diameters, containing identical gas, should be used. Consequently, the three 
tubes used in this research were permanently connected by 15 mm. tubing, and 
mounted so that one tube could be removed and replaced by another without 
disturbing the remainder of the system. 
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In Figs 17-20, the three readings corresponding to any given pressure were 
obtained within an interval of three minutes. The readings were taken in the 
sequence 81-46-30, 30-46-81, 81-46-30, etc. 


There is no close agreement between the values of C (Table V) obtained. 
However, the maxima of these curves are reached at different pressures, and 
the conclusion is reached that B and C do not vary with diameter in the same 
way. 


Dependence of B on Diameter 

Figs. 17 and 18 illustrate the manner in which the constant B was determined 
from the experimental data. It is assumed that the end effect is negligible in a 
discharge tube 70 cm. long. To determine if B is inversely proportional to the 
square of the discharge tube diameter in a long tube, values of Ba? are tabulated 
in Table III. 


Fic. 17. %P. level. The determination of the constants B and C of Equation (9). 
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TABLE III 
EXPERIMENTAL VALUES OF CONSTANT B OF EQUATION (9 
3P. Level 
Wave length 5945A 
L = 70 cm. 


Inside 
diameter, 
a, 

mm. 


30 500 45 
30 460 4] 
46 225 48 
46 240 51 
81 90 | 59 
81 75 49 


Probable 
B Ba? X 10-4 error, 


Res ae 


These results indicate a slight increase of Ba? with increase in the discharge 
tube diameter. 


OS pe +0 rs 
Fic. 18. *P2 level. The determination of the constants B and C of Equation (9). 
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FiG. 19. *Ps level. log T plotted as a function of log P. Curves calculated using Equation (9) 
and the values of B and C obtained from Fig. 17. L = 70 cm. 
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TABLE IV 
EXPERIMENTAL VALUES OF CONSTANT B 





Inside 
diameter, 
a, mm. 


Meissner and Graffunder reas 36 
AGGOOO S65). .ao w asle att 10 
Preliminary measurement, present research| 46 


TABLE V 
EXPERIMENTAL VALUES OF CONSTANT C OF EQUATION (9) 
3P2 Level 
Wave length 5945A 
L = 70 cm. 
Inside | 
diameter | 


30 
30 
46 
46 
$l 
81 


3 
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Fic. 20. *Pz level. log T plotted as a function of log P. Curves calculated using Equation (9) 
and the values of B and C obtained from Fig. 18. L = 70cm. 
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